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INTRODUCTION 


» 

'^This  volume  describes  Che  systems  content  of  the  ROSCOE  code. 


The 


systems  model,  as  defined  here,  consists  of  these  major  blocks:  attack 

generation,  object  observables,  sensor  platforms,  track  filter,  radar 
event  logic,  radar  signal  processing,  and  radar  propagation. 


I 


The  objective  in  building  the  systems  model  was  to  provide  a very 
generalized  set  of  routines  which  could  be  used  by  a variety  of  users  for 
a variety  of  applications.  At;  a result,  the  user  can  select  from  among  a 
number  of  options  to  simulate  a particular  engagement.  For  example,  any 
number  of  different  object  types,  radar  types,  and  weapon  types  can  be 
specified  in  a single  run.  In  addition,  launch  and  impact  points  can  be 
mixed  to  generate  complex  attacks.  1^- 

Section  2 of  this  volume  briefly  describes  the  ROSCOE  models  of  at- 
tack generation  and  object  motion,  object  radar  observables,  sensor  plat- 
forms, and  radar  track  filter.  References  are  given  for  fuller  descriptions. 

Section  3 describes  the  radar  logic  in  terms  of  the  radar  "events": 
search/detect/verify,  track  initiation,  track,  and  discrimination.  Assump- 
tions are  explicitly  identified;  procedures  are  described  in  detail;  and 
annotated  flow  charts  are  presented. 

Section  4 describes  the  ROSCOE  models  of  the  signal  processing  that 
is  applied  to  each  radar  pulse  (radar  look  event).  These  models  account 
for  the  effects  of  the  environment  (by  calling  the  propagation  routines, 

Sec.  5)  and  of  radar  characteristics  such  as  monopulse  processing,  range 
gating,  and  signal  design  (ambiguity  function). 

Section  5 describes  the  radar  propagation  models.  These  models  link 
the  systems  part  of  ROSCOE  with  the  environment  p-rt  by  calculating  the 
effects  on  radar  signals  of  absorption,  fireball  noise,  refraction,  clutter, 
and  multipath. 
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} 

i 

Appendixes  supply  the  background  for  some  of  the  signal  processing 
algorithms. 
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MISCELLANEOUS  MODELS 


2 . 1 ATTACK  GENERATOR  AND  OBJECT  MOTION 

The  attack  generator  for  ROSCOE  is  taken  from  General  Research  Cor- 
poration's BAG  14  code.*  It  consists  of  the  following  set  of  programs: 


ATKGEN 

LAUNCH 

PHIMP 

Attack  generation  program 

Creates  objects  and  object  position  arrays 

Creates  burst  arrays  at  impact 

The  programs  use  the  following  subroutines  to  generate  the  object  trajec- 
tories ard  find  specific  points  on  them: 

0RB2 

Determines  the  elements  of  a Keplerian  orbit 
* 

joining  two  points 

0RB1 

Determines  the  orbital  elements  corresponding 
* 

to  a given  state 

ORBP 

Determines  the  state  vector  of  a point  on  a 

A 

Keplerian  orbit 

BLLSTIC 

Generates  an  endo-atmospheric  trajectory  with 
slowdown 

DR0RB2 

De  -orbit  routine  which  modifies  an  orbit  to 
meet  a specified  impact  point 

GENORB 

Corrects  an  orbit  to  meet  a specified  Impact  time 

ADVANCE 

Advances  the  state  and  covariance  of  an  object 
to  the  current  time 

KUTTA 

Trajectory  integration  routine 

PREDLOC 

Returns  a predicted  object  location 

WHERE 

Finds  object  state  at  a specified  time 

* 

These  routines 

are  part  of  the  General  Research  Corporation  TRAID  system, 

11 
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Figure  2.1  Is  a flow  diagram  of  the  main  program  in  this  module 
(ATKGEN) . The  Input  consists  of  a set  of  launch  points  (dataset  LP)  and 
a set  of  target  points  (dataset  TG)  specified  by  the  user.  These  datasets 
may  have  repeated  entries  and  may  be  In  quite  arbitrary  order.  The  attack 
generator  takes  each  target  point  in  order.  Associated  with  the  target 
point  as  an  input  is  the  number  of  reentry  objects  which  the  target  point 
is  to  receive.  This  number  of  objects  is  assigned  to  the  target  point 
one  at  a time  by  a sequential  pass  through  the  list  of  launch  points.  If 
the  end  of  the  list  is  reached  before  all  the  target  points  are  satisfied, 
the  program  cycles  back  to  the  front  of  the  list  of  launch  points  and  con- 
tinues. Thus,  in  a simple  attack  with  a single  launch  point,  one  entry 
on  the  list  is  sufficient.  It  should  be  noted  parenthetically  that  object 
type  is  peculiar  to  launch  point;  thus  to  mix  RVs  and  decoys  requires 
mixing  launch  points  on  the  launch  point  list. 

The  timing  of  the  attack  and  the  type  of  trajectory  are  specified 
with  the  target  point  (dataset  TG) . These  data  Include  the  desired  arrival 
time  for  the  first  booster,  the  desired  time  separation  between  objects, 
the  sigma  of  arrival  time  errors,  and  the  CEP  of  position  errors.  An 
alternative  form  of  timing  specification  allows  the  user  to  specify  a 
fixed  time  Interval  inside  which  the  reentry  objects  assigned  to  a 
particular  target  point  arrive  in  random  order. 

2 

The  trajectory  is  specified  as  in  TRAID  and  is  calculated  with  drag 
on  a round,  rotating  earth.  For  each  reentry  object,  the  attack  generator 
generates  a launch  event  at  the  time  at  which  the  Keplerian  position  of 
the  orbit  leaves  the  launch  point.  The  storage  for  each  object  will  be 
set  aside  when  the  launch  event  is  processed — all  that  is  saved  during 
the  attack  generation  is  the  computed  orbital  elements,  which  are  appended 
to  the  launch  event. 

The  launch  event  creates  an  object  and  names  it.  (The  object 
LP-5  7 is  the  seventh  launch  from  the  launch  point  called  LP-5.)  Object 
and  object  position  arrays  are  created.  Orbital  element  and  reentry 
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Figure  2.1.  Program  ATKGEN 
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position  state  sectors  are  constructed,  as  Is  a summary  output  array  for 
that  object.  Future  events  associated  with  the  object  are  constructed 
and  placed  on  the  event  list.  These  are 

• The  Impact  event  (Program  PHIMP) 

• A potential  acquisition  event  (Program  RADAR)  for  each 

radar  In  the  defense  system 

The  Impact  event  Is  created  at  the  same  time  the  object  Is  created 
and  occurs  at  the  moment  that  the  object  trajectory  passes  through  Its 
designated  target  point.  In  particular,  this  means  that  If  an  alrburst 
(or  a precursor)  Is  desired,  the  target  point  should  be  specified  at  a 
non-zero  altitude. 

2.2  OBJECT  OBSERVABLES 

This  section  describes  the  following  object  observables  models: 

1.  Ballistic  coefficient 

2.  Tumbling  dynamics 

3.  Radar  cross  section 

3 

The  models  are  taken  from  General  Research  Corporation’s  SETS  model. 

In  the  following  sections,  flow  diagrams  are  shown  for  each  routine 
with  explanatory  notes  to  define  the  variables  and  describe  the  computations. 

2.2.1  Ballistic  Coefficient 

The  ballistic  coefficient  (beta)  Is  found  by  calling  Function  BETAGT 
(Fig.  2.2)  with  the  object  state  and  a dataset  pointer  (DSP)  to  the  bal- 
listic coefficient  model  type  desired.  There  are  three  model  types:  type  1 
uses  a constant  value  of  beta  for  all  altitudes;  type  2 uses  an  N-point 
altitude-beta  table,  resulting  in  a piecewlse-llnear  relation  between  beta 
and  altitude;  and  type  3 relates  the  beta  history  of  an  object  to  the 
dynamics  of  a theoretical  conical  vehicle  of  constant  mass  and  shape. 


Ik 
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Figure  2.2.  Function  BETAGT 
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In  type  2,  the  altitude-beta  pairs  are  given  in  order  of  decreasing 
altitude.  For  object  altitudes  above  the  first  in  the  list,  the  beta 
value  for  the  highest  altitude  is  used.  For  object  altitudes  below  the 
last  in  the  list,  the  last  beta  value  is  used. 

In  type  3,  the  quantities  shown  in  Fig.  2.2  are  defined  as  follows: 
m ■ mass  of  vehicle,  lb 

2 

A ■ aerodynamic  reference  area,  ft 
cx  ■ axial  drag  coefficient 
0c  ■ theoretical  cone  angle,  deg 
M ■ raach  number 
c - local  speed  of  sound,  m/s 
v ■ current  vehicle  velocity,  m/s 

At  altitudes  above  150  km,  type  3 uses  an  input  minimum  beta,  with  no 
further  computation  required. 

The  ballistic  coefficient  computed  by  model  type  3 may  be  randomized 
by  adding  a number  drawn  from  a Gaussian  distribution  G(0,  SIGB)  with 
zero  mean  and  an  input  standard  deviation  SIGB. 

2.2.2  Tumbling 

The  tumbling  model  (Subroutine  TUMBLR,  Fig.  2.3)  is  called  with  the 
object  state  vector  and  a DSP  word  denoting  the  model  type,  and  returns 
a vector  (AVEC)  describing  the  motion  of  the  object  about  its  point  mass 
trajectory. 

There  are  three  tumbling  models  available.  They  all  assume  that 
the  object  is  a solid  of  revolution  with  one  axis  of  rotational  symmetry 
(body  axis) . It  is  assumed  that  the  cone  swept  out  by  the  tumbling  body 
axis  differs  negligibly  from  a plane  (tumble  plane).  Type  1 keeps  the 
object  oriented  along  its  current  velocity  vector;  therefore,  no  data  is 
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Figure  2.3.  Subroutine  TUMBLR 
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required  for  Chls  model.  Type  2 orients  the  body  along  its  initial  reentry 
velocity  vector  (computed  at  100  km,  or  at  the  object's  current  altitude 
if  the  initial  call  occurs  below  100  km).  Type  3 is  a stochastic  model 
which  allows  an  object  to  tumble  at  a given  rate  until  its  body  axis  sta- 
bilizes along  the  velocity  vector  at  a given  (input)  stabilization  altitude. 
The  initial  body  orientation  and  tumble  plane  are  derived  as  follows: 

1.  Three  random  components  are  chosen  from  a uniform  random 
distribution. 

2.  A three-vector  formed  by  normalizing  these,  components 
defines  a tumble  axis  (TAXIS) . 

3.  A tumble  plane  normal  to  this  axis  is  formed  and  a random 
vector  in  the  plane  chosen  by  selecting  two  more  random 
numbers  from  the  uniform  distribution.  This  latter  vector 
is  called  the  body  axis  (AVEC). 

A.  After  initialization,  the  body  axis  is  updated  by  using 
the  input  tumble  rate  and  the  definition  of  the  tumble 
plane  (hQ,uo). 

2.2.3  Radar  Cross  Section 

Four  models  are  available  for  modeling  an  object's  RCS  (Function 
RCSMODL,  Fig.  2. A).  Type  1 assumes  the  RCS  is  constant  and  requires  no 
other  input  data.  The  other  three  models  compute  RCS  as  a function  of 
the  aspect  angle  (defined  as  the  angle  between  the  radar  line-of-slght 
and  the  current  orientation  of  the  body  axis  as  given  by  the  tumbling 
model) . 

In  the  type  2 model,  the  RCS  is  found  by  interpolation  of  an  RCS 
vs  aspect  angle  table.  Type  3 is  intended  for  tank-like  objects  (cylinders) , 
while  type  A is  for  RVs  and  decoys. 
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RCS  = CONSTANT 


REQUEST: 

OBJECT  DATASET 


CALL  TUMBLE 
SUBROUTINE 
FOR  BAX  IS 


COMPUTE  ASPECT 

ANGLE : 

rRR  = rRV  ‘ 

rRAD 

-1 

Frr-BAXIS 

« * cos 
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INTERPOLATION 
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Figure  2.4.  Function  RCSMODL 
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In  Fig.  2.4,  the  variables  used  in  type  3 are  defined  as  follows: 
a ■ radius  of  tank,  m 
L ■ length  of  tank,  m 
I - aspect  angle,  deg 
k - 2it/\ 

\ - wavelength,  m 

and  the  variables  used  in  type  4 are  defined  as  follows: 

w - spin  rate,  rad/s 

a - number  of  RCS  peaks  per  revolution 
A,B  ■ constants  for  vehicle  in  question,  as  used  in  the 
equation  shown  in  Fig.  2.4 

2.3  SENSOR  PLATFORMS 

Three  sensor  platform  models  are  currently  available  in  ROSCOE:  fixed, 

orbital,  and  circular-orbit.  They  are  specified  by  placing  the  proper 
dataset  (PI,  P2,  or  P3,  respectively)  in  the  input  deck.  The  platform 
(radar)  location  is  found  by  calling  subroutine  PLTFRM  with  a pointer  to 
the  dataset  containing  the  platform  position  information.  The  routine 
returns  the  position  through  the  calling  sequence. 

The  "fixed"  platform  is  simply  described  by  a position  vector.  The 
"orbital"  platform  is  described  by  a TRAID  ten-vector  giving  its  orbital 
elements.  The  "circular-orbit"  platform  is  described  by  a simplified  set 
of  orbital  elements  defining  a circular  orbit. 

2.4  TRACK  FILTER 

The  GRC  BAG  14  tracking  filter  has  been  suitably  modified  for  use 
in  ROSCOE.  It  is  a fully  coupled  seven-state  Kalman  filter  with  an  option 
for  exponential  memory  decay. 
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The  filter  module  is  made  up  of  five  subroutines: 


FILTER 


CHOLSKI 


JACOBIAN 


XFORM 


KALMAN 


Initializes  covariance  and  updates  filter 
predictions  (Fig.  2.5) 

Inverts  an  N * N matrix 

Sets  up  Jacobian  matrix  for  filter  routine 
Inverts  Gamma  matrix  to  Covariance  and  vice  versa 
Kalman  filter 


The  formulation  of  the  tracking  filter,  expressed  in  radar  face 

coordinates,  is  as  follows.  If  we  define  the  three-dimensional  vector 
-► 

x by 


x^  ■ r sin  a 


x2  - r sin  B 

2 2 1/2 
x3  ■ r(l  - sin  a - sin  8) 


where  r,  a,  6 are  range  and  radar  angular  coordinates,  respectively, 
then  the  seven-dimensional  state  vector  estimated  by  -he  filter  is  given  by 


X “ (x,xA,A) 


where  A is  the  time  between  pulses  for  a given  track  file  and  A is 

the  natural  logarithm  of  the  ballistic  coefficient.  The  equations  of 

-► 

motion  giving  the  time  history  of  X are  assumed  to  be 


-P8(e) 


R - x + R 


A - u(t) 


-A  . 


+ u x x + u x (ojx  Rfi) 


(2.1) 


21 


Figure  2.5.  Subroutine  FILTER 
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where 


Rg  - vector  from  earth’s  center  to  radar 

*♦  . 

u)  ■ earth  s angular  rotation  vector 

u - stationary,  uncorrelated  random  variable  with  mean 
zero  and  standard  deviation  o 

p 

G - universal  gravitational  constant 
o ■ atmospheric  density 
M » mass  of  earth 

This  formulation  differs  from  the  standard  three-degree-of-freedom  equation 
for  reentry  motion  in  that  no  lift  is  considered,  and  the  centripetal 
acceleration  term  is  taken  as  a constant. 

The  ROSCOE  tracking  filter  is  basically  an  extended  Kalman  filter 
and,  as  such,  gives  a running  estimate  of  the  error  covariance  matrix  C 
associated  with  X . The  filter  can  be  separated  into  three  distinct 
operations — initiation,  prediction,  and  correction — of  which  only  the 
last  two  are  repetitive.  To  distinguish  the  estimates  of  the  state  vector 
and  covariance  matrix  resulting  from  the  prediction  and  correction  steps, 
they  are  given  the  subscript  p or  c , respectively. 

The  filter  can  be  initiated  either  through  the  track  initiation  logic 

or  through  handover  from  another  radar.  In  both  cases,  it  is  assumed 

that  the  first  six  components  of  £ are  given  at  the  time  of  the  first 

P 

tracking  pulse,  along  with  some  estimate  of  variance  and  covariance  between 
them.  The  initial  estimate  of  A is  currently  in(100  lb/ft^)  ■ 4.605  , 
and  the  initial  variance  is  equal  to  4.  The  covariances  between  A and 
the  remaining  components  are  initiated  as  zero. 

Having  an  estimate  of  X and  C at  some  time  t , the  filter  must 
be  able  to  predict  the  best  estimate  of  these  quantities  at  some  future 
time,  usually  to  determine  the  radar  gating  and  beam  pointing  for  the  next 
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pulse.  The  prediction  step  is  always  carried  out  over  one  cycle  time 
A . The  subvectors  of  X are  predicted  by  the  equations 

* “ . 

x (t+A)  - X (t)  + X (t) A + * (t)A2/2 
P C c c 

• « •• 
xp(t+A)  - xc(t)  + xc(t)A 

Ap(t+A)  - Ac(t) 

where  xc  is  evaluated  from  Eq.  2.1.  For  the  first  prediction,  the 
initial  values  are  used  in  Eq.  2.2.  The  covariance  matrix  is  predicted 
through  the  use  of  an  approximate  transition  matrix  $(t+A,t)  • The  usual 
approximation  starts  with  the  equations  of  motion  written  in  the  form 

i - ?«) 

Then  the  equations  are  linearized  about  some  nominal  value  of  the  state 
vector 


3? 


6X  - (X  - X ) * — (X  )6X 
nom  nom 


where  X is  usually  taken  to  be  the  best  estimate  X . The  linearized 
nom  ' c 

equations  then  furnish  a transition  matrix  through  the  matrix  formula  (to 
first  order  in  A): 


t+A 

<J>(t+A,t)  ■ exp  J"  F(t')  dt' 

t 


I + F(t)A 


/ 


In  the  ROSCOE  tracking  filter,  the  further  approximation  is  made  in  F 
that 


so  that  only  the  partlals  with  respect  to  \ contribute  to  the  last  four 
rows  of  p . Once  P Is  evaluated,  C Is  predicted  through  the  formula 


Cp(t+A)  - P(t  A,t)Cc(t)^(t+A,t)  + Q 


(2.3) 


where  Q is  a matrix  whose  only  non-zero  element  Is  In  the  last  row  and 

2 * 

column  and  Is  equal  to  o . 

V 

When  a measurement  Is  made,  It  Is  incorporated  Into  the  estimate 
through  the  correction  step.  For  numerical  reasons,  a pseudo-observation 
vector  Y - (x,y,r)  Is  computed  from  the  measured  values  of  r,  sin  a, 
and  sin  3 • The  single-pulse  variance  estimates  for  r,  sin  a,  and 
sin  3 are  computed  from  the  measured  signal-to-nolse  ratio  and  the  pulse 
type.  If  the  matrix  p is  defined  as 


"sin  a 

r 

0- 

p ■ 

sin  3 

0 

r 

r 

. 1 

0 

0- 

then  the  estimated  error  covariance  matrix  for  Y Is  computed  through  the 
formula 


R ■ v 


0 

2 

sina 

0 


0 
0 

.2 

sinSJ 


The  3*7  matrix  M Is  defined  by 


M 


ij 


3Yt(X) 

9x. 


3 - X 


* 2 

Currently  is  taken  as  0.01A  for  altitudes  less  than  250  kft  and 

zero  otherwise.  The  Q matrix  keeps  the  filter  memory  finite  within 
the  atmosphere. 
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Then  the  correction  formulas  are  given  by 


Xc(t)  - Xp(t)  - W[Y(t)  - Mxp (t) ] 


C (t)  - C (t)  - WMC(t) 
c p p 

W - C MT[MC  MT  + R] 
c c 


If  for  some  reason  no  measurement  Is  available,  the  filter  sets  X 

c 

and  C£  * Cp  in  preparation  for  the  next  prediction  step. 


-f 

X 

P 


The  use  of  the  Q matrix  in  Eq.  2 3 keeps  the  filter  memory  finite 
within  the  atmosphere  and  is  employed  only  below  250  kft  altitude.  Another 
useful  technique  for  decreasing  memory  length  in  tracking  filters  is 
through  the  multiplication  of  the  covariance  matrix  by  a "fading  factor" 

c’  - rc  r* 

p p 

where  T is  a diagonal  matrix  having  non-zero  elements  greater  than  or 
equal  to  one.  In  the  ROSCOE  filter,  a variation  of  this  approach  is  used 
wherein  the  fading  factor  is  applied  only  to  the  diagonal  elements  of 
the  covariance  matrix.  The  elements  along  the  diagonal  of  r are  given 
by  (1,1,1,y,y,y,1)  where  y ■ e T and  t is  a constant  which  can  be 
varied  with  altitude. 
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3 


RADAR  EVENT  LOGIC 


The  ROSCOE  radar  model  may  be  divided  Into  two  subsets:  radar  event 

logic  and  radar  signal  processing. 

The  radar  event  logic  is  responsible  for  setting  up  and  simulating 
the  search,  verification,  track  initiation,  track,  and  discrimination 
events.  Each  event  must  contain  the  logic  for  selecting  the  appropriate 
signals,  transmitting  them  to  the  appropriate  points,  calling  for  the 
correct  data  processing  of  any  return,  deciding  (based  on  the  latest 
return)  what  further  events  are  necessary,  and  finally  setting  up  future 
calls  for  these  events  at  the  appropriate  time. 

This  type  of  event  logic  is  essentially  independent  of  how  the  radar 
signal  processing  is  simulated.  The  latter  therefore  constitutes  an 
independent  segment  of  the  radar  model,  and  will  be  discussed  separately 
in  Sec.  4. 

The  event  logic  is  described  in  the  subsections  which  follow.  Each 
subsection  describes  the  assumptions  which  underlie  the  logic,  the  process 
which  is  being  modeled,  and  any  limitations  of  the  model.  Flow  charts 
are  presented  and  discussed. 

3.1  SEARCH/DETECTION/VERIFICATION  EVENTS 

3.1.1  Assumptions 

During  these  events,  the  system  is  assumed  to  operate  without 
monopulse.  The  target  angular  position  is  assumed  to  be  that  of  the 
commanded  axis  of  the  receive  beam.  Many  real  or  proposed  radar  systems 
operate  without  monopulse  in  the  search  mode,  so  this  is  a realistic 
model. 


Bulk  filtering,  which  is  a part  of  the  search/verify  procedure  in 
many  real  radar  systems,  is  not  modeled  in  ROSCOE. 
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It  is  assumed  that  the  "real  target"  is  known  and  distinguished 
from  other  targets.  In  particular,  when  multipath  is  present,  only  one 
of  the  target  images  is  detected  by  this  procedure;  the  method  by  which 
this  target  is  to  be  selected  is  described  later  in  this  discussion.  When 
real  multiple  targets  are  present  in  the  beam(s) , the  association  and 
filtering  necessary  to  distinguish  the  correct  new  target  are  not  modeled, 
but  are  assumed  to  be  done  perfectly.  Extraneous  targets,  whether  real 
or  Images,  appear  principally  in  the  form  of  "noise"  affecting  the  return 
from  the  "real"  target. 

The  search  sector  is  a sector  (between  specified  azimuths)  of  a 
volume  of  revolution  whose  vertical  cross  section  has  the  shape  shown  in 
Fig.  3.1,  defined  by  maximum  and  minimum  ranges,  altitudes,  and  elevation 
angles.  Note  that  these  quantities  cannot  all  be  independently  specified; 
to  prevent  a gap  at  the  range/altitude  transition,  it  is  required  that 


Figure  3.1.  Vertical  Slice  Through  Search  Sector  with  Range-Height 
Transition 


28 


\ 

t 


t 


•.»W»  ^,,1,4. 


max,  alt.  m max,  range 
min.  alt.  min.  range 

The  range/altitude  transition  occurs  at  an  elevation  angle 

„ . -1  max.  alt. 

E„  - sin  

2 max.  range 

If  this  angle  exceeds  the  maximum  elevation,  or  if  maximum  altitude 
_>  maximum  range,  there  is  no  flat  portion  of  the  search  sector.  Such  a 
search  sector  might  be  used  by  an  area-defense  radar;  compare,  for  example, 
the  radars  in  the  two  scenarios  specified  in  the  appendixes  of  Vol.  1 
(Search  Mode  Parameters  Dataset).  An  airborne  or  orbiting  radar  might 
have  a negative  minimum  elevation  angle  for  search.  It  should  go  without 
saying  that  the  angular  extent  of  the  search  sector  for  any  radar  face 
should  lie  within  the  field-of-view  specified  for  that  face.  In  the 
program,  it  is  implicitly  assumed  that  this  is  the  case. 

If  the  target  appears  within  the  search  sectors  of  more  than  one 
radar  face,  detection  is  first  attempted  by  that  face  for  which  the  target 
has  the  smallest  of f-boresight  angle  at  the  time  it  crosses  the  outer 
boundary  of  the  search  sector. 

The  time  of  flight  of  radar  signals  is  not  modeled. 

3.1.2  Procedures 

The  user  may  elect  to  make  the  first  try  at  detection  as  soon  as 
the  target  crosses  the  outer  boundary  of  the  search  sector,  thus  ignoring 
the  random  orientation  of  the  search  frame  relative  to  the  target  position. 
Alternatively,  he  may  select  the  option  of  modeling  this  "frame  randomiza- 
tion" by  making  the  first  try  at  detection  at  some  randomly-selected 
instant  during  the  first  search  frame  period  after  the  target  crosses 
the  outer  boundary. 
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If  Che  first  try  at  detection  fails,  succeeding  tries  are  made  at 

intervals  of  the  search  frame  time,  T , after  the  first  try,  until 

8 

either  a successful  detection  is  made  or  the  target  crosses  the  inner 
boundary  of  the  search  sector.  In  the  latter  case,  detection  fails  to 
occur. 


The  search  frame  time,  Tg  , cannot  be  chosen  completely  at  random. 

It  must  be  short  enough,  compared  with  the  target  velocity  and  the  range 

window  (AR),  so  that  the  radar  gets  at  least  two  chances  at  any  target 

before  it  crosses  the  inner  boundary.  Since  the  first  chance  might 

(because  of  frame  randomization)  occur  as  much  as  T after  the  target 

8 

crosses  the  outer  boundary,  this  means  approximately  that 
2Tg  < AR/Vt 

where  is  the  target  velocity  at  the  time  it  crosses  the  outer  bound- 

ary. This  inequality  is  not  exact  and  is,  in  fact,  rather  conservative; 
in  many  cases,  considerably  larger  values  of  Tg  would  be  quite  satisfactory. 

When  successful  detection  is  accomplished,  a verification  pulse  is 

ordered  to  be  transmitted  at  time  T after  the  successful  detection 

v 

pulse.  The  verification  pulse  uses  the  same  beams  and  waveform  as  the 
search  pulse,  but  the  beam  axis  is  along  the  commanded  position  of  the 
beam  axis  for  the  successful  search  pulse.  If  the  verification  return  is 
above  the  search-return  threshold,  verification  is  assumed  to  be  accom- 
plished and  track  Initiation  is  begun.  If  the  verification  is  unsuccessful, 
search  continues  at  intervals  of  Tg  after  the  previous  "successful" 
detection. 

If  the  search  process  were  being  explicitly  modeled,  with  its  raster 
scan,  it  would  usually  be  the  case  that  the  target  appears  somewhere  off 
the  axes  of  the  transmit  and  receive  beams.  This  is  another  form  of  frame 
randomization,  and  it  is  modeled  in  ROSCOG  by  offsetting  the  axes  of  the 
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search-mode  transmit  and  receive  beams  with  respect  to  the  target  position. 
The  amount  and  direction  of  offset  are  randomly  generated  in  such  a way 
that  the  target  will  be  somewhere  within  the  3 dB  beamwidths  of  these 
beams,  allowing  for  refraction. 

Refraction.  The  term  "allowing  for  refraction"  is  interpreted  in 
the  following  ways  under  the  stated  conditions: 

1.  When  only  one  ray  path  exists.  The  ray  path  from  the  true 
target  position  to  the  radar  is  computed.  The  true  azimuth 
and  elevation  at  the  radar  are  used  for  the  apparent  values 
A'  and  E'  ; the  apparent  range,  R'  , is  computed  with 
allowance  for  both  the  added  length  of  the  curved  ray  path 
and  the  decreased  value  of  c In  plasma. 

2.  Multiple  paths  through  a single  fireball  or  stria.  When 
multiple  paths  are  present,  one  of  them  will  always  be  the 
analytic  continuation  of  the  single  ray-path  occurring  at 
large  distances  from  the  plasma.  This  ray-path  is  the 
one  to  use  for  both  outgoing  and  incoming  signals;  its 
azimuth  and  elevation  8t  the  radar  give  A'  and  E1  , and 
the  apparent  range  along  it  gives  R*  ; this  is  called  the 
"apparent  position  of  the  true  target." 

3.  Multiple  paths  through  a striated  region.  In  ROSCOE,  a 
"fuzz-ball"  approximation  is  used  (see  Sec.  4.3  for  further 
details) , with  R' , A' , E'  being  essentially  the  center  of 
the  apparent  target  image  distribution  which  is  the 
"fuzz-ball." 

3.1.3  Flow  Charts 

Definitions.  The  flow  charts  refer  to  several  sets  of  object  posi- 
tion coordinates.  In  general,  R,  A,  E refer  to  range,  azimuth,  and 
elevation  with  respect  to  the  radar,  and  R,  u,  v refer  to  "sine-space" 
coordinates,  where  u and  v are  the  sines  of  two  orthogonal  angles 
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measured  from  the  radar  boresight  (see  Appendix  A) . These  coordinates 
are  subscripted  or  primed  to  indicate  the  following: 


R ,u 

,v 

s s 

s 

R ,u 

,v 

V V 

V 

R'  ,A' ,E'  or  R' ,u* ,v*  Apparent  target  coordinates  allowing  for 

refraction  but  not  including  radar  measure- 
ment errors 

Measured  coordinates  from  search  return 

Measured  coordinates  from  verification  return 

Sine-space  coordinates  of  transmit  beam  axis 

Sine-space  coordinates  of  receive  beam  axis 

Other  symbols  used  in  the  flow  charts  are: 

6 Angle  between  line-of-sight  to  target  and  boresight 

acU(x,y)  Random  number  chosen  from  uniform  distribution  between 
x and  y 


(VVo>T 

'vVr 


Ii  ltlallzation  (Fig.  3.2).  The  initialization  procedure  is  not  an 
event,  but  takes  place  before  any  actual  events  have  occurred.  Its 
purpose  is  to  set  up  the  first  search  event  for  each  target.  An  implicit 
assumption  is  that  search  and  verification  are  being  performed  by  the 
radar,  rather  than  by  some  other  kind  of  sensor. 


In  the  initialization  procedure,  the  time  for  the  first  attempt  at 
detection,  tQ  , is  computed  in  vacuum,  without  reference  to  refractive 
effects.  This  of  course  will  not  give  the  strictly  "correct"  time,  but 
it  permits  the  search  events  to  be  time-ordered  in  advance  without  calling 
upon  nuclear-effects  and  propagation  subroutines.  The  time  error  involved 
is  not  likely  ever  to  be  significant  in  the  context  of  program  applica- 
tions, and  the  approximation  makes  it  possible  to  maintain  a completely 
event -based  structure. 


Search  (Fig.  3.3).  The  search  event  is  first  set  up  by  the  ini- 
tialization procedure,  and  itself  sets  up  the  appropriate  succeeding  event. 
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Figure  3.2.  Search  Initialization  Flow  Chart 
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Figure  3.3.  Search  Flow  Chart 
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Verification  (Fig.  3.4).  The  verification  event  is  set  up  by  a 
successful  search  event. 

3.2  TRACK  INITIATION  EVENT 

3.2.1  Assumptions 

Only  one  radar — the  same  one  that  was  used  for  search — is  used  for 
track  initiation. 

Because  of  the  event  structure,  any  number  of  targets  may  pass  (one 
at  a time)  through  the  track  initiation  procedure. 

The  track  initiation  procedure  uses  its  own  pre-designated  waveform 
and  beamshape,  which  may  (but  need  not)  be  the  same  as  those  for  the 
search  event.  For  track  initiation,  the  beam  size,  shape,  and  axis  are 
the  same  for  both  transmit  and  receive  beams.  The  waveform  is  usually 
chirped. 

Monopulse  i£  used  in  the  track  initiation  procedure. 

A general  rule,  applicable  to  this  as  well  as  all  other  radar 
events,  is  that  each  separate  transmitted  radar  signal  is  a separate  event. 
(A  "signal"  may,  of  course,  consist  of  more  than  one  pulse,  provided  that 
they  are  sent  as  a unit  under  a single  command — e.g.,  pulse  pairs  or 
pulse  bursts.)  This  event  structure  may  appear  to  be  less  efficient  than 
lumping  together  several  transmissions  which  are  a logical  unit  (such  as 
the  search  and  succeeding  verification  pulse,  or  the  various  track  ini- 
tiation pulses  to  be  described  below) , but  it  avoids  a possible  problem 
arising  from  such  lumped  events.  This  problem  would  occur  when  an  event 
related  to  a different  target  or  radar  should  be  scheduled  for  a time 
in  the  midst  of  the  lumped  period  occupied  by,  for  example,  the  five 
possible  track  initiation  pulses.  This  would  distort  the  proper  time- 
ordering of  events  so  necessary  to  the  smooth  operation  of  ROSCOE. 
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Another  general  rule  is  that  each  event  is  coincident  in  time  with 
the  instant  of  transmission  of  the  relevant  radar  signal.  The  tirae-of- 
flight  of  the  signal  between  the  radar  and  the  target  is  ignorea. 

3.2.2  Procedure 

The  first  track  initiation  pulse  is  transmitted  T seconds  after 
the  successful  verification  pulse,  where  T is  the  period  between  track 
pulses  (an  input  quantity).  Succeeding  track  initiation  pulses  will  be 
transmitted  at  intervals  of  T , whether  the  first  one  was  successful 
or  net. 

Two  good  track  initiation  returns  are  needed  to  start  track.  If 
the  first  track  initiation  pulse  is  not  successful,  a second  attempt  is 
made  Tc  seconds  later.  These  two  pulses  are  referred  to  as  Til  pulses. 

If  neither  of  them  is  successful,  the  program  reverts  to  the  search  process 
in  the  regular  schedule. 

If  there  is  a successful  Til  return,  another  track  initiation  pulse 
is  sent  to  the  appropriate  place  (see  below)  ; this  pulse  is  referred  to 
as  a TI2  pulse,  although  it  is  exactly  the  same  in  waveform  and  beamshape 
as  the  Til  pulses.  If  a successful  return  is  not  obtained  from  the  first 
TI2  pulse,  a second  is  sent  after  T ; if  that  is  not  successful,  a third 
attempt  is  made  after  the  same  interval.  At  this  point,  failure  to  obtain 
a successful  return  causes  the  program  to  revert  to  the  search  procedure 
in  the  regular  schedule.  If  there  is  a successful  TI2  return,  a track 
event  is  scheduled  T^.  seconds  later. 

Thus,  under  the  best  circumstances  one  Til  pulse  and  one  TI2  pulse 
are  transmitted,  whereupon  the  track  process  can  begin.  Under  the  worst 
circumstances,  two  Til  pulses  and  three  TI2  pulses  could  be  transmitted 
before  admitting  defeat  and  reverting  to  the  search  procedure. 
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In  a real  system,  the  radar  would  continue  to  detect  the  object  with 
search  pulses  as  long  as  It  remained  within  the  search  sector,  even  after 
track  had  commenced.  Such  detections  would  be  compared  with  existing 
track  files,  and  ignored  when  an  association  had  been  established.  In 
ROSCOG,  it  is  assumed  that  such  an  association  process  is  always  perfectly 
rerformed,  and  no  additional  search  events  are  scheduled  for  an  object 
once  track  initiation  has  been  completed.  In  other  words,  possible  errors 
in  association  leading  to  multiple  track  files  on  a single  object  are  not 
modeled. 

On  the  other  hand,  if  two  track  initiation  returns  are  not  success- 
fully completed,  ROSCOG  reverts  to  the  search  process. 

The  first  Til  pulse  is  transmitted  in  the  direction  of  the  com- 
manded axis  of  the  verification  pulse,  since  (in  the  absence  of  monopulse 
for  the  verification  procedure)  that  is  the  best  estimate  of  the  target's 
angular  coordinates. 

The  receiver  range  gate  for  this  pulse  has  width  (an  input 

quantity)  and  is  centered  on  the  estimated  range 


whe'  -»  R * measured  range  from  the  verification  return 

v 5 

R - measured  range  from  the  preceding  successful 
8 

search  return 

T^  • time  interval  between  the  search  and  verify  pulses 

T ■ time  interval  between  the  verify  pulse  and  the 
first  Til  pulse 
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The  second  Til  pulse.  If  needed,  has  Its  range  gate  centered  at  the 
estimated  range 


The  gate  width  remains  the  same. 

The  first  TI2  pulse  Is  transmitted  in  the  direction  Indicated  by 
the  monopulse  return  of  the  successful  Til  pulse.  The  receive  gate  again 
has  width  , and  is  centered  at  the  estimated  range 


where  - measured  range  from  the  successful  Til  pulse 

Rg  - measured  range  from  the  last  successful  search  pulse 

t^  - transmission  time  of  the  successful  Til  pulse 

t ■ transmission  time  of  the  last  successful  search  pulse 

If  needed,  the  second  and  third  TI2  pulses  have  range  gates  centered 
at 


respectively,  where  n ■ 2 or  3 as  required. 

'n  computing  the  return  from  trick  initiation  and  track  pulses,  all 
applic  Me  object,  radar,  environmental,  and  propagation  effects  are  taken 
into  account.  These  include  refraction,  multipath,  focusing  or  defocusing, 
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absorption,  dispersion,  noise,  and  clutter.  The  "signal-to-noise"  ratio, 
S/N,  referred  to  in  several  places  in  the  flow  charts,  should  more  properly 
be  called  the  "signal-to-noise-plus-clutter"  ratio,  S/(N  + C)  „ since 
clutter  is  included;  see  Sec.  4.5. 

Initializing  the  Tracking  Filter.  The  tracking  filter  is  initialized 
after  a successful  TI2  pulse.  The  state  of  the  target  as  handed  over  to 
the  filter  is  in  the  Geocentric  Rectangular  Cartesian  (G.R.C.)  coordinate 
system  (see  Appendix  A).  It  consists  of  the  following  elements: 

t j * time  for  which  the  target  state  was  determined 
(transmission  time  of  the  successful  TI2  pulse) 

-*■ 

T ■ measured  position  vector  of  the  target  at  time  t2 , 
with  components  (T^.T^.T^) 

i 

T - calculated  rate-of-change  of  target  position  at 
time  t2>  with  components  (T^.^.T^) 

aT’  °x  “ estimates  of  the  measurement  errors  in  position  and 
velocity  of  the  target  at  time 

It  must  be  emphasized  here  that  the  values  of  o are  decidedly  not 
supposed  to  be  the  actual  radar  measurement  errors  computed  by  ROSCOE 
for  the  TI2  pulse,  but  rather  are  estimates  which  a system  would  make 
on  the  basis  of  available  Information.  The  computation  of  these  esti- 
mates is  discussed  in  Sec.  3.2.3. 

Setting  Up  Discrimination  Events.  Once  the  tracking  filter  has 
been  initialized  and  a track  file  has  been  opened  on  the  target , it  is 
possible  to  set  up  future  discrimination  events  (see  Sec.  3.4). 

3.2.3  Details 


3.2. 3.1  Interactions  with  the  Tracking  Filter  Subroutine 

The  tracking  filter  subroutine,  FILTER,  has  a number  of  variant 
calling  sequences  depending  on  its  application  at  the  calling  time.  The 
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subroutine  is  assumed  to  know  what  the  value  of  the  calling  time  is;  in 
the  present  case,  it  is  t^  . The  calling  sequence  for  the  initialization 
call  is 

(PV,  SIGMAT,  DUMMY,  DUMMY,  DSPWRD,  DUMMY) 

The  dataset  pointer  DSPWRD  is  set  to  zero,  indicating  that  this  is 
the  initial  call  to  the  tracking  filter  for  this  target. 

The  Cartesian  vector  PV  is  a six-dimensional  vector  containing  the 
position  and  velocity  of  the  target  at  time  t2  » as  measured  and  cal- 
culated from  the  two  TI  pulses  and  transformed  into  G.R.C.  coordinates; 
i.e.,  the  vector  ,T2  .T^T^T^  ,T3)  . 

The  covariance  matrix  SIGMAT  is  a 6x6  matrix  which  contains  the 
initial  estimate,  in  G.R.C.  coordinates,  of  the  errors  in  PV.  For  reasons 
that  will  be  discussed  in  Sec.  3. 2. 3. 2,  only  a rough  estimate  of  the 
diagonal  elements  of  this  matrix  can  be  made;  the  off-diagonal  elements 
are  not  known  at  all  and  are  set  equal  to  zero.  The  components  of  SIGMAT 
are: 

SIGMAT(I , J)  - 0 for  I l1  J 

SIGMAT(l.l)  - SIGMAT (2  ,2)  - SIGMAT  (3,3)  - o* 

SIGMAT(4 ,4)  - SIGMAT(5 ,5)  - SIGMAT(6,6)  - oj 

The  track  initiation  procedure  also  initializes  two  parameters  that 
are  maintained  in  the  track  event:  TLAST,  which  indicates  the  time  of 

the  last  successful  track  return  for  this  radar/target  combination, 
and  ILT,  which  indicates  the  number  of  successive  unsuccessful  track 
pulses  for  this  radar/target  combination.  This  information  is  needed 
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in  case  of  lost  track  on  the  target,  when  It  is  desired  to  keep  a record 
of  how  long  it  has  been  since  the  target  was  last  seen.  Both  parameters 
are  initialized  at  zero  by  the  track  initiation  procedure. 

After  the  initial  call  to  subroutine  FILTER,  the  data  needed  for 
the  track  file  to  be  established  for  this  target  has  been  prepat  *?d,  and 
their  location  is  indicated  by  the  returned  value  of  the  DSP  word  DSPWRD. 

Finally,  a first  track  event  is  set  up  for  this  radar /target  com- 
bination at  time 


The  symbol  "t^"  will  be  used  consistently  to  represent  the  time  of  a 
track  event. 

3.2. 3.2  System  Estimates  of  Measurement  Errors 

The  radar  modeling  in  ROSCOE  of  course  Introduces  measurement 
errors,  both  bias  and  random,  in  each  of  the  radar  coordinates.  In  cal- 
culating the  values  of  these  errors  and  applying  them  to  the  "apparent" 
coordinates  to  yield  the  "measured"  coordinates,  the  ROSCOE  model  is  playing 
the  part  of  Nature,  whose  errors  are  completely  known  to  herself  but  can 
only  be  imperfectly  estimated  by  the  radar  system.  This  modeling  of 
Nature's  errors  is  a part  of  the  radar  measurements  simulation,  and  will 
be  fully  discussed  in  Sec.  4.7. 

But  ROSCOE  needs  another  kind  of  measurement  error  model  also:  this 

is  the  model  which  the  radar  system  would  employ  in  estimating,  given 
imperfect  Information,  what  the  measurement  errors  might  be  for  any 
particular  return.  Here  ROSCOE  plays  the  part  not  of  Nature  but  of  the 
radar  engineer.  He  has  certain  formulas,  developed  according  to  various 
theoretical  models,  for  estimating  the  probable  error  inherent  in  a given 
measurement.  Such  estimated  errors  are  used  by  the  track  filter  to  decide 
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how  Co  weight  each  new  measurement  in  comparison  with  the  old  ones.  They 
are  also  needed  for  properly  Initializing  the  track  filter.  This  section 
discusses  a simple  model  which  is  used  in  ROSCOE  for  making  these  error 
estimates. 

The  track  filter  requires  estimates  of  the  measurement  errors  in 
each  of  the  primary  coordinates:  R,  u,  v . Range  or  angular  rates  are 

not  directly  used  in  this  filter,  so  their  errors  need  not  be  estimated. 
Furthermore,  since  the  filter  uses  them  In  the  same  coordinate  system 
in  which  they  are  measured,  there  is  no  need  to  transform  the  error 
estimates  into  another  coordinate  system. 

For  track  initiation,  however,  things  are  not  so  simple.  Not  only 
the  position  but  also  the  velocity  of  the  target  is  needed  to  start  off 
the  tracking  filter,  and  these  must  be  expressed  in  the  G.R.C.  coordinate 
system  (see  Appendix  A).  The  corresponding  error  estimates  must  therefore 
also  be  calculated  and  transformed  into  that  frame.  Procedures  for  doing 
this  are  described  below. 

Position  Error  Estimates.  Three  types  of  errors  are  usually  dis- 
tinguished In  each  measured  coordinate:  a bias  error,  a jitter  error 

component  which  is  independent  of  signal-to-nolse  ratio,  and  a jitter 
error  component  which  is  slgnal-to-nolse  dependent.  The  bias  error  is 
represented  in  the  "Nature's  model"  part  of  ROSCOE,  but  cannot  be  esti- 
mated by  the  system  (if  it  could,  being  constant,  it  would  be  removed). 

The  jitter  error  in  each  radar  coordinate  is  modeled  as  being  Gaussian 
distributed  with  zero  mean  and  a standard  deviation,  o , of  the  general 
form 


~2 

a 


-2  . -2 
°l  + 02 


where 


is  a constant  and  o^ 


K_ 

/s7n 
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Here  K is  a constant,  dependent  on  the  beamshape  or  waveform  parameters, 
which  is  given  in  all  standard  texts  on  radar  theory. 

Although  there  are  a number  of  good  theoretical  formulas  for  the 
constant  K , there  are  no  general  theoretical  models  for  the  constant 

o.  , because  this  is  supposed  to  be  dependent  on  specific  mechanical  and 

A 4 

electronic  limitations  in  the  radar  itself.  However,  it  has  been  observed 

that  in  many  existing  or  planned  radars  the  ratio  K/o^  falls  within 

relatively  narrow  limits  (possibly  because  the  designers  all  use  the  best 

current  state-of-the-art).  By  using  a mean  value  for  this  ratio,  y 

(which  is  different  for  linear  and  angular  coordinates),  can  be 

estimated  given  K . The  resulting  formula  for  a accordingly  takes 

the  form 


The  constant  K is  usually  quoted  in  terms  of  the  altazimuth  coordinate 
system,  i.e.,  in  range  and  angle.  The  angular  error,  however,  then 
depends  on  the  off-boreslght  angle  and  the  consequent  effective  beam 
broadening  as  well  as  beamshape  distortion.  These  effects  do  not  arise 
in  the  sine-space  (R,  u,  v)  coordinate  system,  which  is  one  of  the  main 
reasons  why  ROSCOE  uses  it.  In  R,  u,  v coordinates,  the  following  values 
of  K are  used: 


*R 


K 


u 

K 


v 


c 


B 


speed  of  light  in  vacuum 
effective  signal  bandwidth 


where 


6w^  ■ sin(68^/2) 

66^  ■ one-way  3-dB  beamwldths  of  receive  bean  in  the 
u and  v directions 


Velocity  Error  Estimates.  The  velocities  R,  u,  v are  determined 
for  track  initiation  purposes  by  taking  differences  of  the  measurements 
made  with  the  Til  and  TI2  pulses,  as  previously  discussed. 


The  error  in  computing  the  rate-of-change  of  a measured  coordinate 

from  a series  of  n independent  measurements  of  that  coordinate  has  been 
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estimated  by  Manasse.  Although  his  derivation  is  not  strictly  applicable 

to  the  present  case  where  n ■ 2 , ROSCOE  uses  a modified  form  of  his 
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result  and  allows  for  the  smallness  of  n by  dropping  the  factor  n 
from  his  equation.  The  result  is 


where  p is  the  coordinate  Involved  (R,  u,  or  v as  the  case  may  be), 

T is  the  time  interval  between  the  two  measuring  pulses,  and  5 is  the 

P 

estimated  error  in  measurement  of  p . Accordingly,  ROSCOE  uses 


- _ 2/3  - 

°ft  Tt  °R 


where  Tfc  , 


as  before,  is  the  interval  between  track  initiation  pulses. 


Transformation  to  G.R.C,  Coordinates.  It  is  not  in  general  possible 
to  make  a meaningful  transformation  of  error  sigmas  from  one  coordinate 
system  to  another  when  the  coordinate  transformation  is  nonlinear.  The 
actual  error  ellipsoid  remains  the  same,  of  course,  so  in  effect  the 
error  sigmas  in  the  transformed  coordinates  can  have  any  values  between 
the  smallest  and  the  largest  linear  radar  measurement  errors.  The  obvious 
simple  and  safe  course  is  to  make  the  most  conservative  estimate,  which 
is  to  assume  that  the  error  sigmas  in  G.R.C.  coordinates  (W^,  W0 , U\) 
are  all  the  same  and  equal  to  the  largest  linear  measurement  sigma  (and 
similarly  for  the  velocities).  This  leads  to  the  following  expressions: 

\ ' 5W2  ■ \ * ’“‘x[SR’R8u’RSv] 

\ ■ \ ' a*3  2 •“[5ft’RVR5o] 

Thus  only  the  diagonal  elements  of  the  6x6  error  covariance  matrix  can 
be  estimated,  and  those  only  roughly.  The  off-diagonal  elements  are  not 
even  guessed  at,  but  are  set  equal  to  zero. 

3.2.4  Flow  Charts 

Flow  charts  for  the  Til  and  TI2  events  are  given  in  Figs.  3.5  and 

3.6. 


3.3  TRACK  EVENT 

3.3.1  Assumptions 

3. 3. 1.1  Logical  Flow 

ROSCOE  generally  implements  a simplified,  straight-line  logical 
flow  from  one  procedure  to  the  next.  Only  one  radar  face  does  Search/ 
Detection/Verification;  if  the  radar  has  more  than  one  face,  the  S/D/V 
initialization  routine  selects  for  each  target  the  radar  face  which 
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Figure  3.6.  Track  Initiation  (TI2  Pulses)  Flow  Chart 
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would  detect  it  first.  The  user  can  also  control  which  radar  face  does 
the  Initial  detection  by  specifying  that  only  one  of  them  can  search. 

The  logic  next  assumes  that  the  same  radar  that  did  the  initial 
detection  also  performs  the  track  initiation  on  the  detected  target. 

This  is  a realistic  assumption,  provided  that  one  is  only  concerned  with 
targets  for  which  no  track  file  has  been  established  until  detected  by 
a radar.  Ultimately,  however,  ROSCOE  will  have  to  consider  situations  in 
which  initial  detection  and  track  initiation  were  performed  by  optical 
sensors  before  any  radar  detection  could  take  place.  The  radar  logic 
would  then  have  to  bypass  the  radar  search,  verification,  and  track  ini- 
tiation routines  and  start  right  in  with  tracking  on  this  target  whose 
track  file  has  already  been  established.  This  would  constitute  a "multiple- 
origin,"  or  non-straight-line  logical  flow.  The  system  here  described 
does  not  yet  handle  such  a situation. 

3. 3. 1.2  The  Tracking  Filter 

Tracking  filters,  in  general,  can  suffer  from  stability  problems; 
filters  which  work  well  in  a normal  environment  can  become  unstable  in  a 
nuclear-disturbed  environment.  In  such  situations,  certain  filter  param- 
eters may  have  to  be  "tuned"  in  order  to  get  a filter  to  perform  reason- 
ably well.  Or  one  can  turn  off  the  filter  (l.e. , refuse  to  accept  further 
data)  for  a certain  period  after  a nuclear  burst  in  the  vicinity  of  the 
line-of-slght.  Or  one  can  allow  the  filter  time  constant  and/or  data 
weighting  factors  to  be  continuously  adjustable  parameters  in  response 
to  some  observable  which  (it  is  hoped)  is  correlated  with  the  degree  of 
environmental  degradation  being  encountered.  Ultimately,  one  may  attempt 
to  design  special  tracking  filters  for  use  in  refractive  environments — 
for  example,  filters  which  Include  refraction  error  parameters  among  their 
fitted  variables. 
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All  of  these  schemes,  and  others,  have  been  under  active  considera- 
tion in  the  community.  Some  users  may  wish  to  use  ROSCOE  as  a vehicle  to 
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test  out  the  effectiveness  of  their  own  filters.  If  so,  the  modular 
arrangement  of  ROSCOE  will  permit  them  simply  to  substitute  their  filter 
for  the  standard  ROSCOE  filter  with  a minimum  of  difficulty.  But  it  must 
be  clearly  understood  that  the  standard  ROSCOE  filter  is  not  intended  to 
be  general  enough  to  permit  all  these  contemplated  adjustments.  The 
ROSCOE  filter  is  a more-or-less  standard  Kalman  filter,  with  only  such 
adjustable  parameters  as  are  described  below. 

Only  the  requirements  on,  and  the  external  features  of,  the  track 
filter  are  given  in  this  subsection;  the  internal  operational  details  of 
the  filter  were  discussed  in  Sec.  2.4. 

Filter  Inputs.  ROSCOE' s filter  is  a Kalman  filter  operating  with 
seven  target  state  variables:  position,  velocity,  and  ballistic  coeffi- 

cient. The  input  data  may  be  two-dimensional  measurements  (e.g.,  angles- 
o.ily  from  an  optical  sensor),  three-dimensional  measurements  (e.g.,  R,  u, 
v from  a radar),  four-dimensional  measurements  (e.g.,  R,  u,  v,  and  R)  , or 
whatever  is  available.  Any  convenient  measurement  coordinate  system  can 
be  used;  its  nature  must  be  known  to  the  sensor  but  need  not  be  known  to 
the  filter.  It  is  necessary,  however,  that  the  sensor  module  know  how  to 
transform  a predicted  target  position  given  in  G.R.C.  coordinates  into 
the  measurement  system.  The  present  discussion  assumes  that  all  radar 
measurements  are  made  in  radar  face  coordinates  (R,  u,  v) . 

In  the  filter,  the  input  measurer-  nts  are  weighted  in  accordance 
with  their  estimated  error  sigmas,  which  must  also  be  input  to  the  filter. 
Section  3. 2. 3. 2 explains  how  to  calculate  5^,  5u,  and  oy. 

The  time  at  which  the  measurement  was  made  must  also  be  an  input. 

Filter  Outputs.  The  filter  outputs  needed  by  the  track  event 
routines  are: 

• The  estimated  target  state  vector  at  the  time  of  the 
measurement,  in  G.R.C.  coordinates. 
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The  filter's  current  estimates  of  the  errors  in  these 
coordinates,  derived  from  the  filter's  internal  covariance 
matrix. 


Filter  Memory  Decay.  It  is  convenient,  though  not  absolutely 
necessary,  that  the  filter  has  an  exponential  memory  decay.  Without 
such  iecay,  older  data  may  be  given  too  much  weight  in  comparison  with 
new  data,  leading  to  unnecessarily  prolonged  propagation  of  bias  or 
start-up  errors,  errors  resulting  from  temporary  refractive  effects,  or 
errors  due  to  imperfectly  known  atmospheric  models  and  ballistic  coeffi- 
cients. On  the  other  hand,  if  the  decay  time  is  too  short  the  target 
may  be  lost  during  prolonged  periods  of  missing  returns  (such  as  those 
due  to  blackout)  or  refractive  errors.  Thus  the  filter  decay  constant 
t (if  present)  is  one  of  the  more  sensitive  parameters  in  terms  of  which 
the  filter  can  be  "tuned"  to  function  best  under  given  conditions. 


The  following  two-stage  scheme  is  used  for  defining  the  filter 
decay  time  constant:  the  user  inputs  an  altitude,  , and  two  values 

of  the  decay  constant,  x^  and  x^  - The  filter  uses 


T ■ 


when  H > H 

p - T 


when  H < H 
P 


where  H is  the  target  altitude  as  determined  from  the  filter  output 
(not  from  raw  measured  data) . 


Note  that  the  user  can  set  ■ 0 and  * 00  (i.e.,  some  very 

large  number)  if  he  wishes  to  use  no  filter  decay.  Similarly,  setting 
■ 0 permits  the  user  to  specify  only  a single  decay  constant  for  the 
entire  track.  It  might  also  be  pointed  out  that  the  above  decay  constant 
definition  scheme  is  for  the  filter  itself  and  hence  applies  uniformly 
to  all  radars  in  the  system. 
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Other  Adjustments.  The  standard  ROSCOE  filter  weights  input  data 
in  accordance  with  the  estimated  measurement  error  sigmas,  but  does  not 
apply  any  sort  of  special  weighting  based  on  measured  or  assumed  environ- 
mental conditions. 

Rv'SCOE  does  not  have  any  provision  for  "turning  off"  the  filter  in 
response  to  nuclear  bursts:  neither  an  "exclusion  time"  nor  an  "exclusion 

volume"  nor  any  other  such  process  is  employed. 

3.3.2  Procedure 


3. 3. 2.1  Multiple  Radars 

In  a real  system,  several  radars  may  detect  the  same  target  during 
search.  Each  one  will  start  a separate  track  file  on  it,  but  presumably 
the  system  logic  will  at  some  point  make  the  appropriate  associations 
and  coalesce  them  all  into  a single  track  file.  Thereafter  all  radar 
d-.ta  on  this  target  would  go  into  the  same  tracking  filter  file,  regardless 
of  which  radar  took  the  measurements. 

Environmental  effects  could  certainly  have  an  adverse  effect  on 
this  association  process,  probably  resulting  in  the  prolonged  maintenance 
of  multiple  track  files  on  a single  target.  But  although  this  is  Indeed 
an  interesting  problem  and  one  related  to  nuclear  effects,  its  solution 
lies  In  the  realm  of  system  logic  and  is  outside  the  purview  of  ROSCOE. 
Accordingly,  ROSCOE  ignores  the  track  association  problem  and  maintains 
a separate  track  file  for  each  radar/object  pair. 

Other  radars  (if  present)  may  join  in  the  track  process,  all  main- 
taining separate  track  files.  The  track  event  logic  is  structured  in  such 
a way  that  any  radar  may  be  making  the  present  measurement,  according  to 
a previously-set-up  track  event  which  told  it  where  and  when  to  transmit 
its  signal.  The  logic  also  schedules  the  next  track  event  for  this  radar/ 
target  combination. 
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3. 3. 2. 2 Track  Cut-Off  Range 

Some  radars  have  a minimum  range  inside  of  which  they  cannot  track 
a target.  This  option  is  included  in  the  model:  the  user  specifies 
Rtmin  f°r  each  radar.  Then,  as  soon  as  the  range  of  a target  (as 
obtained  from  the  tracking  filter  projection)  is  less  than  Rttn^n  » no 
further  track  events  are  scheduled  for  that  radar  on  that  target. 

3. 3.2. 3 Tracking  Range  Gates 

A real  radar  system  listens  for  a return  only  during  some  specified 
time  interval  spanning  the  estimated  time  for  the  track  pulse  return. 

This  Interval  is  called  the  range  gate,  and  is  usually  expressed  in  range 
units  rather  than  time  units.  It  needs  to  have  sufficient  width  to  achieve 
a high  probability  of  receiving  the  target  return,  yet  be  small  enough  to 
eliminate  many  extraneous  echoes  and  minimize  the  amount  of  time  the  radar 
has  to  spend  in  listening.  To  facilitate  this,  the  width  of  the  range 
gate  is  usually  given  as  a function  of  the  estimated  range  variance  as 
obtained  from  the  tracking  filter. 

Before  the  filter  has  been  started  (i.e.,  during  track  initiation), 
the  range  gate  width  is  a specified  input  value,  . Thereafter,  the 
range  gate  is  defined  as  follows.  The  user  inputs  parameters  k^  and 
for  each  radar.  On  the  basis  of  all  previous  track  pulses,  an  esti- 
mated (Cartesian)  state  vector  for  the  target  has  been  computed  by  the 
tracking  filter,  along  with  an  estimate  of  the  covariance  matrix  of  the 
state  parameters.  These  are  projected  ahead  to  the  time  tQ  for  the 

current  track  event,  and  the  predicted  range  R and  range  variance 
o P 


Oa  are  calculated  for  t 
P c 

the  range  window 


Then  the  track  return  is  listened  for  in 


where 


R < R < R + ~ 


H Vr 
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If  the  measured  range,  after  allowing  for  refraction  and  noise  errors,  does 
not  fall  within  this  window,  the  return  is  recorded  as  "not  received." 

3. 3. 2. A Lost  Targets 

In  a real  radar  system,  it  may  happen  that  a number  of  consecutive 
track  pulses  fail  to  produce  a return — a "lost  track."  If  all  radars  that 
were  tracking  the  target  lost  track  on  it,  then  the  system  is  said  to  have 
"lost  target."  In  this  event,  a real  radar  system  will  initiate  some 
special  course  of  action  (such  as  a limited  search  scan  about  the  predicted 
position)  in  an  attempt  to  recover  the  target.  Such  system  logic  is  not  a 
proper  concern  for  ROSCOE,  and  such  specific  responses  are  not  modeled. 

In  ROSCOE,  the  following  procedure  is  used.  If  no  successful  track 
returns  have  been  received  on  a given  target  for  a time  period  AT^T  (an 
input  parameter) , the  target  is  assumed  to  be  "lost"  and  a new  search  pulse 
is  scheduled.  These  will  continue  until  either  the  target  is  reacquired 
or  it  ceases  to  exist  (impacts  or  detonates). 

The  method  described  above  gives  the  user  a step-by-step  picture  of 
what  the  nuclear  environment  is  doing  to  his  target  tracking  attempts. 

3.3.3  Definitions 

(R,u,v)  represent  true  vacuum  coordinates  of  the  target  relative  to 
the  radar  face,  with  no  propagation  effects  or  radar  measurement  errors 
Included.  They  are  computed  by  transforming  the  true  target  state  vector 
from  G.R.C.  to  radar  face  coordinates. 

(R'.u'.v')  represent  the  apparent  target  coordinates  relative  to  the 
radar  face,  allowing  for  all  propagation  effects  but  not  including  radar 
measurement  errors.  They  are  computed  by  applying  appropriate  propagation 
interactions  to  (R,u,v)  and  then  calculating  what  the  radar  "thinks"  it 
sees,  taking  into  account  the  beamshape  and  the  commanded  position  of  the 
beam  axis. 
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(R  ,u  ,v  ) represent  the  measured  target  coordinates  relative  to 
m m m ________ 

the  radar  face,  including  all  radar  measurement  errors.  They  are  computed 
from  (R'.u’.v1)  by  taking  into  account  the  current  value  of  the  signal- 
to-noise-plus-clutter  ratio. 

(Rp,Up,Vp)  represent  the  predicted  target  coordinates  relative  to 
the  radar  face,  calculated  from  the  estimated  target  state  vector  from 
the  track  file,  updated  to  time  t^  . The  calculation  includes  no  pro- 
pagation effects  or  radar  measurement  errors;  it  is  made  by  transforming 
the  updated  estimated  state  vector  from  G.R.C.  to  radar  face  coordinates. 

(R,u,v)  represent  the  "expected11  target  coordinates  relative  to 

the  radar  face,  for  use  by  the  track  filter.  After  the  filter  has  absorbed 

the  current  measurements  (R  ,u  ,v  ) it  returns  a Cartesian  state  vector 

m m m 

CARTEZ  giving  its  estimate  of  where  the  target  should  now  be;  (R,u,v)  are 

calculated  directly  from  this  vector  with  no  propagation  effects  or  radar 

measurement  errors  included,  as  above.  This  follows  the  procedure  laid 

down  for  using  the  track  filter;  the  coordinates  (R,u,v)  are  fed  back 

into  the  filter  to  be  compared  with  (R  ,u  ,v  ) in  order  to  compute  the 

m m m 

next  iteration  of  the  filter  variables.  This  may  have  to  be  done  several 
times  before  the  filter  is  satisfied  with  the  result.  For  further  details, 
see  below,  and  also  Sec.  2.4. 

3.3.4  Inputs 

The  following  data  are  required,  either  as  inputs  or  in  datasets, 
for  the  track  event  procedure: 

• The  radar  identification  (IR)  and  location 

• The  radar  face  number  (IFACE) 

• The  target  identification  (IT) 

• The  transmission  time  of  this  track  pulse  (t^) 

• Dataset  pointers  to 

• The  true  trajectory,  ballistic  coefficient,  and 
radar  cross  section  parameters  of  the  target 
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• The  track  filter  estimates  of  the  target  trajectory 
and  ballistic  coefficient  parameters 

• The  current  track  filter  data 

• The  radar  waveform  and  beamshape  parameters  for 
this  radar 

• The  boresight  vectors  for  all  faces  of  this  radar, 
and  their  fields-of-view 

• The  track  detection  threshold  (K^,) 

• The  minimum  tracking  range  for  this  radar  (R  . ) 

tmin 

• Parameters  for  computing  the  range  gate  width  (k^^) 

• The  current  values  of  the  parameters  TLAST(IR,IT)  and 
ILT(IR,IT)  for  this  radar/target  combination 

• Decay-time  parameters  (t^,t2,Ht>  for  the  track  filter 

• The  time  period  between  successive  track  pulses  on  a given 
target  (Tt> 

• The  lost-target  parameter  for  this  system  (ATLT) 

• Flags  to  indicate 

• That  this  target  has  been  "killed" 

• That  this  radar  cannot  see  this  target  with  any 
of  it 8 faces 

3.3.5  Outputs 

Possible  outputs  of  this  track  event  procedure  are: 

• An  updated  track  file  on  this  target  and  the  scheduling  of 
the  next  track  event  for  this  radar/target  combination 

• The  scheduling  of  the  next  track  event  for  this  radar/ 
target  combination  without  updating  the  track  file  on 
this  target 

• Indication  that  this  radar  has  lost  track  on  this  target 
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3.3.6  Logical  Outline  of  the  Flow  Charts 

3. 3. 6.1  Basic  Procedure  (Fig.  3.7) 

The  basic  flow  through  this  event,  assuming  a successful  track 
return,  is  found  in  Boxes  1 through  24. 

If  the  target  has  been  "killed"  between  the  time  this  event  was 
scheduled  and  the  time  the  track  pulse  was  to  be  transmitted,  the  proce- 
dure is  aborted,  no  pulse  is  sent,  and  no  new  track  event  for  this  target 
and  radar  is  set  up  (Box  1) . 

If  the  predicted  target  cannot  be  seen  within  the  fields-of-vlew 
of  any  of  the  radar's  faces,  no  pulse  is  sent  at  this  time  (Boxes  25 
through  30) , and  the  procedure  goes  into  the  "no  return"  mode  (Boxes  32 
through  38;  see  below). 

If  the  predicted  target  range  is  less  than  the  minimum  tracking 
range  (if  any)  for  this  radar,  no  pulse  is  sent  and  no  future  pulses  are 
scheduled;  a "dropped  track"  flag  is  set  (Box  31). 

If  a pulse  is  transmitted  but  no  successful  return  is  received, 
either  because  the  return  did  not  exceed  the  track  threshold  (Box  11)  or 
because  the  measured  return  (including  refraction  and  errors)  fell  outside 
the  fleld-of-view  (Box  13)  or  because  the  measured  range  fell  outside  the 
receive  range  gate  (Box  14) , the  procedure  goes  into  the  "no  return"  mode 
(Boxes  32  through  38;  see  below). 

The  interaction  between  this  event  and  the  track  filter  subroutine 
is  represented  in  Boxes  16  through  23.  This  procedure  is  described  more 
fully  below. 

When  no  return  is  received,  for  one  reason  or  another,  the  logical 
flow  branches  to  the  column  headed  by  the  symbol  @ (Boxes  32  through 
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Figure  3.7.  Track  Flow  Chart 
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Figure  3,7  (Cont.) 
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38).  This  Is  Che  "no-return"  mode.  It  decides  whether  to  define  a 
"lost  target"  state*  then  continues  to  attempt  track  as  previously 
described. 

3. 3. 6. 2 Interactions  with  the  Tracking  Filter 

Data  Set  Pointer  DSPWRD.  DSPWRD,  set  in  Box  16,  is  the  DSP  word 
for  the  track  file  into  which  the  current  measurement  is  to  be  incorporated. 
It  must  not  be  zero,  since  the  FILTER  routine  will  assume  this  to  be  an 
initialization  call  if  zero. 

Integer  M (Box  16) . The  Integer  M is  the  dimensionality  of  the 
measurement  vector;  it  is  always  3 for  radar  measurements  in  the  present 
ROSCOE. 


Vector  YMES  (Box  16)  . Vector  YMES  is  the  three-vector  giving  the 

current  radar  measurement  in  radar  face  coordinates  (R  ,u  ,v  ). 

m m m 

Vector  SIGMA  (Box  16) . Vector  SIGMA  is  the  three-vector  giving 
the  estimated  radar  measurement  errors  (a  ,a  ,a  ) associated  with  the 

K U V 

current  measurement.  These  are  computed  (in  Box  12)  in  radar  face 
coordinates  and  stored  in  SIGMA. 

r 

j • ' 

First  Call  to  Subroutine  FILTER  (Box  20).  After  setting  up  the 
input  variables  as  described  above,  and  setting  the  value  of  t , sub- 
routine FILTER  is  called.  It  returns  an  "expected  state  vector"  CARTEZ 
and  an  Integer  KFLAG,  which  are  used  as  described  below. 

Integer  KFLAG  (Box  22a) . The  integer  KFLAG  indicates  when  the 
FILTER  subroutine  is  finished  absorbing  the  current  measurement.  As 
long  as  it  is  returned  by  the  subroutine  with  the  value  zero,  a new 
Iteration  of  the  "expected"  state  vector  (CARTEZ)  is  being  returned  to 
be  converted  into  radar  face  coordinates  and  fed  back  into  FILTER.  When 
the  filter  is  satisfied  with  the  current  iteration  it  feeds  back 
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KFLAG  - 1 , whereupon  the  current  estimated  state  vector  is  ready  to 
be  put  into  the  track  file  (Box  23).  Several  iterations  may  be  required. 

Vector  CARTEZ  (Box  21) . Vector  CARTEZ,  a standard  10-dimensional 
state  vector  in  G.R.C.  coordinates,  is  returned  by  subroutine  FILTER. 

It  represents  the  current  iteration's  estimate  of  the  target  state  vector, 
which  must  be  compared  with  the  current  radar  measurement  in  order  to 
perform  the  next  Iteration.  Since  the  FILTER  subroutine  does  not  know 
the  position  of  the  radar,  this  vector  must  be  converted  by  the  radar 
module  into  a form  matching  that  of  the  measurements.  This  is  done  by 
letting 

W(I)  - CARTEZ(I  +1),  I - 1,2,3 

be  the  expected  target  position  vector  in  G.R.C.  coordinates,  and  con- 
verting to  radar  face  coordinates.  The  resulting  radar  face  coordinates 
(R,u,v)  are  stored  in  vector  YMES,  and  subroutine  FILTER  is  called  again 
(Box  22). 

3. A DISCRIMINATION  EVENT 
3.4.1  Assumptions 

1.  Only  one  radar — the  same  one  that  was  used  for  search  and 
track  initiation — is  used  for  discrimination. 

2.  The  initial  discrimination  event  is  created  by  the  track 
initiation  event  and  thus  follows  it  in  the  time  ordered 
processing  of  events. 

3.  Subsequent  discrimination  events  are  created  in  Subroutine 
DISCRIM  (see  flow  diagrams)  at  time  intervals  t^  over  a 
period  Tid  (input  quantities) . 

4.  The  discrimination  process  will  proceed  independently  of 
the  other  radar  events  once  it  has  been  initialized. 
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5.  Here,  as  In  the  other  radar  events,  the  transmission  time 
of  the  signal  Is  ignored;  the  transmission  and  reception 
of  the  signal  occur  at  the  event  time. 

3.4.2  Procedure 

Two  types  of  discriminants  are  available  for  simulation  in  ROSCOE: 
Fine-Frequency  Length  and  Wide-Bandwldth  Length. 

The  Fine-Frequency  Length  discriminant  involves  making  measurements 
of  the  target's  body  length  by  using  a special  waveform.  For  reasons 
which  will  not  be  discussed  here,  this  measurement  cannot  be  made  if  the 
target  is  below  a specified  altitude.  While  the  actual  limiting  altitude 
depends  somewhat  on  the  nature  of  the  target,  which  cannot  be  assumed 
to  be  known,  it  is  possible  to  define  an  altitude  Hpp  above  which  the 
measurement  can  be  made  for  all  targets.  This  quantity  is  an  input  to 
the  program. 

The  DISCRIM  routine  first  checks  whether  the  user  has  specified  the 
use  of  this  discriminant.  If  so,  the  target  state  vector  is  projected 
ahead  to  the  current  time  for  this  event  and  the  target  altitude,  H , 
is  computed.  Then,  if  H > Hpp  , the  Fine-Frequency  Length  discrimination 
is  performed.  If  H <_  Hpp  , this  measurement  cannot  be  made. 

Subroutine  DISCRIM  then  checks  whether  the  user  has  specified  the 
use  of  the  Wide-Bandwldth  Length  measurement  (independent  of  whether  a 
Fine-Frequency  Length  measurement  has  been  made).  This  discriminant  also 
requires  a special  waveform,  and  cannot  be  performed  if  the  target  is  too 
far  away.  Once  again,  although  the  limiting  range  depends  somewhat  on 
the  target  and  radar  characteristics,  a value  of  R^g  can  be  defined  in 
the  input  such  that  a measurement  can  be  made  whenever  the  target's  slant 
range  is  less  than  R^g  . If  the  user  has  specified  this  discriminant, 
the  target  range  at  the  current  time  is  computed  and  compared  with  R^g  . 
If  R < R^g  , the  Wide-Bandwldth  Length  calculation  is  made  and  a summary 
output  dataset  created  and  put  on  the  output  list. 
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Finally,  the  total  elapsed  discrimination  time  (T^)  is  computed  and 
checked  against  the  Input  time  limit  (T^) • If  the  elapsed  time  has  not 
exceeded  the  limit , a new  discrimination  event  is  created  at  a time 
interval  tD  in  the  future.  If  the  limit  has  been  reached,  no  new 
discrimination  events  will  be  set  up. 

Figure  3.8  is  a flow  diagram  of  the  discrimination  event.  The 
following  subsections  describe  the  specific  models  used  to  implement  the 
Fine-Frequency  Length  and  Wide-Bandwidth  Length  discriminants. 

3. A. 3 Fine-Frequency  Length  Model 

If  we  imagine  a simple  "dumbbell"  target  with  point  scatterers  of 
cross  section  a ^ and  at  the  front  and  rear,  then  the  total  radar 

cross  section  is  given  by 


o - o,  + o + 2 \la ,a  cosf-^— ^ Llcos  0 1 ) 
net  f r Vfr  \ c ' '/ 

where  L is  the  length  of  the  body,  8 is  the  look  angle,  f is 
the  radar  frequency,  and  c is  the  speed  of  light. 

If  the  cross  section  is  measured  at  four  frequencies,  f^,  f^,  f^, 
f^  , such  that  f^  • f^  + Af,  f^  ■ f^  + 2Af,  and  f^  ■ f^  + 3Af  , then 
one  can  calculate  the  quantity 


°A  ~ C1 
o3  - a2 


sin  3 


27TAfLlcOB  el 
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! f2  7T  A f L | 

sin  L 

r, 

— 


cos  6 


which  reduces  to 
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Figure  3.8.  Discrimination  Flow  Chart 
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This  can  be  solved  for  the  body  length, 


L 


c 

2nAf | cos  0 


L : 

1/2 


3.4.4  Wide -Bandwidth  Length  Model 

If  a body  Is  illuminated  by  a short  pulBe,  the  amplitude  of  the 
reflected  signal,  f(t),  can  provide  a measure  of  the  length.  Specifi- 
cally, twice  the  standard  deviation  of  f(t)  , i.e. , 


is  called  the  pulse  stretching  in  this  method,  and  the  length  xs  approxi- 
mated as  cT/2. 


As  an  example  consider  a target  with  front  and  rear  point  scatterers 
of  cross  sections  o^  and  o^  , separated  by  a distance  L . If  the 
inherent  spread  of  the  radar  pulse  is  such  that  the  return  from  a perfect 
point  scatterer  is 


f(t)  « e 


-t2/2 a2 


then  the  above  approximation  will  give  a "length" 


"L"  -Lx 


2 (o.o  ) 
f r 


1/4 


/jT  + 

f r 


(^7  + srf  , ,2 
1 + l r/  ^ 


1/2 


4/o,a 
f r 
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The  true  length  L is  modified  by  two  factors,  the  first  an  error 
caused  by  the  fact  that  the  two  scatterers  have  different  cross  sections 
and  the  second  caused  by  the  non-zero  variance  of  the  radar  pulse. 


If  the  radar  waveform  has  no  inherent  spread  (i.e.,  the  variance 

2 

o * 0) , the  square  root  term  goes  away,  but  the  estimated  L is  still 
slightly  in  error  by  a factor 


2*af0r^4 

Jo~c  + </a~ 
f r 


If  the  two  scatterers  are  equal,  this  factor  is  unity  and  there  is 
no  error.  If  they  are  unequal  Jt  differs  from  1,  but  only  very  weakly. 

If  one  scatterer  has  a cross  section  10  dB  larger  than  the  other  the 
error  is  only  15%;  Indeed,  if  one  cross  section  is  30  times  the  other 
(15  dB  difference),  the  error  is  only  about  28%.  It  is  doubtful  whether 

cross  sections  with  that  great  a dB  difference  could  be  measured  anyway. 

Because  of  the  finite  bandwidth  of  the  radar  signal  the  variance 
is  not  zero,  but  as  long  as  the  Inherent  width  of  the  pulse  (co)  is  much 
less  than  the  target  spread  (L)  , the  term  in  the  square  root  is  negligible. 
If  plasma  dispersion  Increases  the  variance  so  thAt  co  is  on  the  same 
order  as  L (or  greater) , the  measured  length  "L"  depends  on  the  length 
of  the  pulse  rather  than  the  target  characteristics;  if  ac  >>  L , then 

"L"  - ca  . In  this  case,  a short  decoy  would  appear  to  be  as  long  as 

an  RV. 

3.4.5  Flow  Diagrams 

The  discrimination  models  described  above  have  been  coded  into  two 
sub.,  utines,  FFLD  and  WBLD.  Flow  diagrams  for  these  routines  are  given 
in  Figs.  3.9  and  3.10. 


66 


Figure  3.9.  Fine-Frequency  Length  Discrimination  Flow  Chart 
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Figure  3.10.  Wide-Bandwidth  Length  Discrimination  Flow  Chart 
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RADAR  SIGNAL  PROCESSING 


4.1  OVERALL  STRUCTURE 

The  signal  processing  applied  to  each  radar  pulse  (radar  look  event) 

Is  simulated  by  a group  of  subroutines  called  by  a controlling  subroutine, 
RADMODS.  The  order  In  which  these  subroutines  are  called  is  designed  to 
minimize  the  required  computation  by  calculating  the  simplest  effects  first. 
For  example,  a target  that  cannot  be  detected  because  of  absorption  is  de- 
leted from  the  list  for  further  processing,  so  that  neither  the  other 
environmental  effects  nor  the  effects  of  radar  characteristics  need  be  cal- 
culated for  that  target. 

Figure  4.1  shows  the  hierarchy  of  subroutines  called  by  RADMODS.  The 
processing  begins  when  RADMODS  is  called  by  one  of  the  radar  event  sub- 
routines (Sec.  3)  to  calculate  the  apparent  position  and  signal-to-noise 
ratio  for  a particular  target.  RADMODS  calls  subroutines  to  carry  out 
these  calculations  in  four  steps: 

1.  Targets  that  might  interfere  with  the  desired  target  are 
listed,  taking  into  account  the  possible  effects  of  refrac- 
tion. 

2.  The  effects  of  absorption,  noise,  refraction,  dispersion, 
and  Faraday  rotation  are  computed  for  each  target  on  the 
list. 

3.  A revised  target  list  is  prepared,  taking  into  account  the 
foregoing  environmental  effects  as  well  as  any  false  targets 
due  to  reflection  from  fireballs. 

4.  For  each  target  on  the  revised  list,  the  apparent  position 
R* , u',  v'  and  the  sign^l-to-noise-plus-clutter  ratio  ZQ 
are  calculated,  taking  into  account  the  radar  monopulse  pro- 
cessing and  range  gating  and  the  effects  of  clutter.  As 
explained  earlier,  the  apparent  position  R',  u',  v'  does 
not  include  radar  measurement  errors. 
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Figure  4.1.  Radar  Signal  Processing  Subroutine  Hierarchy 
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AN-4/759A 


RADMODS  then  returns  to  the  calling  event  with  the  values  R* , u’, 
v',  and  ZD  . Two  further  subroutines  are  called  by  the  event  rather  than 
by  RADMODS:  XTHRSHS  compares  ZG  with  the  threshold  signal-to-noise-plus- 

clutter  ratio  specified  for  that  event  type,  and  MEASERR  adds  radar  measure- 
ment errors  (bias  and  random)  to  R' , u' , v'  to  give  the  measured  target 

position  R , u , v 
m m m 

Figure  4.2  is  a flow  chart  of  subroutine  RADMODS,  showing  the  sequence 
in  which  it  calls  the  other  subroutines.  The  first  decision  block  in  Fig. 
4.2  chooses  subroutine  POSSV  or  subroutine  POSLIST  to  construct  the  list  of 
possx  'le  targets,  POSLST,  depending  on  whether  the  event  type  is  search/ 
verify  (no  range  gating)  or  track  initiate/track  (with  range  gating).  In 
the  latter  case,  it  may  happen  that  there  is  no  target  (not  even  the  desired 
one)  that  might  be  in  the  range  gate;  in  that  case,  the  length  of  POSLST  is 
zero  (NLST  ■ 0) , and  RADMODS  sets  the  "no  detection"  flag  (IFDET  ■ 0)  and 
returns. 

Subroutines  RADIS  and  REF1S  are  next  called  to  compute  the  environ- 
mental effects;  in  each  case,  if  none  of  the  targets  on  the  list  can  be 
detected,  a flag  IFAIL  is  set  to  a value  other  than  zero  and  the  "no  detec- 
tion" return  is  made. 

The  revised  target  list  is  prepared  in  two  different  ways  (LISTA  or 
LISTB) , depending  on  whether  there  is  a single  (direct  or  refracted)  path 
from  radar  to  target  (Np  ■ 1)  or  many.  In  the  former  case,  unless  the  re- 
vised list  is  empty  (MLST  ■ 0),  the  subroutines  to  compute  R',  u',  v', 
and  Z0  are  called,  the  choice  of  subroutine  depending  on  whether  there 
is  more  than  one  target  on  the  revised  list.  In  the  latter  case,  the  mul- 
tiple refraction  is  assumed  to  produce  "spread  target  clouds"  of  images, 
which  are  treated  statistically  in  computing  R',  u’,  v' , and  Z0  . 

Inputs:  an  indicator  of  the  calling  event  type 

Outputs:  IFDET,  (R' , u',  v'),  ZQ  * 

A 

Because  many  of  the  signal  processing  subroutines  have  numerous  inputs  and 
outputs,  they  are  listed  in  symbolic  form  only.  See  Appendix  G for  defi- 
nitions of  symbols  used  here  and  in  flow  charts. 
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Figure  4.2.  Subroutine  RADMODS 


72 
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Before  going  on  to  describe  the  processing  steps  in  detail,  a few 
comments  ought  to  be  made  concerning  features  that  are  not  simulated  (some 
of  which  have  been  mentioned  earlier) . 

• Only  a single  radar  beamshape  is  used:  circular  with  a 

constant  or  tapered  sidelobe  level.  No  beam  stacking  is 
modeled:  transmit  ar.d  receive  beams  have  identical  para- 

meters and  axis  locations. 

• No  monopulse  is  used  in  search  and  verification  events. 

The  angular  position  output  by  these  procedures  is  that 

of  the  center  of  the  radar  beam.  The  range  output  is  that 
of  the  highest  peak  within  the  resolution  cell  centered  on 
the  target  of  interest. 

• For  all  events,  pulse-compression  (chirped)  signals  may  be 

used.  A single  model  of  the  matched  filter  output  (ambiguity 
function)  is  provided  for  pulse-compression  signals:  this 

model  is  an  approximation  falling  between  the  shapes  for 
Hamming-weighted  and  Hanning-weighted  waveforms,  and  has  a 
constant  range  sidelobe  level. 

• Dispersive  distortion  of  the  ambiguity  function  is  taken 
into  account,  but  the  model  used  is  strictly  correct  only 
for  nnchlrped  waveforms.  For  pulse-compression  signals, 
the  uodel  used  is  based  on  a heuristic  approach. 

• A single  range  marking  system  is  provided  for  events  other 

than  search  and  verification:  the  split-gate  range  tracker. 

• It  is  assumed  that  when  refraction  is  present,  only  two  cases 

exist:  either  there  is  a single  ray  path  between  the  radar 

and  the  target,  or  there  are  very  many  direct  paths  (striated 
debris  case)  producing  a "spread  target  cloud"  of  many  images. 
The  former  c»se  is  treated  deterministically,  the  latter  case 
stochastically.  No  intermediate  case  (e.g.,  three  direct 
paths)  is  provided  in  the  present  version  of  ROSCOE. 


• The  method  of  dealing  with  monopulse  and  range-tracking  in 
the  "spread  target"  case  is  a simple  random  phase  approxi- 
mation, which  is  considered  adequate  for  most  practical 
cases. 

• False  alarms  are  not  modeled. 

4.2  POSSIBLE-TARGETS  LIST 

4.2.1  Subroutine  POSSV  (Fig.  4.3) 

This  subroutine  performs  two  tasks  needed  to  initiate  the  signal  pro- 
cessing for  search  and  verification  events.  It  is  used  only  with  those 
events,  since  it  makes  certain  assumptions  which  are  only  correct  for  them: 
that  there  is  one  particular  (specified)  target  of  concern  to  the  event  at 
this  time,  that  the  event  is  not  interested  in  detecting  any  other  targets 
except  those  which  might  cause  direct  interference  with  the  given  target, 
and  that  no  range  gate  is  defined  or  used  by  this  event. 

Subroutine  POSSV  first  finds  the  current  parameters  and  true  state 
vector  of  the  given  target,  and  determines  itB  doppler  shift,  fd  . It 
then  finds  the  value  of  fd  , the  center  frequency  of  the  doppler  "bin"  in 
which  the  return  from  this  target  will  be  found;  here  it  is  assumed  that 
for  search/verify  functions  the  receiver  has  an  evenly  spaced  bank  of  doppler 
filters  so  that  all  targets  can  be  detected  even  though  no  a priori  informa- 
tion on  their  velocity  is  available;  these  doppler  bins  will  be  sufficiently 
narrow  and  numerous  that  any  reasonable  target  velocity  will  cause  the  target 
to  appear  in  one  of  the  bins  with  no  more  than  1 dB  of  doppler  mismatch  loss. 

The  subroutine  then  goes  on  to  prepare  a list  of  all  the  real  targets 
which  might  be  of  concern  to  the  present  event  and  beam  position:  the  list 

[called  POSLSTl  developed  by  this  subroutine  contains  only  the  given  target 
of  Interest  plus  any  other  targets  which  might  later  turn  out  to  be  within 
the  radar  beam  and  in  the  same  resolution  cell.  The  region  of  interest  is 
conservatively  defined  with  the  possibility  of  later  refractive  effects  in 
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Figure  4.3.  Subroutine  POSSV 
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mind,  by  including  all  targets  within  AJ{  half-beamwidths  of  the  desired 
target.  Currently  A ■ 4 . 

Inputs:  t,  (u  ,v  ),  <5w,  6R  , NT  , B 

v * o o o’  o n 

Outputs:  ?d,  NLST,  POSLST 

4.2.2  Subroutine  POSLIST  (Fig.  4.4) 

This  subroutine  prepares  a list  of  those  real  targets  which  might  be 
of  concern  to  the  present  event,  beam  position,  and  range  gate.  It  is  called 
by  the  events  (track  initiation  and  track)  which  specify  a range  gate  center 
and  gate  width  in  addition  to  the  beam  position.  The  list  (also  called 
POSLST)  developed  by  this  subroutine  contains  all  real  targets  which  might 
later  turn  out  to  lie  within  the  radar  beam  and  within  one  range  resolution 
distance  on  either  side  of  the  range  gate  limits.  This  region  may  or  may 
not  contain  the  target  of  particular  concern  to  the  present  event  at  the 
given  instant.  The  region  is  conservatively  defined  with  the  possibility 
of  later  refractive  effects  in  mind,  by  using  the  constant  defined  for 

POSSV  and  a range  constant  (currently  3000  m). 

Inputs:  t,  (Rc,uo,vo),  W,  Sw,  <5Rq,  NTo 

Outputs:  NLST,  POSLST 

4.2.3  List  POSLST 

This  list  is  created  (before  refraction  calculations)  in  either  POSSV 
or  POSLIST,  depending  on  the  calling  event.  The  list  contains  NLST  entries 
(if  NLST  » 0 there  is  no  list),  each  of  which  is  a dataset  pointer  (DSP) 
to  a dataset  containing  the  following  parameters  for  one  target: 

IT  serial  number  of  thiB  target  in  the  list  (1  <_  IT  <_  NLST) 

NT  corresponding  target  identification 

R,u,v  current  true  coordinates  of  the  target 

o current  target  cross  section 

f^  current  target  doppler  shift 
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Figure  4.4.  Subroutine  POSLIST 
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£ current  true  angle  off  the  beam  axis  for  this  target  In 
sine-space: 

£2  = (u  - uq)2  + (v  - VQ)2 

4 

X range-weighted  cross  section,  o/R 

IFIT  flag  to  indicate  whether  this  is  (IFIT  ■ 1)  or  is  not 

(IFIT  * 0)  the  target  of  primary  interest  to  this  event 

4 . 3 ENVIRONMENTAL  EFFECTS 

4.3.1  Subroutine  RADIS  (Fig.  4.5) 

This  subroutine  begins  the  job  of  calculating  the  radar  return. 
Selecting  an  appropriate  target  from  among  those  in  POSLST,  it  conserva- 
tively estimates  the  maximum  amount  of  absorption  along  the  path  to  this 
point,  and  the  maximum  amount  of  noise  in  the  beam  position,  which  could 
exist  and  still  allow  a threshold  crossing  to  be  possible.  The  actual  ab- 
sorption and  noise  are  then  computed  by  calling  the  propagation  subroutines 
ABSORB  and  NOISE  (see  Sec.  5).  If  these  exceed  the  estimated  maxima,  the 
computation  is  aborted  with  a message  to  indicate  the  reason  for  no  return. 
If  the  absorption  and  noise  are  less  than  these  limits,  their  values  are 
recorded  and  certain  parameters  of  the  radar  equation  are  calculated  and 
output. 

The  "conservative  estimate"  of  maximum  allowable  noise  and  absorption 
Is  accomplished  by  multiplying  the  target  cross  section  by  a factor  a . 
This  factor  is  a predesignated  constant;  the  value  a - 10  is  presently 
specified  in  a data  statement  in  subroutine  RADIS.  This  allows  for  focus- 
ing and  other  enhancement  effects. 

This  subroutine  operates  in  terms  of  a form  of  the  radar  equation 
which  calculates  the  pre-detection  signal-to-noise-plus-clutter  ratio, 

, in  terms  of  four  input  radar  parameters  and  several  factors  computed 
in  this  subroutine  and  others : 
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Figure  4.5.  Subroutine  RADIS 


where 


- (s/n)t1(rq/r) 


(cos  0)FtFr0 


L [1  + (T  /T  ) + (C/N  )] 
o x n n 


* 


(S/N) 


Ti 


pre-detection  signal-to-noise  threshold 

2 

range  at  which  aim  target  will  just  exceed 
this  threshold  under  nominal  conditions** 


R * current  range  to  target 
0 ■ off-boresight  angle  of  beam  axis 

F^,  ■ relative  gain  factor  due  to  target's  position 

off-axis  of  transmit  beam 

F ■ relative  gain  factor  due  to  target's  position 
off-axis  of  receive  beam  (may  differ  from  F-p 
because  transmit  and  receive  paths  to  the  target 
need  not  be  the  same) 

o - target  cross  section 

t 

■ nominal  system  noise  temperature 

Nn  - nominal  system  noise  power  (kB^^) 

i 

■ noise  bandwidth  for  the  transmitted  signal 
(Bn  - 1/pulse  length) 

Lq  * total  outside  loss  factor,  including  absorption 
(L2),  dispersive  loss  (Lp) , and  Faraday  rotation 
loss  (Lp) , if  applicable 

Tx  ■ total  effective  outside  noise  temperature 
C - total  effective  clutter  power 


# 


» 


Factors  marked  by  an  asterisk  are  ir put  radar  parameters. 

Target  on  the  beam  axis  at  radar  boresight,  and  no  outside  losses,  noise, 
or  clutter. 
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For  some  applications,  it  is  convenient  to  exhibit  directly  the  de- 
composition of  into  the  ratio  of  a signal  power,  S , to  a noise-plus- 

clutter  power,  Xp  : 


(cos  9)FtFr° 


L R 
o 


3C  - Nil  + (T/T  ) + (C/N  )] 
in  x n n 


where 


<S/N)n  Ro»n 


is  computed  in  advance  for  each  radar  waveform. 


IFA1L  la  an  index  to  indicate  whether  a threshold  crossing  may  be 
possible,  and  if  not,  why  not. 


Value  of  IFAIL 

0 

1 

2 

3 

4 

5 


Meaning 

Possible 

Impossible-  -range  too  great 

Impossible — too  much  absorption 

Impossible — too  much  outside  noise 

Impossible — absorption  plus  noise 

Impossible — dispersive  and  Faraday 
rotation  losses 

Impossible — no  targets  in  LISTA  (Sec.  4.4.1) 


Inputs:  NLST,  POSLST,  (u  ,v  ),  6w,  R , R , S , T , N , a 

■ o o o c o n n 

Outputs ! IFAIL,  ITO,  cos  0,  -vt*  ^2 


4 - 
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4.3.2  Subroutine  REF1S  (Fig.  4.6) 

This  subroutine  calculates  further  propagation  effects:  refraction, 

dispersion,  and  Faraday  rotation. 


To  begin  with,  the  user  may  define  a limiting  altitude,  , below 

which  no  refractive  effects  are  to  be  computed.  The  idea  is  that  for  ray 
paths  lying  entirely  below  Href  » significant  refraction  will  always  be 
accompanied  by  overwhelming  absorption,  so  it  would  be  a waste  of  time  to 
calculate  the  refraction.  Note  that  both  the  target  and  the  radar  must  lie 
below  in  order  for  this  condition  to  be  applicable.  The  user  may, 

of  course,  select  Href  ■ 0 , in  which  case  this  test  becomes  irrelevant. 


Refraction.  It  is  assumed  in  ROSCOE  that  the  number  of  propagation 
paths  between  the  radar  and  the  target  (Np)  is  either  one  or  much  greater 
than  one.  The  latter  case  occurs  only  when  there  is  an  intervening  striated 


region.  When 
N2  . 


Np  >>  1 , the  number  of  apparent  target  images  is  given  by 


Subroutine  REF1S  computes  only  one  set  of  refraction  parameters: 
those  affecting  the  previously-defined  target  whose  coordinates  are  (R^, 
Ul,Vl).  These  aame  parameters  are  then  issumed  to  apply  unchanged  to  all 
targets  within  the  region  of  interest.  In  other  words,  the  range  and  angular 
extent  of  the  region  is  assumed  to  be  small  enough  so  that  the  same  refrac- 
tion affects  all  targets  within  it.  This  does  not  necessarily  mean  that  all 
these  targets  suffer  the  same  refractive  displacements,  since  there  may  be 
uncertainties  associated  with  each  bias  component. 

Subroutine  REFRCT  (Sec.  5.3)  is  called  to  compute  the  refraction  para- 
meters. It  first  determines  whether  there  is  a striated  region  between  the 
radar  and  the  target  point  (R^,u^,v^).  If  so,  it  estimates  the  number  of 
propagation  paths,  Np  . If  not,  it  sets  Np  ■ 1 . In  either  case,  it 
computes  the  integrated  electron  density, 
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CALL  SUB  “REFRCT"  TO  COMPUTE 
SIMPLIFIED  REFRACTION  ERRORS, 
NUMBER  OF  PROPAGATION  PATHS, 
Np,  AND  INTEGRATED  ELECTRON 
DENSITY.  IE*  TO  TARGET  AT 

(1*1  > V V,) 


lf 

* i 

< 

COMPUTE  DISPERSION 
PARAMETER:  , ,, 

A c P A T * " 

aD  fPc  *Dn  ‘E 

CALL  SUB  "DISPERS"  TO 
COMPUTE  SIMPLIFIED 
DISPERSIVE  LOSS,  Lp 


/ IS  N.  MO 

< POLARIZATION  >_ 

\LIMEAR  / 

JUS 


DETERMINE  THE  MEAN  ANGLE, 

V,  BETWEEN  THE  RADAR  LINE- 
OF-SIGHT  AND  THE  GEOMAGNETIC 
FIELD  IN  THE  REGION  OF 
GREATEST  N_  ALONG  PATH 


COMPUTE  THE  FARADAY 
ROTATION  ANGLE, 

HF  * 104  I£  COS  F/f2. 

AND  LOSS, 

2 

Lp  * sec  «p 


Figure  4.6.  Subroutine  REF1S 


I_  * / N ds 

E J e 

along  the  direct  path  between  the  radar  and  this  target.  This  value  of 
1 j7  (conventionally  given  In  units  of  electrons /cm^)  is  assumed  in  this 
version  to  apply  to  all  targets  in  the  region  of  interest. 

When  Np  - 1 , the  properties  of  the  refraction  are  given  by  six 
parameters : 

AR,  AA,  AE,  aAR,  aAk,  o^g 

The  first  three  parameters  give  the  refractive  bias  errors  which  must  be 
added  to  the  true  target  coordinates  (R,A,E)  to  give  the  apparent  target 
coordinates  (R',A' ,E').  The  last  three  parameters  give  the  uncertainties 
(if  any)  of  these  bias  errors,  in  terms  of  their  standard  deviations.  The 
use  of  these  parameters  is  described  in  Sec.  4.4.1,  Subroutine  REFLST1. 

When  Np  > 1 (striated-debris  refraction),  there  are  many  false 
target  images  scattered  over  a three-dimensional  region  of  space.  In  this 
case,  subroutine  REFRCT  provides  statistical  information  on  the  spatial 
distribution  of  these  images,  and  on  their  effective  cross-section  distri- 
bution. 

2 

The  spatial  distribution  of  the  set  of  target  images  is  assumed 

to  have  the  form  of  a probability  ellipsoid.  The  dimensions  and  orienta- 
tion of  this  ellipsoid  are  approximately  determined  from  the  phenomenology. 
For  the  target  whose  true  position  is  (R^,u^,v^),  the  centroid  of  this 
ellipsoid  is  located  at  (R^  + AR,  uj_  + Au,  v^  + Av);  the  values  of  AR,Au,Av 
are  among  the  outputs  of  REFRCT.  The  3x3  covariance  matrix  defining  the 
error  ellipsoid  in  Cartesian  space  is  projected  onto  the  radar  line-of-sight 
to  produce  the  range-error  standard  deviation  . The  matrix  is  also 

projected  onto  the  sine-space  plane,  using  a linearized  approximation  based 
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on  Che  assumption  chat  Che  dimensions  of  Che  ellipsoid  are  absolutely  small; 
the  result  is  an  error  ellipse  in  sine-space,  defined  by  its  semimajor  axis 
(ai),  its  semiminor  axis  (02),  and  Che  angle  (y)  between  Che  major  axis  and 
Che  u-axls  of  Che  sine-space  coordinate  system  (to  eliminate  possible  am- 
biguities, the  model  selects  0 y < it). 

The  effective  cross  sections  of  the  target  images  will  be  affected  by 
focusing  and  defocusing,  so  that  they  differ  from  the  true  cross  section  of 
the  real  target  and  in  fact  will  have  a distribution  of  values.  The  sub- 
routine provides  the  mean  (c)  and  standard  deviation  (o^)  of  this  cross- 
section  distribution  under  the  assumption  that  it  is  Gaussian. 

Thus,  for  Np  > 1 , the  output  of  subroutine  REFRCT  is  the  following 
set  of  parameters: 

AR,  Au,  Av,  oAR,  olf  02,  y,  0,  Oq 

In  using  these  parameters,  ROSCOE  makes  the  following  simplifying 
assumptions : 

1.  The  errors  in  R,u,v  are  uncorrelated 

2.  The  standard  deviations  are  all  sufficiently  small  that 
linearization  is  permissible  and  sine-space  coordinates 
can  be  treated  additively. 

3.  The  random  phase  approximation  can  be  employed  when  com- 
bining the  signal  returns  from  this  set  of  false  target 
Images . 

Dispersion.  This  is  treated  by  a method  derived  from  an  analysis 
which  is  described  in  Appendix  B. 

The  dispersion  parameter,  Ap  , for  a signal  of  bandwidth  B , is 
defined  to  be 
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where 


Ajj  “ 4yB 


2 ■ -4  k - <^2>r2  - 

C<0 


and  where 


to 


• 2irf 


wf  = 3.182  x 10*N  rad/s 

N e 

In  most  practical  cases,  y can  be  approximated  in  terms  of  the  In- 
tegrated electron  density,  I£  , by  the  formula 


2 

y 


4.28  x 10"A  IE/f 3 


where  f Is  the  radar  frequency  In  Hz.  Therefore,  we  write 

-2  3 1/2 

Ap  * 8.28  x 10  * (I  /r)  B 

This  formula  was  actually  derived  only  for  rectangular  unmodulated  pulses, 

where  B = B “ 1/x  . For  these  we  can  write 
n 


^ ~ ADnIE 


1/2 


where 


A a 0.0828  B f 
un  n 


-3/2 


For  chirped  pulses,  the  effective  bandwidth  is 


B = F B 
pc  n 


where  Fpc  is  the  pulse-compression  ratio.  Therefore,  for  such  signals 
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we  use  for  the  dispersion  parameter  the  form 


The  parameter  A^n  is  precomputed  for  each  waveform  and  stored  in  the  ap- 
propriate dataset. 

In  terms  of  this  dispersion  parameter,  the  dispersive  mismatch  loss 
(Lq)  at  the  peak  of  the  output  of  the  receiver  filter  is  given  by  a com- 
plicated expression  Involving  Fresnel  integrals  and  trigonometric  functions. 
Since  this  Is  only  approximately  valid,  however,  especially  in  the  case  of 
chirped  signals,  subroutine  DISPERS  (Sec.  4.3.3)  computes  an  approximation 
to  Lq  by  using  a piecewise  c.  rve  fit  to  a plot  of  the  complicated  expres- 
sion mentioned  above. 

Faraday  Rotation.  If  the  radar  uses  a linearly  polarized  signal, 
Faraday  rotation  has  to  be  taken  into  account.  Because  of  this  rotation, 
some  of  the  return  signal  may  arrive  in  the  crossed  polarization  and  hence 
not  contribute  to  the  received  power;  this  results  in  the  "Faraday  rotation 
loss,"  Lp  . 

These  two  losses,  LD  and  Ly  , are  combined  with  the  other  known 
losses  and  noise  factors  (excluding  clutter,  which  will  be  calculated  later) 
and  compared  with  a conservatively  estimated  value  of  the  return  signal 
power  (as  in  RADIS)  to  see  if  a threshold  crossing  is  still  possible.  If 
not,  the  subroutine  terminates  the  computation  and  returns  with  an  appro- 
priate diagnostic  index. 

Finally,  the  degraded  range  resolution  widths,  6R^  and  (if  chirped 
signals  are  used)  6R|  , are  estimated.  The  estimate  is  made  by  assuming 
that  the  ambiguity  function  in  the  range  dimension  is  triangular  (as  it  is 
for  an  unchirped  signal  in  the  absence  of  dispersion) , with  its  peak  ampli- 
tude lowered  by  and  its  width  correspondingly  widened  so  as  to  satisfy 
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Che  principle  of  energy  conservation  before  and  after  dispersion.  By  this 
method,  the  degraded  range  resolution  turns  out  to  be  simply  Lq  times  the 
nominal  resolution  width.  This  is  a good  approximation  for  an  unmodulated 
rectangular  pulse,  for  which  comparisons  have  been  made  with  the  exact  cal- 
culation. ROSCOE  also  uses  the  same  technique  for  a chirped  pulse,  though 
no  estimates  of  the  goodness  of  the  approximation  have  been  made  for  this 
case. 


The  nominal  range-resolution  width  for  an  unchlrped  rectangular  pulse 
is 


6R  ■ ct/2 
o 


This  is  also  the  half-extent  of  the  region  in  which  the  range  sidelobes  are 
important  for  a chirped  pulse. 


The  nominal  range-resolution  width  for  a chirped  pulse  is 


26R  /F 
o 


pc 


The  numerical  factor  takes  roughly  into  account  the  spreading  which  results 
from  the  weighting  usually  applied  to  reduce  the  magnitude  of  the  range 
sidelobes. 

Inputs:  (Rj^Uj^),  H^.,  H^,  Fpc>f»  6Rq,  SR^ 

polarization  indicator,  S, , L , Y,  IFAIL 

i om 


Outputs: 

V v v 

lf* 

fiF’ 

V 

<5r[,  IFAIL 

Refraction 

JAR, 

AA, 

AE, 

°AR’ 

a , a if  N - 
AA’  AE  p 

1 

Parameters  | 

Iar, 

Au, 

Av, 

— 

if  N 

1 

aAR’ 

0^ , o 2 , P»  of  aQ : 

4.3.3  Subroutine  DISPERS  (Fig.  4.7) 

This  subroutine  computes  the  loss  factor,  Lp  , due  to  dispersive 
effects  along  the  line-of-sight  between  the  radar  and  the  target.  This 
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loss  factor  represents  the  power  loss  at  the  peak  of  the  ambiguity  function, 
and  is  computed  by  a fit  to  a simplified  model  which,  while  strictly  correct 
only  for  unchirped  waveforms,  is  approximately  applicable  to  chirped  wave- 
forms. The  input  to  this  subroutine  is  a dispersion  parameter,  Ap  , which 
depends  on  the  radar  and  signal  parameters  and  on  the  integrated  electron 

i 

density  along  the  line-of-slght.  See  Appendix  B for  the  derivation  of  this 
parameter. 

I 

Input : Ap 

Output : Lp 

A. 4 REVISED  TARGET  LIST 

4.4.1  Subroutine  REFLST1  (Fig.  4.8) 

This  subroutine  is  called  to  prepare  a list  of  targets  which  will 
affect  the  radar  return  in  those  cases  wiiea  there  is  only  a single  refracted 
(or  direct)  ray  path  between  the  radar  and  each  real  target  (Np  - 1).  This 
list,  called  LISTA,  is  constructed  by  applying  the  refractive  parameters 
derived  in  subroutine  REF1S  to  the  targets  in  list  POSLST,  and  making  use 
of  the  known  properties  of  the  radar  waveform,  beam  parameters,  and  receiver 
range  gate  (if  any)  to  determine  which  of  the  target  Images  will  contribute 
to  the  return  for  the  particular  radar  event  and  target  presently  under  con- 
sideration. 

f 

i 

As  shown  in  Fig.  4.8,  the  procedure  for  adding  the  refractive  para- 
meters is  to  select  random  numbers  from  Gaussian  distributions  with  the 
means  (AR,  AA,  AE)  and  standard  deviations  (a^R,  o^,  , which  may  be 

zero)  calculated  by  REFRCT,  and  add  them  to  the  previously  computed  target 
coordinates.  Note  that  conversion  from  sine-space  coordinates  (R,u,v)  to 
angular-space  coordinates  (R,A,E)  is  required  because  angles  do  not  add  in 
sine  space. 

The  apparent  range  of  a target  is  also  affected,  when  a chirp  waveform 
is  used,  by  range-doppler  coupling  (see  Appendix  B).  This  is  accounted  for 
by  the  term  a <f ^ - f^)  in  the  range  expressions  of  Fig.  4.8. 
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The  resulting  refracted  target  coordinates  are  tested  to  see  whether 
they  are  within  a range  resolution  distance  (6R^)  of  the  desired  target 
(or  of  the  range  gate  if  one  Is  used)  and  within  the  25-dB  one-way  radar 
beamwidth  (F^  > 0.003).  (The  off-axis  gain  parameter  is  computed  by 

subroutine  GOA,  Sec.  4. A. 2.)  Targets  within  these  limits  are  placed  on 
LISTA.  Note  that,  when  a range  gate  is  used,  the  desired  target  itself 
may  not  be  on  the  list. 

To  this  list  are  then  added  all  significant  "false  target"  images 

caused  by  signal  reflection  from  low-altitude  fireballs;  the  properties  of 

these  multipath  images  are  computed  in  subroutine  MLTPATH  (Sec.  5).  Once 

these  "false"  images  have  been  inserted  into  LISTA,  all  distinction  between 

"real"  and  "false"  targets  is  thenceforth  disregarded  as  being  irrelevant; 

all  images  are  treated  the  same  way  in  the  succeeding  subroutines.  A final 

step  in  subroutine  REFLST1  indexes  all  target  linages  in  LISTA  into  subgroups 

which  are  mutually  unresolvable  by  the  current  radar  waveform;  this  grouping 

is  not  used  in  the  present  ROSCOE  signal  processing  model,  but  may  be  needed 
* 

in  later  versions. 

Inputs l NLST,  PO S L S T , F^^,  ITO,  (R^,u^,v^), 

{AR,  AA,  AE,  oAR,  o^A,  0AE),  6Rq,  5Rlt  B 

Outputs:  MLST,  NGRPS,  MGRP.  LISTA 

4.4.2  Subroutine  GOA  (Fig.  4.9) 

This  subroutine  computes  the  one-way  gain  factor,  F^  , for  a target 
that  may  be  located  off  the  axis  of  the  radar  beam.  The  factor  is  norma- 
lized to  unity  on  the  axis.  Only  a circular  mainlobe  with  constant  or  tapered 
sldelobes  is  provided  for.  The  calculation  is  done  in  sine-space  coordinates. 


For  example,  it  is  not  inconceivable  that  some  user  may  wish  to  examine 
his  simulated  radar  returns  with  a precision  which  demands  doing  actual 
phasor  sums  on  each  group  of  unresolved  Images. 


Figure  4.9.  Subroutine  GOA 


Besides  , this  subroutine  provides  the  value  of  £ , the  target's  angle 
off  the  beam  axis  In  sine  space. 

Inputs:  (u,v),  (u  ,v  ) , p,  S , D 

o o s s 

Outputs:  F^,  £ 

4.4.3  Subroutine  REFLSTN  (Fig.  4.10) 

This  subroutine  is  called  to  prepare  a list  of  targets  which  will 
affect  the  radar  return  in  those  cases  when  there  are  many  refracted  ray 
paths  between  the  radar  and  each  real  target  (Np  >1).  As  mentioned 
earlier , the  radar  model  assumes  that  this  case  involves  a "spread  target 
cloud"  of  many  Images,  produced  by  an  intervening  striated  debris  region, 
and  these  images  are  defined  and  treated  statistically.  The  list,  called 
LISTB,  is  constructed  by  applying  the  statistical  refractive  parameters 
derived  in  subroutine  REF1S  to  the  targets  in  list  POSLST,  and  making  use 
of  known  properties  of  the  radar  waveform,  beam  parameters,  and  receiver 
range  gate  (if  any)  to  determine  which  of  the  target  image  clouds  will 
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Figure  4.10.  Subroutine  REFLSTN 


9h 


AM  59496 


contribute  to  the  return  for  the  particular  radar  event  and  target  presently 
under  consideration.  Because  there  are  so  many. target  images  that  they  must 
be  dealt  with  statistically  in  this  case,  it  is  assumed  that  a few  additional 
images  produced  by  low-altitude  fireball  multipath  would  be  of  negligible 
importance,  so  these  are  not  included  in  the  calculation  in  this  subroutine. 
"Close-target"  groupings  are  also  not  included.  The  effective  parameters  of 
the  image  clouds  as  viewed  by  the  radar  are  computed  in  subroutine  MULTOAR 
(Sec.  4.4.4). 

Inputs:  NLST.  POSLST,  W,  R , (R, ,u.  ,v.  ),  N , 6r' , o.„,  o,  o 

__  C 111  P OAK  0 

Outputs:  MLST,  LISTB 

4.4.4  Subroutine  MULTOAR  (Fig.  4.11) 

This  subroutine  performs  a similar  function  for  a "spread  target 
cloud"  to  that  which  subroutine  GOA  performs  for  a single  target  point. 

It  combines  the  statistical  parameters  of  the  image  cloud  with  the  beam 
position  and  waveform  parameters  to  compute  the  effective  position  of  the 
weighted  center  of  the  cloud  (R,u,v)  with  respect  to  the  beam,  and  then 
calculates  the  target's  apparent  off-axis  angle  (in  sine-space),  £ , and 
the  effective  one-way  relative  gain  factor,  . 

Inputs:  (R,u,v) , AR,  Au,  Av,  c^,  y,  fj,  Fpc>  (u0.v0>»  P» 

6R1,  B,  f.,  S 
o’  ’ d’  s 

Outputs:  (R,u,v),  £,  F^ 

4.4.5  LISTA  and  LISTB 

These  are  the  lists  of  target  points  which  contribute  to  the  received 
signal.  LISTA  applies  in  the  "single  direct  path"  case  (Np  ■ 1),  and  LISTB 
in  the  "striation"  case  (Np  > 1) . False  target  images  due  to  low-altitude 
fireball  reflections  are  Included  in  LISTA  (but  not  in  LISTB)  in  the  same 
way  as  "real"  targets.  Either  list  contains  MLST  entries  (if  MLST  - 0 
there  is  no  list),  each  of  which  is  a DSP  to  a dataset  containing  the  fol- 
lowing parameters  for  one  apparent  target  (or  target  cloud). 
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Figure  4.11.  Subroutine  MULTOAR 
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4. 4. 5.1  LISTA 
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4. 4. 5. 2 LISTB 
IT 
NT 


R* ,u’ ,v’ 


serial  number  of  this  target  In  the  list  (1  <_  IT  <_  MLST) 

corresponding  target  Identification  (may  be  the  same 
for  several  "targets"  if  some  of  them  are  reflected 
Images  of  the  same  real  target) 

current  apparent  coordinates  of  the  target 

current  effective  target  cross  section 

current  target  doppler  shift 

current  apparent  angle  off  the  beam  axis  for  this 

target  in  sine  space:  C'  ■ [ (u * - u )2  + (v'  - v )2] 

o o 

one-way  relative  gain  off-axis  for  this  apparent  target 

2 4 

cross-section  parameter  for  this  target,  oF^/(R') 

flag  to  indicate  whether  this  is  (IFIT  - 1)  or  is  not 
(IFIT  - 0)  the  target  of  primary  interest  to  this  event 

index  number  which  is  the  same  for  all  targets  belong- 
ing to  a given  dose-target  group 


serial  number  of  this  "target"  in  ^he  list  (1  IT  <_  MLST) 

identification  of  the  real  target  of  which  this  cloud 
is  the  refracted  image 

current  apparent  coordinates  of  the  effective  center 
of  this  cloud  (R'  ■ R,u'  ■ u,v’  ■ v) 

effective  cross-section  of  this  entire  cloud  as  calcu- 
lated using  the  random  phase  approximation 

current  effective  doppler  shift  of  the  effective 
center  of  this  cloud 


i 


I 

\ 
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current  angle  off  the  beam  axis  in  sine-space  for  the 
effective  center  of  this  cloud  (£’  = O 

F^  effective  one-way  relative  gain  off-axis  for  this  cloud 

2 4 

X cross-section  parameter  for  this  cloud,  oF^/(R') 

IFIT  flag  to  indicate  whether  this  cloud  is  (IFIT  ■ 1)  or 
is  not  (IFIT  ■ 0)  the  image  of  the  target  of  primary 
interest  to  this  event 

4.5  TARGET  POSITION  AND  SIGNAL-TO-NOISE-PLUS-CLUTTER  RATIO 

The  final  portion  of  the  radar  signal  processing  calculations  con- 
trolled by  subroutine  RADMODS  is  to  calculate  the  "appa<  ■nt"  target  coor- 
dinates (R'.u'.v*)  and  output  signal-to-noise-plus-clutter  ratio  (ZQ), 
using  the  list  of  interfering  targets  and  their  properties  (LISTA  or  LISTB) 
that  has  resulted  from  the  previous  steps. 

One  of  four  subroutines  is  used,  depending  on  the  situation: 

TARGS1  Np  « 1,  only  one  target  on  LISTA 

TARGM1S  N ■ 1,  more  than  one  target  on  LISTA 

p ° 

TARGMSV  Np  > 1,  no  range  gating 

TARGMTS  Np  > 1,  range  gated  signals 

These  subroutines  in  turn  call  other  subroutines  (Sec.  4.6)  to  cal- 
culate the  effects  of  monopulse,  range  gating,  and  targets  in  the  range 
sidelobes,  and  the  propagation  routine  (Sec.  5.4)  that  computes  clutter. 

4.5.1  Subroutine  TARGS1  (Fig.  4.12) 

This  subroutine  is  called  to  compute  the  apparent  position  (R'.u’.v') 
and  output  signal-to-noise-plus-clutter  ratio  (ZQ)  for  the  case  when  there 
is  single  refraction  (Np  ■ 1)  and  only  one  real  target  is  found  to  be  in 
LISTA.  Although  this  case  is  basically  simple,  monopulse  and  split-gate 
simulations  (MN0PLS1,  SPLTGAT)  must  be  called  into  play  when  appropriate, 
and  the  clutter  must  be  computed  (FBCLTR,  Sec.  5.4). 
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S/V  EVENT 

\ ^ 


SET: 

u,  ' u'(l) 

v,  * VO) 
FROM  L1STA 


CALL  SUB  "MNOPLS1 " TO 
COMPUTE  APPARENT  SINE- 
SPACE  COORDINATES, 

(u1  , V),  & A SINGLE 
TARGET  AT  (u, , v, ) 


SET: 

R*  = R'(1)  FROM  LISTA 
AND 

u‘  - u 


CALL  SUB  "SPLTGAT"  TO 
COMPUTE  INDICATED 
RANGE,  R',  AND  AMPLITUDE 
FACTOR,  FM,  OF  OUTPUT 
SIGNAL 


S * S,  X(l) 


S * S,  (FM)' 


P = 2.0 

*-r  ‘"oiVV'i 


r « o.oi 
\ ^ 


1 

• 1.0 

NO 

LDOP  = 1 

sin  Y , ”P 
Y 1 

s * s/ldop 

'pc  > «' 

\ ^ 


CALL  SUB  "FBCLTR"  TO 
COMPUTE  FIREBALL  AND 
DUST  CLOUD  CLUTTER,  C, 
FOR  THIS  WAVEFORM  AND 
BEAM  POSITION 


XT  * XT  + C 


zi  " S/XT 


Zi  +1 


Zo  1 Z1/Lx 


RETURN 


Figure  4.12.  Subroutine  TARGS1 


Inputs : 


MLST,  LISTA,  (u^),  Rq,  F^,  fd,  XT 

Outputs:  (R',u',v'),  Zq 

4.5.2  Subroutine  TARCM1S  (Fig.  4.13) 

This  subroutine  is  called  to  compute  R',  u',  v',  and  ZQ  for  the 
case  when  there  is  single  refraction  (Np  - 1)  but  more  than  one  target  is 
found  to  be  in  LISTA.  Direct  interference  between  targets  in  the  same 
resolution  cell,  as  well  as  equivalent  noise  due  to  targets  in  the  range 
sidelobes,  where  applicable  (subroutine  SLNINT),  is  taken  into  account. 

Inputs:  MLST,  LISTA,  (u  ,v  ),  S, , F , X„ 

v ’ o’  o 1 pc  T 

Outputs:  (R'.u'.v1),  Zq 

4.5.3  Subroutine  TARGMSV  (Fig.  4.14) 

This  subroutine  is  called  to  compute  R',  u* , v',  and  ZQ  for  the 
case  of  refraction  through  a striated  region  (Np  > 1);  it  is  applicable 
only  to  search  and  verification  events  (no  range  gating),  since  it  makes 
the  same  assumptions  as  were  mentioned  in  Sec.  4.2.1  in  connection  with 
subroutine  POSSV.  A simplified  statistical  model  of  the  "spread  target 
cloud"  is  used  to  compute  the  apparent  target  range  R'  ; since  monopulse 
is  not  used,  the  system  estimate  of  the  target's  angular  position  is  just 
that  of  the  beam  axis.  A special  subroutine  (SLNINTM)  is  used  to  compute 
equivalent  noise  from  spread  targets  in  the  range  sidelobes  (applied  only 
to  those  cases  where  pulse  compression  is  used  in  search  and  verification 
events)  with  a statistical  model. 

Inputs:  (uo,Vo),  fp^,  *5R^,  ^AR*  LISTB,  ITO,  S^,  X^, 

Outputs:  (R'.u'.v'),  Zq 

4.5.4  Subroutine  TARGMTS  (Fig.  4.15) 

This  subroutine  is  called  to  compute  R' , u',  v',  and  ZQ  for  the 
case  of  striated  refraction  (Np  > 1)  when  the  event  type  is  one  that  de- 
fines a range  gate  position  and  width  (track  initiation,  track).  As 
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Figure  4.13.  Subroutine  TARGM1S 
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Figure  4. 14,  Subroutine  TARGMSV 


AN  405M 


CALL  SUB  "MLTSPLT"  TO 
COMPUTE  INDICATED  RANGE, 
R\  AND  AMPLITUDE  FACTOR, 
FM,  OF  OUTPUT  SIGNAL 


CALL  SUB  "FBCLTR"  TO 
COMPUTE  FIREBALL  AND 
DUST  CLOUD  CLUTTER,  C, 
FOR  THIS  WAVEFORM  AND 
BEAM  POSITION 


CALL  SUB  "MNOPLSM"  TO 
COMPUTE  APPARENT  SINE- 
SPACE  COORDINATES,  (u\ 
v').  OF  TARGET 


CALL  SUB  "SLNINTM"  TO 
COMPUTE  EQUIVALENT  NOISE 
INCREMENT,  AXp  DUE  TO 

RANGE  SIOELOBE  TARGET 
INTERFERENCE 


S * Sj (FM)* 


xT  + axt  + C 


Zo  " VLX 


RETURN 


Figure  A. 15.  Subroutine  TARGMTS 
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above,  a simplified  statistical  model  of  the  interaction  between  the  "spread 
target  cloud"  and  the  radar  is  used.  Monopulse  and  split-gate  range  marking 

are  simulated  in  this  way,  and  range -side lobe  interference  is  modeled  by  an 

equivalent  noise  increment. 

Inputs:  S1§  XT 

Outputs:  (R',u',v'),  Zq 

A. 6 SUBROUTINES  CALLED  BY  TARGET-POSITION  SUBROUTINES 

As  Fig.  4.1  and  the  preceding  flow  charts  showed,  the  "TARG"  sub- 
routines call  a number  of  other  subroutines  to  help  compute  the  apparent 
target  position  and  signal-to-noise-plus-clutter  ratio.  These  may  be 
grouped  as  follows: 

• Monopulse  angle  output:  MN0PLS1,  MNOPLS,  MNOPLSM 

• Split-gate  range  output:  SPLTGAT,  MLTSPLT 

• Range  output  for  ungated  signals:  SVPEAKS 

• "Noise"  due  to  targets  in  range  sidelobes:  SLNINT , SLNINTM 

Functions  AMBGN  and  QINV  are  in  turn  called  by  some  of  these  subroutines. 

The  following  sections  describe  these  subroutines  and  functions  in  the 
order  Just  given.  Subroutine  FBCLTR,  which  computes  the  clutter  contribu- 
tion to  Zq  , is  described  with  the  propagation  routines,  in  Sec.  5.4. 

4.6.1  Subroutine  MN0PLS1  (Fig.  4.16) 

This  subroutine  computes  the  apparent  sine-space  coordinates  output 
by  the  monopulse  system  when  only  a single  target  is  present  in  LISTA.  If 
the  target  point  lies  within  the  antenna  pattern  mainlobe,  the  monopulse 
indication  is  identical  with  the  target's  angular  position.  However,  if 
the  target  point  lies  outside  the  mainlobe,  the  monopulse  system  would 
indicate  a false  position  within  the  mainlobe,  and  this  subroutine  computes 
this  apparent  indicated  position. 

Inputs:  (u^v^,  p,  (uq,vo) 

Outputs:  (u',  v') 
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A. 6. 2 Subroutine  MNOPLS  (Fig.  4.17) 

This  subroutine  simulates  the  operation  of  the  radar  monopulse  when 
there  is  more  than  one  target  in  LISTA.  It  computes  the  in-phase  and 
quadrature  components  of  the  sum  and  difference  channel  signals  contributed 
by  each  target,  adds  them,  and  calculates  the  apparent  sine-space  coordinates 
of  the  indicated  "target."  (Appendix  C gives  details.)  Only  one  beam  shape 
is  provided  for:  a circular  beam  with  tapered  sidelobe  level. 

Inputs:  MLST,  LISTA,  R^,  (5Rj,  Sw,  (U0»VQ)»  P» 

Outputs:  (u',v') 

4.6.3  Subroutine  MNOPLSM  (Fig.  4.18) 

This  subroutine  performs  the  same  function  in  the  case  of  refraction 
through  a striated  region  (Np  > 1)  as  is  performed  by  subroutine  MNOPLS  for 
the  single-refraction  case  (Np  ■ 1).  Using  a simplified  statistical  model 
of  the  "spread  target  cloud,"  it  simulates  the  operation  of  the  radar  mono- 
pulse on  all  the  target  clouds  in  LISTB.  Because  of  the  statistical  nature 
of  the  target  model,  the  sum  and  difference  channel  signals  are  calculated 
by  using  a random  phase  approximation  for  the  individual  contribu.. xuns  of 
the  target  elements  (see  Appendix  C).  The  outputs  are  the  apparent  sine- 
space  coordinates  of  the  indicated  "target."  As  with  MNOPLS,  only  one  beam 
shape  is  provided  for:  circular  with  tapered  sldelobes. 

Inpucs : MLST , LISTB , R_ , 6 R., , 6 w , o A n , p,  (u  _ , v _ ) 

' T O J.  O O 

Outputs:  (u',v') 

4.6.4  Subroutine  SPLTGAT  (Fig.  4.19) 

This  subroutine  simulates  the  operation  of  a split-gate  range  tracker 
on  the  set  of  individual  target  images  in  LISTA.  Since  split-gate  range 
tracking  is  used  only  in  those  events  which  also  employ  a pulse-compression 
waveform,  the  subroutine  assumes  that  the  signal  processing  is  matched  to 
the  known  form  of  the  received  signal  return  from  a single  target  unaffected 
by  propagation  errors.  However,  if  there  are  several  target  images  in  LISTA, 
their  contributions  to  the  received  video  signal  output  are  all  computed, 
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Figure  4.17.  Subroutine  MNOPLS 
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COMPUTE  RANDOM  MONOPULSE 
INDICATION: 
p ■ 2.2  a Sw 
4 * 2»b 

u'  * Uq  ♦ p COS  4 

v'  ■ »Q  ♦ c sin  4 


INITIALIZE : 

X.  • X / 100 

L ffl 

6 * P/2 
Ju  • Jv  0 
Su  * Sv  * 0 
Du  • Dv  ■ 0 

b * [1*12(^)Z]1/2 


LOOP  ON  ALL  TARGETS 
IN  LI  STB  (IT)  |7| 


R,  < R'(|T)  « R? 
\ T ^ 


X( IT)  * X. 


0,  • tin (b[u‘ (IT)-Uj,]! 
Qy  ■ tin  f fi[v 1 ( IT)-yQ]i 


|0  | « 100 


lu 

■ 1 

Ju 

• Ju  4 1 

L ■ 

[R£-R'(IT)]/b 

A • 
B • 

*N  (C) 

X(IT)*A2 

Su 

■ Su  4 P 

Du 

* Du  4 B»q‘ 

END  TARGET  LOOP 


SELECT 
» cU  (-1,1) 
d ■ 2.2  i iw 
u'  • «0  ♦ d 


Y,  ■ 2 (Dv*Sv)’/? 
i2  ■ Sv  - Dv 
6 • ATAN2  (Y,,Y2) 
»'  * v„  ♦ e/p 


10,1  1 100 


tu  • 0^ 

Trts 


Jv  * Jv  ♦ I 
Sv  * Sv  ♦ B - 
Dv  ■ Dv  4 B*Q‘ 


t ■ [*£-*’ (JT)]/b 

* • X(l  U) 

B • X(IT)*A2 


Figure  4.18.  Subroutine  MNOPLSM 


108 


Caking  into  account  their  proper  phase  relationships  and  any  dispersive 
effects  which  may  have  distorted  the  received  waveform.  The  resulting 
video  envelope  may  therefore  no  longer  be  matched  to  the  split-gate  re- 
ceiver function,  leading  to  an  error  in  the  indicated  range  of  the  "target." 
Since  this  is  what  would  be  happening  in  a real  radar  under  the  same  con- 
ditions, the  indicated  range  output  is  precisely  what  one  desires  to  simu- 
late, errors  and  all.  The  subroutine  also  computes  the  indicated  signal 
amplitude  which  would  be  output  by  the  split-gate  system,  which  can  be 
immediately  converted  into  an  indicated  "peak"  signal  power  for  use  in 
calculating  the  effective  SNR. 

Inputs:  MLST,  LISTA.  W,  Rc»  6R|,  f , 6R 

Outputs:  R',  FM 

4.6.5  Subroutine  MLTSPLT  (Fig.  4.20) 

This  subroutine  performs  the  same  function  in  the  case  of  refraction 
through  a striated  region  (Np  > 1)  as  is  performed  by  subroutine  SPLTGAT 
for  the  single-refraction  case  (Np  - 1).  Using  a simple  statistical  model 
of  the  "spread  target  cloud,"  it  simulates  the  operation  of  the  split-gate 
range  tracker  on  all  the  target  clouds  in  LI STB.  Because  of  the  statisti- 
cal nature  of  the  target  model,  the  random  phase  approximation  is  used  in 
combining  the  contributions  of  different  targets  to  the  video  signal. 
Dispersive  effects  on  the  return  waveform  are  taken  into  account.  The 
outputs  of  this  subroutine  are  the  indicated  range  and  amplitude  of  the 
"target"  return. 

Inputs:  MLST,  LISTB,  W,  Rc,  6R^,  6R^ 

Outputs:  R' , FM 

4.6.6  Subroutine  SVPEAKS  (Fig.  4.21) 

This  subroutine  computes  the  indicated  target  range  and  peak  ampli- 
tude for  search  and  verification  events,  which  do  not  use  a range  gate. 

In  this  case,  the  video  signal  output  due  to  the  contributions  of  all  the 
targets  in  LISTA , combined  in  their  correct  phase  relationships  and  taking 
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into  account  any  dispersive  distortions  of  the  received  waveform,  is  com- 
puted and  sampled  at  intervals  within  one  range  resolution  distance  on 
either  side  of  of  the  primary  target  of  current  interest.  The  position 
and  amplitude  of  the  highest  peak  thus  discovered  are  the  outputs. 

The  subroutine  also  computes  the  positions  and  amplitudes  of  all 
other  video  signal  peaks  within  the  indicated  region;  these  outputs  are 
not  used  in  the  current  version  of  ROSCOE  but  may  be  required  by  more 
sophisticated  versions  to  follow. 

Inputs:  ML  ST,  LiSTA,  ITO.  6R,\  6R  . F . f , f . 

— - — j.  o pc  o a 

Outputs:  R',  FM 

4.6.7  Subroutine  SLNINT  (Fig.  4.22) 

This  subroutine  computes  the  equivalent  noise  increment,  AXj  , which 
represents  the  effects  of  all  individual  target  images  in  the  range  side- 
lobes.  Such  targets  are  considered  to  contribute  too  little  to  be  worth 
including  them  in  the  calculation  of  the  video  signal;  Instead,  they  are 
each  assumed  to  add  to  the  noise  a term  proportional  to  their  peak  ampli- 
tude factor,  modified  by  the  signal's  current  mean  sidelobe-to-mainlobe 
ratio  (taking  dispersive  effects,  if  any,  into  account).  It  is  only  used, 
of  course,  for  signals  which  have  range  sidelobes. 

Inputs:  R'.  MLST.  LISTA.  L_.  S_.  . S.  . R . W,  6R' , 6R,' 

O U bL  1C  OX 

Outputs:  AX.J. 

4.6.8  Subroutine  SLNINTM  (Fig.  4.23) 

This  subroutine  performs  the  safl*  function  in  the  case  of  refraction 
through  a striated  region  (Np  > 1)  as  is  performed  by  subroutine  SLNINT  for 
the  single-refraction  case  (Np  ■ 1) . Using  a simple  statistical  model  of 
the  "spread  target  cloud,"  it  computes  the  equivalent  noise  Increment, 

AXt  , which  represents  the  effects  of  all  those  clouds  or  parts  of  clouds 
which  lie  in  the  range  sidelobes.  It  is  only  used,  of  course,  for  signals 
which  have  range  sidelobes. 
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Figure  4.22.  Subroutine  SLNINT 
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Figure  4.23.  Subroutine  SLNINTM 
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Inputs: 


MLST,  LISTS,  1^.  S^.  o^R, 


<5R* , 6R‘ , W,  R 
o 1 c 


*;• 

Outputs:  AX^ 

4.6.9  Function  AMBCN(X)  (Fig.  4.24) 

This  function  computes  the  normalized  ambiguity  function,  XfjW  , 
for  a given  waveform,  as  a function  of  range,  X , from  the  center  of  the 
function.  Dispersive  effects  are  included,  and  it  is  normalized  so  that 
Xn(0)  - 1 . Two  different  kinds  of  waveforms  are  defined:  a simple  rec- 

tangular pulse  without  frequency  modulation  (FpC  - 1),  and  a pulse- 
compression  signal  (FpC  >>  1)  whose  characteristics  fall  roughly  between 
those  of  Hamming-  and  Hanning-weighted  waveforms.  It  should  also  be  noted 
that  the  ambiguity  function  is  computed  only  as  a function  of  range:  in 

affect,  it  is  assumed  that  <J>  (the  doppler  frequency  difference)  - 0; 
wherever  doppler  mismatch  effects  are  required  they  are  calculated  separately. 
T',e  reason  for  this  is  that  function  AMBGN  attempts  to  take  into  account 
the  waveform  distortions  caused  by  dispersion  along  the  propagation  path(s). 
Although  this  can  be  done  for  the  range  coordinate,  an  exact  formulation 
for  the  dispersive  effects  on  the  doppler  coordinate  has  yet  to  be  derived  t 
an  approximation,  which  separates  the  two  dimensions  of  xN  > is  currently 
used. 

Inputs:  X,  Fp<;,  5R;,  5RJ,  SgL 

Outputs:  Xjj 

4.6.10  Function  QINV(s)  (Fig.  4.25) 

This  function  is  used  by  the  split-gate  range  tracking  subroutines, 
SPLTGAT  and  MLTSPLT , to  locate  the  apparent  target  point  with  respect  to 
the  center  of  the  range  gate.  This  function  is  specifically  matched  to  the 
weighted  waveform  (FpC  >>  1)  mentioned  above,  and  would  give  the  correct 
range  of  the  target  point  if  there  were  only  a single  target  in  the  resolu- 
tion cell  and  if  there  were  no  dispersive  distortion  of  the  signal.  In  the 
presence  of  either  or  both  of  these  "anomalies,"  an  incorrect  position  is 
computed  for  the  target,  which  is  just  what  a real  split-gate  range  tracker 
would  do  under  the  same  circumstances.  See  Appendix  C for  details. 
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Figure  A. 24.  Function 


FIND  Z(X)  BY  INTERPOLATION  OR  OTHERWISE  IN 
THE  FOLLOWING  TABLE: 


X 

Z 

X 

Z 

0 

0 

1.2714 

1.6000 

0.1000 

1.3138 

1.7000 

0.1993 

0.2000 

1.3514 

1.8000 

0.2975 

0.3000 

1.3843 

1.9000 

0.3940 

0.4000 

1.4127 

2.0000 

0.4884 

0.5000 

1.4368 

2.1000 

0.5801 

0.6000 

1.4569 

2.2000 

0.6686 

0.7000 

1.4733 

2.3000 

0.7536 

0.8000 

1.4865 

2.4000 

0.8347 

0.9000 

1 .4966 

2.5000 

0.9115 

1.0000 

1.5042 

2.6000 

0.9837 

1.1000 

1.5095 

2.7000 

1.0512 

1.2000 

1.5130 

2.8000 

1.1138 

1.3000 

1.5151 

2.9000 

1.1714 

1.4000 

1.5160 

3.0000 

1.2239 

1.5000 

1.5162 

3.1000 

1.5163 

3.1416 

1 

QINV 

= Y*Z 

^ RETURN  ^ 


Figure  4.25. 
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Function  QINV 
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Inputs:  s 

Outputs:  QINV 

4.7  THRESHOLDING  AND  RADAR  MEASUREMENT  ERRORS 

As  explained  earlier,  the  output  of  subroutine  RADMODS  is  the  "apparent" 
target  coordinates  (R' ,u',v'),  which  include  the  effects  of  the  environment 
on  the  radar's  measurements  but  do  not  Include  the  inherent  errors  of  the 
radar  itself.  These  are  added  by  another  subroutine,  MEASERR,  called  not  by 
RADMuDS  but  by  the  event  subroutine  that  called  RADMODS. 

A final  subroutine,  XTHRSHS  (also  called  by  the  event  subroutine), 
compares  the  RADMODS  output  ZQ  with  a threshold  signal-to-noise-plus- 
clutter  ratio. 

4.7.1  Subroutine  MEASERR  (Fig.  4.26) 

This  subroutine  just  adds  radar  measurement  and  bias  errors  t j the 
apparent  target  coordinates  which  are  output  by  subroutine  RADMODS.  The 
calculation  involves  selecting  random  samples  from  appropriate  zero-mean 
distributions  having  variances  given  by  a standard  radar  formula,  and  adding 
fixed  biases  as  required. 

Since  the  angular  coordinates  are  not  measured  in  search  and  verifica- 
tion events,  but  are  merely  set  equal  to  the  commanded  receive  beam  axis 
coordinates  (u0,v0),  subroutine  MEASERR  does  not  compute  any  measurement 
errors  for  these  coordinates  for  these  events.  Then  we  have  just 

u = u'  and  v = v* 
m m 

6v^  , remain  constant  for  each  radar 
throughout  a given  computer  run.  They  may  be  input  directly  by  the  user. 

An  alternative  procedure  is  to  input  the  standard  deviations  of  their  dis- 
tributions for  each  radar  type;  then  at  the  beginning  of  each  run  the  values 
6v.  to  be  used  for  each  radar  during  that  run  may  be  drawn 


The  bias  errors. 


6V 


o 


119 


as  random  samples  from  the  appropriate  zero-mean  Gaussian  distributions. 
In  either  case,  the  values  of  these  bias  errors  for  each  radar  are  stored 
in  appropriate  datasets  for  use  in  subroutine  MEASERR. 


The  random  measurement  errors,  6Rf,  6uf,  6vf  , are  sampled  each 
time  from  zero-mean  Gaussian  distributions  with  standard  deviations  oR, 
ou,  a which  are  computed  within  subroutine  MEASERR.  The  variances  of 
these  distributions  are  found  from  a formula  of  standard  type: 


2 

a 

q 


2 _2 
o P.B 

A2  + -a-SL 

q Zo 


where  "q"  stands  for  R,  u,  or  v , respectively.  Here  Z0  is  the  out- 
put signal-to-noise-plus-clutter  ratio  from  RADMODS.  The  six  constants 
{A^,B^}  are  input  by  the  user  for  each  radar  type;  their  values  are 
usually  known  as  a result  of  routine  error  analyses  of  a given  radar  type. 


The  factor  Pq  in  the  above  equation  is  related  to  the  spread  In 
effective  range  resolution  due  to  dispersion.  It  obviously  does  not  apply 
to  angular  measurements;  hence,  ■ pv  = 1 . For  range  measurements,  it 

may  be  shown  that  pR  ® , where  Lg  is  the  dispersive  loss  factor  cal- 

culated in  subroutine  REF1S  for  the  current  radar  waveform  under  present 
conditions. 


The  outputs  of  this  subroutine  are  (R  ,u  ,v  ).  Note  that  in  Sec.  3, 

m m m 

these  outputs  were  denoted  in  different  ways  for  different  event  types: 


Search: 
Verification: 
Track  Initiation: 
Track: 


(VvV 

(Rj^u^v^)  and  (R2,u2,v2)  as  appropriate 
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Inputs:  Zq,  Lp,  (R'ju'jv'),  {A^,A^,A^,Bj^,B^,B^} , {6R^»6u^|6v^} 

Outputs:  (R  ,u  ,v  ) 

m m m 

4.7.2  Subroutine  XTHRSHS  (Fig.  4.27) 

This  subroutine  determines,  In  a very  simplified  manner,  whether  the 
output  signal-to-noise-plus-clutter  ratio  (ZQ)  computed  by  subroutine 
RADMODS  does  or  does  not  cross  a predefined  detection  threshold.  The  de- 
termination is  made  simply  by  computing  the  output  value  equivalent  to  the 
input  threshold  (S/N)^,^  which  Is  one  of  the  parameters  associated  with  the 
basic  wave  form/ event  dataset,  and  then  comparing  Z0  with  this  value.  The 
output,  XTH  , is  a Hollerith-valued  constant,  either  NO  or  YES. 

Inputs:  IFDET , Zq,  (S/N)Ti 

Outputs : XTH 
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ENTER 


Figure  4.27.  Subroutine  XTHRSHS 
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5 


PROPAGATION 


The  following  propagation  models  have  been  included  in  ROSCOG: 

1.  Absorption 

2.  Fireball  Noise 

3.  Refraction 

4.  Clutter 

5.  Low-Altitude  Fireball  Multipath 

5.1  ABSORPTION 

Absorptions  due  to  fireballs,  betas,  delayed  gammas  and  neutrons, 
and  prompt  X-rays  and  neutrons  are  computed  when  subroutine  ABSORB  is 
called  in  the  program.  This  routine  in  turn  calls  a series  of  subroutines 
which  perform  the  individual  absorption  calculations.  The  calculations 
are  terminated  if  at  any  time  the  integrated  absorption  (DBSUM)  exceeds 
a specified  threshold  (DBMAX) . 

5.1.1  Assumptions 

The  attenuation  calculation  consists  of  a numerical  integration  of 
the  non-deviative  incremental  absorption  equation 


3A 

3r 


7.3  x 10"8N  (v  + v .) 
e em  ei 

f2  + (0.775v  + v ,)2 

em  ei 


dB/cm 


where  Ng  ■ electron  density,  cm 

f ■ radar  frequency,  MHz 

v ,v  . ■ average  electron-neutral  and  electron-ion  collision 
em  el  ^ 

frequencies  divided  by  2tt  * 10°  . 


The  non-deviative  assumption,  although  not  strictly  valid,  is  made  to 
reduce  computation  time.  It  permits  one  to  decouple  the  absorption  and 
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refraction  calculations  so  that  a relatively  fast  absorption  computation 
can  be  made  along  the  path  of  interest. 

The  error  introduced  by  this  assumption  is  expected  to  be  small 
except  under  severe  refractive  conditions;  under  these  conditions,  the 
system  is  probably  seriously  degraded  by  refractive  errors  to  the  extent 
that  a precision  absorption  calculation  is  unwarranted. 

5.1.2  Procedure 

Subroutine  ABSORB  computes  absorption  between  two  points  due  to  all 
significant  radiation  sources.  This  subroutine  calls  FBA3S  to  compute 
fireball  and  dust  absorption,  and  DELABS  to  compute  absorption  due  to  all 
other  sources. 

Fireball  absorption  is  computed  by  simply  integrating  the  incremental 
absorption  as  the  radar  line-of-sight  intersects  each  fireball  region.  A 
more  complex  integration  procedure  is  used  in  DELABS.  The  scheme  follows 
these  steps: 


1.  An  array  of  mandatory  points  on  the  path  of  Interest  is  found 
where  absorption  may  be  significant  (e.g.,  for  fireball  ab- 
sorption the  point  of  closest  approach  to  the  fireball  center 
is  used) . 

2.  The  mandatory  points  are  ordered  in  range  from  the  radar — 
shortest  range  first. 

3.  The  integration  is  started  at  the  first  mandatory  point  and 
proceeds  down  in  range  until  one  of  the  following  occurs: 

(a)  range  ■ 0;  (b)  incremental  absorption  falls  below  some 
minimum  threshold;  (c)  Incremental  absorption  is  less  than 
5%  of  the  total  integrated  absorption. 

4.  The  integration  is  then  continued  u£  in  range  from  the  first 
mandatory  point  until  either  the  maximum  range  is  reached 
(e.g.,  the  range  to  the  target)  or  3(b)  or  3(c)  is  satisfied. 
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5. 


For  subsequent  mandatory  points  that  have  not  been  included 
in  the  above  integration,  the  same  procedure  is  followed, 
except  that  when  integrating  down  the  path  the  minimum  or 
limiting  range  is  the  last  integrated  point  rather  than  0. 

In  the  integration  procedure,  a variable  integration  step  size  is 
used  depending  on  the  magnitude  of  the  current  and  last  computed  values 
of  incremental  absorption.  The  increment  in  slant  range  DELINC  is  given 
by: 


DELINC  - AMINKAMAX1  (DELINC*  * DBINC1/DBINC,DSMIN)  ,DSMAX) 

where  DSMIN  - .5  * DSNOM 

DSMAX  « 2 * DSNOM 

DSNOM  « AMIN1 (AMAX1 (DSNOM ' , .1  x 105) , 5 * 105) 

DSNOM'  ■ An  internally  set  value  (e.g.,  .05  x target  altitude 
for  fireball  absorption) 

DBINC1  - Incremental  absorption  for  last  calculation  (taken 
as  zero  initially) 

DBINC  » Current  calculated  incremental  absorption 
DELINC'  * DSNOM  for  first  integration  step 

5.1.3  Flow  Charts 

Flow  charts  for  the  absorption  subroutines  are  presented  in  Figs.  5.1 
to  5.3. 

There  are  a number  of  flags  and  variables  used  in  the  code  to  in- 
dicate integration  direction  and  integration  limits.  Some  of  these  are: 

* 

! 

i 

i 

i 


SI  ■ slant  range  to  the  _ith  mandatory  point 
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4 


* 

» 


S2  - maximum  slant  range  already  integrated  (“0,  initially) 

k - Integration  direction  (denoted  DIRECTN  in  code)  * +1,  up 
path;  -1,  down  path 

L ■ flag  to  Indicate  current  Integration  step  is  starting  at 
point  of  closest  approach  (denoted  ATPOINT  in  code)  - 1, 
if  a mandatory  point;  0,  otherwise 

kl  - flag  to  indicate  integration  has  reached  lower  integration 
limit  (S2)  (denoted  ATSMIN  in  code)  - 0,  if  at  S2;  1, 
otherwise 


5.1.4  Inputs  and  Outputs 

Subroutine  ABSORB  is  called  by  subroutine  RADIS  of  the  radar 
signal  processing  module.  RADIS  feeds  ABSORB  the  radar  position. 


current  calculation  time,  radar  frequency,  and  maximum  one-way  absorption 
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Figure  5.2.  Subroutine  FBABS 
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Figure  5.3.  Subroutine  DELABS 
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of  interest  (equal  to  the  total  dB  required  to  reduce  the  signal-to-noise 
ratio  to  a value  Just  below  the  detection  threshold).  ABSORB  returns  the 
total  one-way  integrated  absorption  over  the  current  path  of  interest. 

The  inputs  and  outputs  of  the  absorption  routines  are  summarized  as 
follows: 


ABSORB  Inputs: 

DBMAX 

Maximum  one-way  absorption  of 
interest,  dB 

TIME 

Current  time,  seconds 

FREQ 

Radar  frequency,  MHz 

R(3) 

Radar  position,  cm 

ABSORB  Output: 

DBSUM 

Total  integrated  absorption,  dB 

FBABS  and  DELABS  Inputs: 

DBMAX 

Maximum  absorption  of  interest,  dB 

R(3) 

Radar  position,  cm 

UV(3) 

Unit  vector  along  radar  line-of- 
sight 

SMAX 

Maximum  range  of  interest,  cm 

ELEV 

Current  radar  beam  elevation,  deg 

MP 

Pointer  to  mandatory-points  dataset 

TIME 

Time,  seconds 

FREQ 

Radar  frequency,  MHz 

FBABS  and  DELABS  Output: 

DBSUM 

Total  integrated  absorption,  dB 

5.2  FIREBALL  NOISE 

All  fireballs  in  the  main  beam  of  the  radar  antenna  are  assumed  to 
be  potential  noise  sources  in  ROSCOE.  Subroutine  NOISE  computes  the 
effective  noise  temperature  along  the  radar  line-of-sight  from  each 
fireball  satisfying  this  condition. 
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The  noise  temperature  calculation  performed  in  ROSCOE  attempts  to 
account  for  non-uniform  temperature  and  emissivlty  distributions  through 
the  fireball  edge  by  numerically  integrating  through  this  region.  The 
form  of  the  expression  for  the  noise  temperature  contribution  from  a 
particular  fireball  is 


-0* 1Adb 

Tef  " 10  x Teff  * GAIN(M) 


where 


“0* 

Teff  " /T  x e x 10  d£ 


.Jt 


GAIN  (9 


•*>  mfJ- 


4tt 


sin  $ dcp  d0 


and 


Apg  - cumulative  absorption  from  outer  edge  of  fireball  to 

radar  (or  next  fireball  if  there  is  more  than  one),  dB 

T ■ temperature  at  point  of  interest  along  path  l , °K 

e ■ emissivity  at  point  of  interest 

■ 1 - io'"'1*DBI,,c 

ABBS  " cumu^atdve  absorption  from  outer  £dge  of  fireball  to 
point  of  interest  Inside  fireball,  dB 

^DBINC  " incremental  absorption  for  point  of  interest,  dB 


5.2.1  Assumptions 

The  following  assumptions  are  made  in  the  NOISE  model: 

1.  The  radar  beamwidth  is  narrow  and  the  main  contributors  to 
the  noise  temperature  are  fireballs  in  the  main  beam. 


2. 


Integration  over  the  main  beam  pattern  is  used  to  find  the 
effective  contribution  to  the  noise  for  each  fireball,  but 
the  noise  is  propagated  along  a path  coincident  with  the 
antenna  pointing  direction.  This  simplifies  the  noise 
absorption  calculation  sincp  only  a single  path  has  to  be 
considered  for  all  fireballs. 

3.  A numerical  integration  is  performed  through  the  fireball 
edge  to  account  for  non-uniform  temperature  and  emissivity 
distributions  as  discussed  above. 

5.2.2  Procedure 

The  noise  calculation  follows  these  steps  (see  Fig.  5.4): 

1.  An  array  of  mandatory  points  is  found  where  noise  temperature 
may  be  significant. 
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2.  The  points  are  range  ordered — shortest  range  fiom  the  radar 
f irst . 

3.  The  calculation  is  started  at  the  first  point  and  proceeds 
down  in  range  until  the  temperature-emissivity  product  drops 
below  rO  of  the  svstetn  temperature.  At  each  step  along  the 
path,  the  temperature-emissivity  product  and  incremental 
absorption  are  stored — step  (c). 

4.  After  reaching  this  point,  the  calculation  reverses  direction 
and  begins  integrating  out  in  range  over  the  fireball  until 
the  starting  point  is  reached  (halfway  over  the  fireball 
intersection  path)— step  (d). 

5.  The  absorption  is  then  computed  from  the  fireball  edge  (SPT) 
to  the  radar  for  to  the  next  fireball  down  in  range  (S2)) — 
step  (e). 

6.  The  procedure  then  returns  to  Mandatory  Point  No.  1,  and 
continues  integrating  out  over  the  fireball  until  the  upper 
edge  is  reached  (i.e.,  temperature-emissivity  again  falls 
below  10"  of  the  system  temperature)  . This  |a  int  is  S2 
for  the  next  sequence  of  calculations — step  .1). 

7.  The  procedure  is  repeated,  beginning  at  the  next  mandatory 
point.  This  is  continued  until  either  all  mandatory  points 
are  covered  or  the  cumulative  absorption  along  the  path  is 
larger  than  30  dB. 

5.2.3  Flow  Chart 

A flow  chart  of  subroutine  NOISE  is  given  in  Fig.  5.5.  A definition 
of  some  of  the  flags  and  variables  used  in  the  routine  follows: 


slant  range  to  the  point  of  closest  approach  to  the  _ith 
fireball  (see  Fig.  5.4). 
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I igure  5.5,  Subroutine  NOISE 
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Figure  5.5  (Cont.) 
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S2^  ■ range  to  the  upper  edge  of  the  (i-l)th  fireball  already 
integrated  over  (see  Fig.  5.4). 

K ■ flag  indicating  the  integration  direction  (denoted  DIRECTN 
in  code)  - +1,  if  integrating  up  the  path;  -1,  if  integrat- 
ing down. 

K1  ■ flag  indicating  that  the  lower  integration  limit  (S2)  has 
been  reached  (denoted  ATSMIN  in  code)  ■ 0,  if  at  S2;  1, 
otherwise. 

L ■ flag  indicating  that  the  integration  step  is  currently  at 
a mandatory  point  (denoted  ATPOINT  in  code)  ■ 1,  if  at 
mandatory  point;  0,  otherwise. 

J ■ number  of  integration  steps  when  integrating  over  fireball. 

M * counter  to  loop  over  J integration  steps. 

5.2.4  Inputs  and  Outputs 

The  NOISE  routine  is  called  by  subroutine  RADIS  of  the  radar  signal 
processing  module.  The  inputs  and  outputs  of  the  NOISE  routine  are 
defined  as  follows: 


Inputs : 

TSYS 

System  temperature,  °K 

• 

TIME 

Current  time 

R(3) 

Radar  position,  cm 

FREQ 

Radar  frequency,  MHz 

Output: 

TEF 

Effective  noise  temperature 

j 
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5.3  REFRACTION 

Roscoe  makes  the  basic  assumption  that  the  number  of  possible 
propagation  paths  to  the  target  In  the  mainbeam  is  either  one  (correspond- 
ing to  the  smoothly  varying  electron  density  gradient  region) , or  much 
greater  than  one  (when  strlations  are  present).  These  two  cases  are 
discussed  separately. 

5.3.1  Gross  Refraction 

In  this  case,  a ray  trace  is  performed,  starting  at  an  altitude  of 
60  km  along  the  radar  beam-pointing  direction  and  ending  at  a point  that 
may  be  displaced  from  the  target  but  is  at  the  same  range  from  the  radar. 
The  incremental  range  and  angular  bias  errors  at  each  step  are  given  by: 


3R  _ 
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where 


u ■ real  part  of  the  Index  of  refraction 

g 1 1/2 

3.18  x io*Ne 

(2.  « lof[f2  + (vra  + »el)2 

_3 

Ng  ■ electron  density,  cm 

VNg^  ■ electron  density  gradient  in  the  _ith  direction, 

-3, 

cm  /cm 

f » radar  frequency,  MHz 

vem,vei  " avera8e  electron-neutral  and  electron-ion  collision 
frequencies  divided  by  2tt  x 10® 


It  is  assumed  that  the  total  angular  error  obtained  from  this  procedure 

is  not  large,  so  that  if  the  beam  is  rotated  just  this  amount,  the 

* 

refracted  beam  will  hit  the  target. 


5.3.2  Refraction  due  to  Strlations 

For  refractive  effects  due  to  ionized  strlations  (the  latter  case 
mentioned  above),  the  assumption  is  made  that  the  electron  density  does 
not  lead  to  excessive  departures  of  a ray  in  the  beam  from  its  nominal 
smooth  trajectory.  With  this  assumption,  it  is  possible  to  develop  a 
theory  which  describes,  in  statistical  terms,  the  effects  of  the 
electron  density  fluctuations. 


A refinement  of  this  procedure  would  be  to  iterate  on  the  beam  pointing 
direction  until  the  refracted  beam  does  hit  the  target. 
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The  theory  presented  here,  termed  the  "thin  screen  approximation," 
differs  from  previous  models^’ ^ in  that  it  systematically  develops  the 
"dogleg"  range  error  and  phase  terms.  Because  of  these  differences,  a 
rather  complete  description  of  the  model  is  given  in  Appendix  F. 


5.3.3  Procedure 

Both  computations  described  above  are  calculated  within  subroutine 
REFRCT  in  ROSCOE.  The  procedure  follows  these  steps  (see  Fig.  5.6): 

1.  The  radar  beam  is  assumed  to  be  straight  to  an  altitude  of 
60  km. 

2.  At  this  point,  an  incremental  range  step  size  (DSINC)  is 
selected  and  the  deflection  of  the  ray  over  this  range 
segment  is  computed. 

3.  In  addition,  random  erro.'s  due  to  strlatlons  over  this 
range  segment  are  computed  and  projected  to  the  target 


range  (SMAX) . 

4.  The  procedure  continues  in  this  way  until  the  range  from 
the  radar  to  the  point  on  the  refracted  ray  (PT^)  is  equal 
to  the  target  range. 


5.  When  this  condition  is  satisfied,  the  program  outputs  bias 
and  random  errors  in  range,  azimuth,  and  elevation  at  the 
final  point,  the  integrated  electron  density  over  the  path, 
and  the  angle  between  the  line-of-sight  and  the  magnetic 
field  at  the  point  of  maximum  electron  density  along  the 
path. 

5.3.4  Flow  Chart 

A flow  chart  of  subroutine  REFRCT  is  presented  in  Fig.  5.7.  This 
routine  requires  Information  from  several  phenomenology  routines  and 
calls  other  propagation  models  as  shown.  A description  of  some  of  these 
subroutines  follows. 
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TARGET 


SMAX  * SLANT  RANGE  TO  TARGET 

RNGBIAS  = TOTAL  RANGE  BIAS  ERROR  DUE  TO  SMOOTHLY  VARYING  REFRACTION 
PT(J)  * COORDINATE  OF  END  POINT  OF  REFRACTED  RAY  PATH  AT  RANGE  OF  TARGET 
RNGFUZ,  ELFUZ,  AZFUZ  « RMS  UNCERTAINTY  IN  END  POINT 
DSINC1  * INCREMENTAL  RANGE  STEP 

* AMAX1  (0.5  * DSNOM,  AMIN1  (2.  * DSNOM,  0.05/DEFL  * DSNOM)) 

DSNOM  * AMIN1  (50.,  RANGE  CORRESPONDING  TO  10  km  ALTITUDE  STEP) 

OEFL  * AMAX1  (INCREMENTAL  DEFLECTION  IN  AZIMUTH  DIRECTION,  INCREMENTAL 
DEFLECTION  IN  ELEVATION  DIRECTION) 


Figure  5.6.  Ray  Trace  Procedure 
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POINT  TO  EARTH 
CENTERED 
COORDINATE 


CfMPuTE  PROMPT 
ENERGY  DEPOSITION 
(ED)  AT  POINT 


CALL  FIELD  (PT.  f) 

CALL  UNI  TV  (F,  FUVUND) 


FIND  UNIT  VECTOR 
(V?)  PERPENDICULAR 
TO  CURRENT  LlNE-OF - 
SIGHT 


GET  ELECTRON  DENSITY 
GRADIENTS  IN  ELEVATION 
AND  AZIMUTH 


USE  GRADIENTS  TO 
COMPUTE  SMOOTH 
BIAS  REFRACTION 


CALL  bias 


COMPUTE  INCREMENTAL 
STEP  SIZE  (DINC)  FOR 
STRIATION  CALCULATION 


COMPUTE  FUZZ  BALL 
AROUND  BIASED  POINT 
AT  TARGET  RANGE 
(AZFUZ.  ELFUZ. 
RNGFUZ) 


FINO  RANGE  TO 
NEXT  frOINT(S) 
AND  RANGE  FROM 
THIS  POINT  TO 
TARGET  (SPT) 


COMPUTE  RMS  POSITION 
UNCERTAINTIES  DUE  TO 
REFRAjCTION 

"call  strTtur 

CALL  FUZINC 


COMPUTE  RANDOM 
ERRORS  IN  RANGE, 
ELEVATION  AND 
AZIMUTH 


FINO  NEXT  INTEGRATION 
POINT  IN  RAY  TRACE 
USING  BIAS  ERRORS 


Figure  5.7.  Subroutine  REFRCT 
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Phenomenology  Subroutines 


ELDENS  Provides  electron  density  and  temperature  at  a 

point  due  to  all  sources 

GRADNE  Computes  electron  density  gradient  at  point  P 
in  a specified  direction  V 

STRCTUR  Computes  structure  of  a striated  region 

BFIELD  Returns  direction  of  B-field  at  a point 

Other  Propagation-Related  Routines 

COLLF  Calculates  electron-neutral  and  electron-ion 

collision  frequencies 

BIAS  Computes  smooth  bias  refractive  errors 

FUZINC  Computes  RMS  position  uncertainties  due  to 

refraction  in  striated  medium  over  the  path 
length 


5.3.5  Inputs  and  Outputs 

Subroutine  REFRCT  is  called  by  subroutine  REF1S  of  the  radar 
signal  processing  module.  The  Inputs  and  outputs  of  REFRCT  are: 


Inputs : 


TIME  Engagement  time,  s 

SMAX  Slant  range  to  target,  cm 

R(3)  Radar  position,  cm 

FREQ  Radar  frequency,  MHz 


Outputs : NP 

XNE 

REFRAC (1-3) 
REFRAC (4-6) 
REFRAC (7) 


Number  of  possible  paths  to  target 
Integrated  electron  density 
Bias  errors  in  R,A,E  coordinates 
Random  errors  in  R,A,E  coordinates 
Angle  between  LOS  and  B-field  at 
point  of  maximum  electron  density 
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SIGBAR 


Amplitude  gain  factor  for  each 


return 

XMU  Angle  between  major  axis  of  uncertainty 

ellipse  and  the  u-axis  in  sine  space 

The  refraction  errors  are  used  to  compute  the  apparent  target  position 
that  the  radar  "sees."  The  Integrated  electron  density  is  used  in  the  esti- 
mation of  the  pulse  dispersion,  and  the  angle  between  the  LOS  and  B-field 
at  the  point  of  maximum  electron  density  is  used  to  determine  the  Faraday 
rotation  angle,  if  appropriate. 

5.4  CLUTTER 

Radar  backscatter  or  clutter  can  occur  from  a variety  of  regions 
following  a nuclear  burst;  in  general,  whenever  the  beam  traverses  a 
region  in  which  the  dielectric  constant  varies.  Clutter  sources  such  as 
dust,  fireballs,  and  field-aligned  regions  are  treated  in  ROSCOE. 

5.4.1  Assumptions 

The  following  assumptions  have  been  made  for  the  clutter  model: 

kk 

1.  The  absorption  of  the  clutter  return  from  the  region  surface 
is  computed  along  the  path  between  the  region's  center  and 
the  radar.  This  allows  a single  path  absorption  Integration 
to  be  used  for  each  clutter  source. 

2.  The  solid  angle  subtended  by  the  clutter  source  is  divided 
into  25  projected  area  segments  for  incremental  clutter 
calculations.  This  makes  the  clutter  Integration  independent 
of  the  geometry  of  the  region. 

k 

Faraday  rotation  is  considered  when  a linearly  polarized  radar  antenna 
is  specified  in  the  input. 

hit 

This  absorption  calculation  includes  the  absorption  through  the  gammas 
surrounding  the  fireball,  if  appropriate. 


3.  The  total  clutter  return  Integrated  over  the  beam  pattern 
Is  multiplied  by  a factor  which  represents  the  fraction  of 
the  clutter  source  in  the  compressed  and  uncompressed 
range  cells;  that  is,  a detailed  integration  is  not  per- 
formed over  the  range  cell. 

5.4.2  Procedure 

The  total  clutter  power  reaching  the  radar  from  all  clutter  sources 
is  computed  in  subroutine  FBCLTR.  This  routine  determines  if  a clutter 
source  is  in  the  compressed  or  uncompressed  range  cell  and,  if  it  is, 
computes  the  fraction  of  the  cell  covered.  The  clutter  power  is  then 
integrated  over  the  antenna  beam  pattern  in  subroutine  CLINT  and 
returned  to  FBCLTR.  Finally,  the  one-way  absorption  is  computed  from 
the  .source  to  the  radar  and  the  clutter  power  is  degraded  by  this  factor 
times  the  fraction  of  the  range  cell  subtended  by  the  source. 

The  clutter  integration  in  CLINT  follows  these  steps: 

1.  Define  the  integration  step  size  (incremental  area)  to  be 
small  compared  to  the  extent  of  the  source.  (A  maximum  of 
25  steps  is  allowed.) 

2.  Direct  a vector  from  the  radar  to  the  center  of  each  of 
these  area  segments. 

3.  Compute  the  clutter  power  per  unit  area  along  the  defined 
vector  intersecting  the  source.  Degrade  the  clutter  power 
by  the  antenna  angular  gain. 

4.  Sweep  in  azimuth  and  elevation.  Integrating  over  the  portions 
of  the  source  which  contribute  significantly  to  the  total 
clutter. 

An  Individual  calculation  is  skipped  when  (1)  the  vector  does  not 
intersect  the  geometric  region  defining  the  source,  or  (2)  the  individual 
segment  contribution  is  a small  fraction  of  the  cumulative  clutter  power 
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(this  is  necessary  for  fireballs  with  diffuse  edges).  The  Integration 
Is  terminated  when  all  25  steps  have  been  considered. 


The  Incremental  clutter  Is  computed  In  subroutine  CLUTINC,  which 
Is  called  by  CLINT.  It  In  turn  calls  routines  to  get  the  surface 
reflection  coefficient  and  applies  a doppler  loss  to  the  incremental 
clutter  power  which  is  computed  by  comparing  the  doppler  shift  of  the 
clutter  source  with  the  doppler  shift  of  the  target.  The  doppler  loss 
Is  computed  from  the  following  equation: 


FDOP  - ABS[ (sin  X) /X] 


3.4 


where 


X - f-  ABS(FDt  - FDC) 
n 

B ■ noise  bandwidth.  Hz 
n 

FDt  ■ expected  doppler  shift  for  the  target  based  on  past 
tracking  measurements,  Hz 

FD_  ■ doppler  shift  of  the  clutter  source 


2f  R 
o c 


, Hz 


f - 

o 

c ■ 
R - 


radar  frequency,  Hz 
speed  of  light,  cm/s 

velocity  of  the  clutter  source  toward  the  radar,  cm/s 


Thus,  the  equation  for  the  incremental  clutter  power  per  unit  area 
(AC/M)  is 


I. 


77  - S * FDOP  x ALBEDO 
AA 
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where  S « signal  for  this  conputatlon 

ALBEDO  - surface  reflection  coefficient 

i 

The  total  clutter  power  for  this  source  is  then  given  by: 

! 

I 

where  GAIN  - antenna  gain  In  the  direction  corresponding  to  the 

_ith  area  segment 

AX  - fraction  of  range  cell  subtended  by  the  source 
DBSUM  - total  abiorption  from  the  source  to  the  radar 

5.4.3  Flow  Charts 

Figures  5.8  and  5.9  are  flow  charts  for  subroutines  FBCLTR  and 
CLINT.  Some  of  the  other  subroutines  and  functions  called  by  these 
routines  are  listed  below. 

ft 

Computes  the  absorption  from  the  fireball 
surface  to  the  radar  (see  Sec.  5.1) 

Computes  the  Incremental  clutter  power  per 
unit  area  for  the  _ith  source 

Computes  the  range  to  the  surface  of  the 
clutter  source  and  the  extent  of  the  source 

Computes  the  angular  gain  for  an  input  vec- 
tor direction  off-beam  center  (see  Sec.  4.4.2) 

Computes  the  solid  angle  subtended  by  the 
clutter  source 


ABSORB 

CLUTINC 

EXTENT 

GOA 

SLDANGL 
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Figure  5.8.  Subroutine  FBCLTR 
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CALL  GOA 
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Figure  5.9.  Subroutine  CLINT 
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5. A. 4 Inputs  and  Outputs 

The  inputs  and  outputs  to  the  three  main  routines  which  make  up 
the  clutter  module  (FBCLTR,  CLINT,  CLUTINC)  are  listed  in  Table  5.1. 

5.5  LOW- ALTITUDE  MULTIPATH 

The  low-altitude  multipath  model  consists  of  four  subroutines. 

The  main  routine  (MLTPATH)  is  called  by  subroutine  REFLST1  (see  Sec.  4.4.1). 
It  determines  if  there  are  any  additional  target  Images  caused  by  low- 
altitude  fireball  multipath  which  might  Interfere  with  the  primary  target 
signal.  The  other  routines  perform  the  following  supporting  computations: 

BOUNCE — Computes  the  apparent  target  image  position  and  signal 
amplitude  for  the  three  additional  returns  that  can  occur  for 
each  signal 

AMPREF — Computes  the  amplitude  reflection  coefficient  at  the 
designated  specular  reflection  point  on  the  fireball  surface 

SPCULAR — Computes  the  specular  reflection  point  on  the  fireball 
surface 

5.5.1  Assumptions 

The  following  assumptions  have  been  made  in  the  MLTPATH  model: 

1.  Only  fireballs  below  60  km  altitude  are  considered  as 
possible  reflecting  surfaces. 

2.  Multiple  reflections  are  ignored;  that  is,  only  a single 
reflection  between  any  object/fireball  pair  is  considered 
(but  all  pairs  are  considered). 

3.  Differences  in  attenuation  along  the  reflected  path  as 
compared  to  the  direct  path  are  ignored  in  computing  the 
received  signal. 

In  addition,  some  simplifying  assumptions  are  made  in  AMPREF 
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and  SPCULAR: 


TABLE  5.1 

INPUT /OUTPUT  DATA  FOR  CLUTTER  ROUTINES 


FBCLTR  Inputs 

(from  calling  sequence): 

RTARG 

XT 

T 

FREQ 

R 

BVEC 

FBCLTR  Inputs 
(from  datasets) : 

DELRU 

DELRC 

FDBAR 

SLRED 

SNMIN 

FBLIS 

FBCLTR  Output: 

SCTRC 

CLINT  Inputs: 

DSPFB 

FPC 

BN 

R 

BVEC 

CLINT  Output: 

CLTR 

CLUTINC  Inputs: 

DSPFB 

V 

R 

FPC 

BN 

FREQ 

CLUTINC  Output: 

CLTINC 

Range  to  target,  cm 
Total  noise  power,  erg/s 
Time  of  calculation,  s 
Radar  frequency,  MHz 

Earth-centered  coordinates  of  radar,  cm 
Boresight  vector  of  radar 

Uncompressed  pulse  range  resolution,  cm 
Compressed  pulse  range  resolution,  cm 
Expected  doppler  shift,  MHz 
Range  sidelobe  reduction  factor 
Signal-to-noise  threshold 
List  of  clutter  sources 

Total  clutter  power  reaching  receiver, 
erg/s 

Pointer  to  clutter  source  of  interest 
Pulse  compression  ratio 
Noise  bandwidth 
Radar  coordinates,  cm 
Boresight  vector 

Clutter  power  output  (at  surface  of 
region),  erg/s 

Pointer  to  clutter  source  of  Interest 
Vector  from  radar  to  source  point  of 
present  calculation 

Earth-centered  coordinates  of  radar,  cm 
Pulse  compression  ratio 
Noise  bandwidth 
Radar  frequency,  MHz 

Incremental  clutter  from  ^th  source, 
erg/s 
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1.  The  amplitude  reflection  coefficient  is  computed,  using  the 
REFC0  routine,  by  inputting  the  proper  angle  of  incidence 
of  the  radar  beam. 

2.  The  fireball  is  spherical  and  the  reflection  point  lies  on 
the  surface  at  the  point  of  closest  approach  to  the  direct 
line-of-sight  (see  Fig.  5.10). 

5.5.2  Discussion 

Flow  charts  are  presented  for  subroutines  MLTPATH  and  BOUNCE  in 
Figs.  5.11  and  5.12.  The  MLTPATH  routine  accesses  the  radar  signal 
processing  dataset  to  get  the  list  of  primary  targets  that  have  been 
acquired.  It  also  accesses  the  propagation  Inputs  dataset  to  get  the 
mandatory  points  list,  which  are  the  points  of  closest  approach  of  the 
fireballs  to  the  current  line-of-sight.  The  routine  then  loops  on  each 
object/point  pair  for  those  mandatory  points  where  the  associated  fire- 
ball is  below  60  km  and  calls  subroutine  BOUNCE. 


ACCESS  RADAR 
SIGNAL  PROCESSING 
DATASET  (NRSP) 

1 

LOOP  ON  ALL  TARGETS 
IN  "LISTA"  (LIST  RS 
(NRSP))  IN  REVERSE 
ORDER 

♦ 

ACCESS  PROPAGATION 
INPUTS  DATASET 
(NPI ) 

J 

LOOP  ON  ALL 
MANDATORY  POINTS 
(MPLIS  MP  (NPI)) 

♦ 

COMPUTE  FIREBALL 
ALTITUDE  FOR  THIS 
MANDATORY  POINT. 
ALTFB  - XMAX  (FBPOS) 
- RBODY 

<^ALTFB  > 60 ^*>1^ 

Tno 

CALL  BOUNCE  (GET 
RUV  COORDINATE  OF 
APPARENT  IMAGES 
AND  AMPLITUDE  OF 
RETURN  "A") 

1 

CREATE  A TARGET 
IMAGE  DATASET 
FOR  1-1,3 
PATHS 

I 

CALL  GOA  (GET 
SIGNAL  D I GRAD 
DUE  TO  TARGET 
0FF-AXIS-“F2" ) 

COMPUTE  TOTAL 
SIGNAL  DEGRAD 
FACTOR  "x" 

* F?/*? 


UNLOCK  RS.  PI 
DATASETS 


— RETURN  J 


^ MANDATORY/ 

\ POINT  J 

v^NO 

PUT  THIS  (1th) 
TARGET  IMAGE  ON 

"LISTA"  IN 
INCREASING  ORDER 
OF  RANGE 

COMPUTE  OFF-AXIS  FACTOR! 
Fc  FOR  EACH  PATH: 

T~T\*¥2*~ 

F2  = F1*F2*A 
F3  = (F2*A)**2 


WHERE:  FI 


DIRECT  PATH 
FACTOR 


Figure  5.11.  Subroutine  MLTPATH 
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Figure  5.12.  Subroutine  BOUNCE 
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BOUNCE  constructs  the  geometry  shown  in  Fig.  5.13  for  the  object/ 
fireball  pair  of  interest.  It  first  calls  subroutine  SPCULAR  to  find 
the  specular  reflection  point  (P0INT1  in  the  figure)  and  then  computes 
the  angle,  0^  , for  the  dogleg  path  and  the  incremental  Increase  in 
range,  DELRANG,  along  this  path.  There  are  three  possible  returns  which 
the  radar  may  receive  in  addition  to  the  direct  or  primary  return.  They 
are: 

1.  Transmit  along  the  dogleg  path,  receive  along  the  direct 
path 

2.  Transmit  along  the  direct  path,  receive  along  the  dogleg 
path 

3.  Transmit  along  the  dogleg,  receive  along  the  dogleg 

I 

In  the  first  case,  the  radar  will  "see"  a target  image  at  a range  equal 
to  the  actual  range  plus  DELRANG  in  the  same  direction  as  the  real  target. 

In  the  second  and  third  cases,  it  will  "see"  the  image  in  the  dogleg 


direction  at  ranges  equal  to  the  range  plus  DELRANG  and  the  range  plus 
twice  DELRANG,  respectively.  This  Information  Is  returned  to  MLTPATH 
along  with  the  amplitude  of  each  signal  relative  to  the  direct-path 
signal. 

In  MLTPATH,  this  information  is  used  to  calculate  the  reduction 
in  signal  power.  The  signal  power  for  these  images  is  reduced  because 
of  reflection  (signal  amplitude)  and  because  the  target  image  is  off 
the  main  beam  axis  in  either  transmit  or  receive  for  the  first  two  cases 
and  in  both  directions  for  the  third  case.  The  resulting  signal  degra- 
dation factor  (see  Fig.  5.11)  for  each  case  is  compared  to  a similar 

expression  for  the  actual  target  (direct  path) . If  the  ratio  of  these 
numbers  is  greater  than  -10  dB,  the  target  image  dataset  associated  with 
this  reflected  path  is  added  to  the  target  list  and  may  ultimately  affect 
the  radar  measurement  for  this  look. 

5.5.3  Inputs  and  Outputs 

The  MLTPATH  subroutine  is  called  in  REFLST1  with  the  radar  position 
and  boresight  vector  in  the  calling  sequence.  In  addition,  the  radar 
signal  processing  and  propagation  inputs  datasets  are  accessed,  as 
mentioned  above.  The  output  of  MLTPATH  is  the  creation  of  new  target 
linage  datasets  which  are  placed  on  the  target  image  list  LISTA. 

Inputs  and  outputs  for  the  other  subroutines  are  carried  through 
the  calling  sequences.  Definitions  for  these  variables  follow: 


BOUNCE  Inputs: 

R(3) 

Radar  position,  earth  centered 

UVB(3) 

Boresight  vector 

T(3) 

Assumed  target  position,  earth  centered 

DSPMP 

Pointer  to  mandatory  point  dataset  of 
interest 

BOUNCE  Outputs: 

A(3) 

Relative  amplitude  of  return  signal  for 
^ith  path 
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RUVX(3,3) 

RUV  coordinates  of  apparent  target 
image  for  ^th  path 

AMPREF  Inputs: 

DSPMP 

Pointer  to  mandatory  point  dataset  of 
Interest 

R(3) 

Radar  position,  earth  centered 

POINTl(3) 

Earth-centered  coordinates  of  specular 
reflection  point  of  fireball  surface 

AMPREF  Output: 

AMP 

Amplitude  reflection  coefficient  (rela- 
tive to  1.0  for  total  reflection) 

SPCULAR  Inputs: 

R(3) 

Radar  position,  earth  centered 

T(3) 

Assumed  target  position,  earth  centered 

DSPMP 

Pointer  to  mandatory  point  dataset  of 
Interest 

SPCULAR  Output: 

POINT (3) 

Earth-centered  coordinates  of  specular 
reflection  point  of  fireball  surface 
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APPENDIX  A 

COORDINATE  SYSTEMS  IN  ROSCOE 

Three  basic  coordinate  systems  are  used  in  the  radar  modeling  for 
ROSCOE,  but  for  two  of  these  the  components  are  sometimes  quoted  in  rec- 
tangular form  and  sometimes  in  spherical  form,  so  in  essence  there  are 
five  sets  of  coordinates: 

1.  Radar  Face  Coordinates 


Rectangular:  p1,p2,p3 

Spherical  (sine-space):  R,u,v 


2.  Radar  Site  Coordinates 


Rectangular  (topocentric) : ri*r2»r,» 

Spherical  (alt-azimuth):  R,A,E 

3.  Geocentric  Rectangular  Cartesian  Coordinates:  R1,R2,R.J 

three  components,  in  each  coordinate  system,  represent  the  llne-of-slght 
vector,  R , from  the  radar  to  the  (apparent)  target. 

The  following  vector  notation  is  used  in  this  appendix.  If  V 
represents  any  vector,  then 

V represents  the  magnitude  of  the  vector  (|^|) 

V represents  the  equivalent  unit  vector  (tf/V) 

represents  the  1th  component  of  the  vector  in  a given  coor- 
dinate system  (i  - 1,2,3) 


j 

1 


j 


i 

1 


I 

i 

i 

1 
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plane  of  Che  radar  face,  and  the  third  Is  normal  to  the  face  in  the  bore- 
sight  (i.e.,  outward)  direction.  As  long  as  the  radar  face  Is  not  locally 

AAA 

horizontal,  the  basis  vectors,  C * h , C » are  defined  as  follows: 
t points  in  the  boreslght  direction 

£ is  locally  horizontal  and  in  the  direction  determined 

A A A 

by  n x c , where  n is  the  local  unit  earth  normal 
at  the  radar 

n is  defined  by  £ * £ 

If  the  face  is  horizontal,  this  definition  becomes  meaningless.  In 
this  case  we  make  a convenient  arbitrary  definition: 

£ is  locally  horizontal  and  points  north 

a 

n is  locally  horizontal  and  points  west 

A A 

; 2 n 

In  this  system,  the  components  of  the  radar  pointing  vector  R are 

- R • £ 

p2  * R • n 

p3  « R • C 

A. 1.2  Spherical 

This  subsystem  is  often  called  "sine-space."  The  components  of  the 
radar  pointing  vector  are 

R - |t| 

A A 

u - p^R  = R • £ 

v = p2/r  = r • n 
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Another  quantity  of  interest  in  this  system  is  0 , the  off-boresight  angle, 
which  is  given  by 


-i  -I  n — 5 

9 ■ cos  (pj/R)  = sin  \u  + v 

A. 2 RADAR  SITE  COORDINATES 


A. 2.1  Rectangular 

This  is  a right-handed  orthogonal  rectangular  Cartesian  coordinate 
system  associated  with  a radar  site.  Because  a site  on  the  surface  of  the 
earth  is  called  a "topos,'  this  is  sometimes  called  the  topocentric  system 
for  short  (of  course  this  does  not  imply  that  the  system  is  restricted  to 
radars  located  on  the  earth's  surface).  Two  of  its  basis  vectors  lie  in 
the  local  horizontal  plane,  and  the  third  points  toward  the  zenith.  These 

A 4 A 

basis  vectors,  x,y,z  , are  defined  as  follows: 

x is  locally  horizontal  and  points  east 

A 

y is  locally  horizontal  and  points  north 


O A A 

z is  locally  vertical  and  points  up  (z  = n) 


In  this  system,  the  components  of  the  radar  pointing  vector  R are 

rj^  - ^ • x 

r2  - t • y 

r3  - 2 . ; 

A. 2. 2 Spherical 

This  subsystem  is  often  called  angular-space  or  the  alt-azimuth 
system. 

The  components  of  the  radar  pointing  vector  ^ are 


R - |^|  = (r^  + r*  + r*) 

A - tan”1(r1/r2) 

E - sin‘1(r3/R) 

A. 3 GEOCENTRIC  RECTANGULAR  CARTESIAN  COORDINATES 

Tnis  Is  a right-handed  orthogonal  rectangular  coordinate  system  with 
origin  at  the  center  of  a spherical  earth;  the  coordinate  system  rotates 
with  the  earth.  It  has  sometimes  been  referred  to  simply  as  "the  Cartesian 
system,"  but  that  term  is  clearly  ambiguous;  here  it  is  called  the  G.R.C. 
system.  Its  basis  vectors,  X,Y,Z  , are  defined  as  follows: 

A 

X lies  in  the  equatorial  plane  and  points  toward  the 
Greenwich  meridian 

A 

Z lies  along  the  earth* s rotational  axis  and  points 
toward  the  geographic  North  Pole 

A A A 

Y - Z x X 

In  this  system,  the  components  of  the  radar  pointing  vector  I are 

^ - St  • X 

R2  - £ • Y 

r3  - it  • Z 

However,  these  are  not  the  components  that  are  usually  desired  or  meant 

when  one  refers  to  "the  GRC  coordinates"  of  a point.  What  are  really  de- 
sired here  are  the  components  of  the  radius  vector  from  the  center  of  the 
earth  to  the  (apparent)  target  point.  Let  V?  represent  this  vector.  Then 

2 - 2 + 2 

where  f is  the  "topos  vector":  the  radius  vector  from  the  center  of  the 
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earth  to  the  radar  site  location.  Thus,  finally,  the  desired  represents- 
tion  of  the  target  point  in  the  G.R.C.  system  consists  of  the  components 


V1  - W • X i Ti  + Ri 


W2  - W • Y = T2  + R2 
W3  - W • Z = T3  + R3 


where  T1*T2»T3  are  the  G.R.C.  components  of  T 
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APPENDIX  B 


DISPERSION  AND  THE  RADAR  AMBIGUITY  FUNCTIOi: 


B.l  UNCHIRPED  RECTANGULAR  PULSES 

In  Sec.  4.3.2,  a sketch  was  given  of  the  derivation  of  the  dispersion 
parameter,  A^  , for  a simple  rectangular  narrowband  pulse  of  bandwidth  B , 
and  with  the  implicit  assumption  that  the  radar  frequency  f is  always  much 
greater  than  the  plasma  frequency  everywhere  along  the  propagation  path. 

This  derivation,  and  the  results  to  be  discussed  below,  are  based  on  some 
previously-performed  but  as  yet  unpublished  work  by  one  of  the  authors. 


It  is  well  known  that  when  a radar  signal  passes  through  a dispersive 
medium  and  then  returns  to  be  passed  through  a receiver  filter  matched  to 
the  originally-transmitted  waveform,  two  major  effects  are  observed: 

1.  The  peak  signal-to-noise  ratio  at  the  output  of  the  filter 
is  reduced  by  a factor  L^  , the  dispersive  mismatch  loss. 

2.  Both  the  received  signal  and  the  ambiguity  function  which 
is  the  envelope  of  the  filter  output  are  distorted  with 
respect  to  the  undispersed  case,  the  distortion  causing  a 
stretching  of  the  function  in  the  range  direction  with  a 
consequent  loss  of  resolution. 

Both  of  these  effects  are  modeled  in  ROSCOE. 


For  an  unchirped  rectangular  pulse  in  the  absence  of  dispersion,  the 
matched-filter  output  is  triangular  in  the  range  dimension.  When  disper- 
sion is  present,  the  shape  of  this  cut  through  the  ambiguity  function 
changes,  becoming  longer,  lower,  and  more  rounded.  An  expression  for  this 
output  function  has  been  derived.  For  <<  1 , it  remains  very  similar 
to  the  undistorted  triangular  function,  only  having  Its  peak  amplitude 


| 


reduced  and  its  width  correspondingly  expanded  so  as  to  satisfy  energy 
conservation.  As  increases,  the  output  function  becomes  approximately 

Gaussian: 


X(R) 


4b* (R  - Rq) 

2 2 
❖ Ro 
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where  Rq  is  the  range  at  the  peak  of  the  return  and  6Rq  is  the  nominal 
range  resolution  for  the  pulse. 


By  requiring  that  energy  be  conserved  despite  the  presence  of  disper- 
sion, and  equating  the  areas  under  the  distorted  and  undlstorted  output 
functions,  the  parameter  b2  is  found  to  be  given  approximately  by 
b2  * 0.884  . 


If  the  dispersion-distorted  range  resulitrm  width,  6R^  , is  defined 
by  finding  where  |x(R)/x(^0)|  * 1/2  , it  turns  out  that  5R^  is  well  ap- 
proximated by 

a vr0 

There  only  remains  to  find  the  dispersive  loss  factor,  , ir.  terms 
of  which  the  ambiguity  function  (in  range)  and  the  range-resolution  width 
are  expressed.  This  quantity  is  given  as  a function  of  A^  by  a compli- 
cated expression  involving  Fresnel  integrals  and  trigonometric  functions. 
For  reasonably  small  values  of  A^  this  function  was  well  approximated  by 
a piecewise  polynomial  fit.  For  larger  values  of  A^  the  function  can  be 
expanded  into  a simple  "asymptotic"  form.  These  expressions  are  used  in 
subroutine  DISFERS  to  compute  L^CA^)  • 
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B.2  CHIRPED  PULSES 

Many  different  kinds  of  chirped  signals  are  used  by  radar  designers. 
Simple  chirped  rectangular  pulses  have  such  high  range  sldelobes  that  they 
are  hardly  ever  used  in  practical  radar  systems.  Instead,  some  form  of 
weighting  is  employed  to  reduce  these  sldelobe  levels  (at  the  expense  of 
slightly  broadening  the  peak  of  the  output  ambiguity  function  in  range). 
Rather  than  complicate  the  radar  model  by  inserting  a selection  of  possible 
weighted  waveforms,  it  was  decided  to  make  a single  reasonable  choice  to 
represent  the  "typical"  chirped  pulse,  defined  by  its  ambiguity  function 
in  the  range  direction. 


It  was  found  that  a waveform  lying  between  the  Hamming  and  the  Hanning 
weightings  (two  popular  types)  could  be  approximated  by  a malnlobe  of  the 
form 


x(x) 

out  to  its  first 

xOO 

for  |X|  > 6RX  . 

X i 


sinOrX/SR^ 

ttxM1 

null  (at  | X | - 

• constant 
Here, 


P 


iR^),  and 


R - R 


26R  /F 
o pc 


The  constant  range  sldelobe  level  is  input  by  the  user,  and  the  best  ap- 
proximation to  the  Hamming  and  Hanning  waveforms  was  achieved  when  P * 1,7  . 


This  is  the  ambiguity  function  (in  the  range  dimension)  which  is 
used  in  the  radar  model  to  represent  a chirped  signal  wherever  required. 


* fTTSW* 1 V- * w 


for  example  In  computing  the  results  of  dispersion,  interfering  targets, 
split-gate  range  tracking,  etc.  When  dispersion  is  present,  6R^  Is  re- 
placed by  6Rj  * . in  analogy  with  the  results  for  the  unchirped 

rectangular  pulse.  The  spreading  of  the  range  mainlobe  Is  thus  taken  into 
account,  and  the  assumption  Is  made  that  the  dispersive  loss  continues  to 
be  given  by  Lp(A^)  . However,  here  the  parameter  is  increased  by 

the  factor  F (the  pulse  compression  ratio)  compared  to  its  value  for 
the  unchlrped  pulse,  to  take  into  account  the  fact  that  the  effective 
bandwidth  is  now  B = . Since  the  ambiguity  function  is  normalized 

to  unity  at  X • 0 , the  sldelobe  level  is  increased  by  a factor  of 
to  compensate  for  the  actual  dispersive  loss  at  the  peak;  the  effective 
sldelobe  level  is  thus  increased  by  dispersion. 


This  approach  to  computing  dispersive  effects  for  chirped  waveforms 
is  largely  heuristic,  since  the  correct  expressions  have  not  yet  been 
evaluated  in  this  case.  The  results  appear  to  be  physically  reasonable. 


B.3  THE  DOPPLER  DIMENSION 

It  should  be  noted  that  the  ambiguity  functions  discussed  above  are 
all  in  the  range  dimension  only.  The  cross-range,  or  doppler,  coordinate 
is  taken  into  account  by  another  heuristic  procedure  which  seems  to  be 
reasonable  for  an  unchirped  rectangular  waveform  and  which,  for  chirped 
waveforms,  is  based  on  a substitution  which  gives  correct  results  in  the 
known  case  of  a simple,  frequency-modulated  rectangular  pulse. 

For  unchirped  pulses,  there  is  of  course  no  range-doppler  coupling 
to  produce  apparent  range  errors,  but  there  will  be  a doppler  mismatch  loss 
whenever  the  doppler  frequency  (f<j)  of  the  target  differs  from  the  frequency 
(f^)  for  which  the  doppler  filter  is  set.  This  mismatch  loss  is  given  by 

t ~ i sin  X 
uop  1 X 
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where 


X 5 — 6R  | f _ f | 

c o'  d d1 


Notice  that  dispersive  broadening  of  the  resolution  length  does  not  affect 
the  doppler  mismatch  loss. 


For  chirped  pulses,  the  doppler  mismatch  loss  is  given  by  a similar 
expression  except  that  the  exponent  -2  is  replaced  by  -2p  , where  p - 1.7 
for  the  approximate  "universal  chirped  waveform"  defined  above.  Actually, 
however,  it  is  not  necessary  to  compute  the  mismatch  loss  explicitly  for 
the  chirped  waveform,  since  this  is  automatically  taken  care  of  when  the 
apparent  range  to  the  target,  R*  , is  calculated  with  the  Inclusion  of  the 
effect  of  range-doppler  coupling.  In  the  absence  of  dispersion  this  is 
given  simply  by 

«R  6R- 

R1  • R'  - -f-®1  <?d  - V 


When  dispersion  is  present  (LD  > 1) , this  expression  becomes 

5R  6R. 

r:  ■ ri  - s -f—  <!a  - fd> 


which  is  justified  bv  a heuristic  argument.  Notice  that  by  using  Lp  - 1 
for  the  non-dispersive  situation,  both  cases  can  be  represented  by  the  same 
equation,  which  may  be  put  into  the  simpler  but  equivalent  form 


- R’  * -r  <fd  - £d> 
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APPENDIX  C 


MONOPULSE 

The  radar  monopulse  (angli -cracking)  system  is  modeled  in  ROSCOE  by 
three  subroutines,  each  one  pertaining  to  a different  situation. 

Since  the  monopulse  subroutines  operate  prior  to  the  introduction  of 
the  radar  measurement  errors  in  the  code  logic,  the  monopulse  output  for  a 
single  point  target  lying  within  the  antenna-pattern  mainlobe  is  the  target's 
correct  angular  position.*  This  is  the  output  produced  by  subroutine  MN0PLS1 
(Sec.  4.6.1).  However,  a monopulse  system  has  the  property  that  a single 
point  target  lying  outside  the  mainlobe  is  incorrectly  assigned  to  some 
position  within  that  lobe;  this  false  Indication  is  also  appropriately 
computed  by  MN0PLS1. 

When  there  are  several  point  targets  close  enough  to  interfere,  but 
only  a single  path  between  the  radar  and  each  target  (Np  ■ 1) , subroutine 
MNOPLS  is  called  to  compute  the  monopulse  angle  indication.  This  subroutine 
was  qualitatively  described  in  Sec.  4.6.2;  details  will  be  given  below. 

When  refraction  through  a striated  region  is  encountered  (Np  > 1) , 
subroutine  MNOPLSM  is  called  to  compute  the  angle  indication.  This  sub- 
routine was  described  qualitatively  in  Sec.  4.6.3;  details  will  be  given 
below. 


C.l  GENERAL  DEFINITIONS  AND  PROCEDURES 

Assume  a circular  mainlobe  centered  on  the  instantaneous  beam  axis 
(uQ,v0).  The  3-dB  one-way  beamwidth,  on  boreslght,  is  the  input  quantity 
66  . Because  the  beam  size  and  shape  change  in  angular  coordinates  as  the 
beam  departs  from  boreslght,  sine-space  coordinates  are  preferred.  The 


That  is,  correct  except  for  possible  errors  due  to  refraction. 


half-beamwidth  In  sine-space  is  defined  to  be 


6w  ■ sin(66/2) 

In  the  model,  the  angular  sldelobes  of  the  radar  beam  are  represented  by  an 
amplitude,  Ss  , the  one-way  power  level  relative  to  the  beam  axis  at  bore- 
sight  (Input  by  the  user),  which  either  remains  constant  or  is  tapered, 
depending  on  the  model  type. 

For  a given  (point)  target  with  sine-space  coordinates  (u,v),  the 
angle  off  the  beam  axis  In  sine-space  is  given  by 

2 2 1''2 
5 - [(u  - uor  + (v  - VQ)  j 

The  shape  of  the  one-way  malnlobe  of  the  sum  beam  is  represented  by 

the  relative  gain  factor,  F(£)  , defined  to  be  unity  on  the  beam  axis. 

In  the  present  circular-beam  model,  the  malnlobe  is  represented  by  a 
2 

I (sin  x)/x]  function  with  the  argument  adjusted  so  that  the  3-dB  point 
comes  at  the  right  place.  Within  most  of  the  malnlobe  this  is  a good  ap- 
proximation to  a wide  variety  of  real  antenna  patterns.  Thus  the  malnlobe 
relative  gain  function  is 

F(£)  - [-'p/*  ] 2 . for  |C|  < it/P 

where  P = 1.392/6w  . 

C.2  POINT  TARGETS 

Let  R'  be  the  indicated  range  of  the  target  currently  being  tracked, 
o 

and  let  6R£  be  the  effective  half-vidth  of  the  range  resolution  cell. 
MNOPLS  considers  only  targets  with  ranges  R^  that  lie  within  interfering 
range  of  R^  , i.e.,  such  that  |R|  - R^|  <_  6R|  . 

Taking  into  account  doppler  mismatch  and  the  waveform  ambiguity  func- 
tion, the  signal  amplitude  in  the  sum  channel  at  due  to  target  1 at 


172 


1 1 '!•!  HM*  1 


. 


R!  is  proportional  to 


F(Ci) 


‘ ^ ^ ^ • * 


where  Xjj  Is  the  (normalized)  waveform  ambiguity  function  and  is  the 

effective  cross  section  of  target  i . 


The  signal  phase  in  the  sum  channel  at  due  to  target  i at 

is,  apart  from  an  additive  constant 

h ‘ T [fo  + u«)*d(»]*l 


provided  that  the  target  lies  within  the  angular  mainlobe.  For  targets  in 
the  sidelobes,  an  extra  term  of  value  7T  should  be  added  whenever  the  target 
lies  in  a "negative"  sidelobe.  However,  the  model  ignores  this  extra  phase 
factor  on  the  grounds  that  such  targets  are  almost  always  greatly  dominated 
by  those  in  the  mainlobe. 


The  total  sum-channel  signal  voltage  (amplitude  and  phase)  at  R^ 


can  be  expressed  in  complex  notation  as 


-j$i 


• £ v 


In  the  difference  channels,  it  can  be  shown  that  for  a typical  kind 
of  monopulse  circuitry,  the  amplitude  due  to  the  ith  target  is  related  to 
that  in  the  sum  channel  by  an  expression 


A±(d) 


-jAi(s)  tan  2 (ui  - uq) 


for  the  u-plane  difference,  with  a similar  expression  for  the  v-plane. 
Thus  the  u-plane  dif 
can  be  expressed  as 


Thus  the  u-plane  difference  voltage  at  R^  due  to  all  contributing  targets 
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Vdu  “ -j£  Ai  tan  | <ui  - uo)e"j^ 
i 

and  similarly  in  the  v-plane. 

The  real  and  imaginary  parts  of  the  expressions  are  the  quadrature 
components  of  the  signals.  For  the  u-plane,  these  are 


SR  - E Ai  C0s  *i 
i 

si  - -Eai  sln  *i 

i 

Dr  3 “E  Ai  tan  I (ui  ■ Uo)  8in  *i 

i 

Dj.  = tan  j - uq)  cos  4> ± 

i 


It  can  be  shown  that  the  indicated  u-coordinate  of  the  "target"  in  the 
beam  is  given  by 


u 

o 


+ 


P 


where 


0(u)  = arc  tan 


2CVi  - V*] 


2 2 2 2 
SR  + SI  - °R  ' °I 


A similar  expression  applies  to  v'  , the  indicated  v-coordinate. 


17k 
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If  there  are.  In  fact,  no  targets  in  the  beam,  the  circuitry  will 
act  on  noise  alone  to  produce  a false  angular  position  estimate  located 
somewhere  at  random  within  the  mainlobe.  Subroutine  MNOPLS  does  this  when 
LISTA  contains  no  targets. 

If  there  is  only  a single  point  target  in  LISTA  which  can  contribute 
to  the  monopulse  determination,  the  above  expressions  become  so  much  sim- 
plified that  a special  subroutine,  MN0PLS1,  was  written  to  deal  with  this 
case.  If  this  target  point  lies  within  the  mainlobe  of  the  beam,  then  the 
Indicated  position  is  just  the  true  position  of  the  target,  i.e.,  u’  = u 
and  v*  = v • However,  if  this  target  point  lies  outside  the  mainlobe, 
the  system  will  produce  a false  (but  determinate,  provided  that  the  noise 
is  sufficiently  low)  angular  position  indication.  This  is  appropriately 
treated  in  subroutine  MN0PLS1. 

C.3  SPREAD  TARGET  CLOUDS 

For  spread  target  clouds  (Np  > 1) , the  monopulse  model  is  contained 
in  subroutine  MNOPLSM.  Its  basic  operations  are  the  same  as  those  used  in 
subroutine  MNOPLS  for  point  targets,  except  for  obvious  required  modifica- 
tions. For  example,  if  any  significant  portion  of  a target  cloud  (not  just 
its  center)  lies  within  the  range  resolution  cell,  that  portion  will  con- 
tribute to  the  combined  interference  signal  which  is  processed  by  the  mono- 
pulse system.  The  amplitudes  of  each  "target"  signal  are  computed  in  a 
manner  similar  to  that  used  for  point  targets,  except  that  the  target  cloud 
particle  density  distribution  must  be  convolved  with  the  waveform  ambigui.y 
function  to  compute  the  amplitude  contribution. 

Furthermore,  because  of  the  statistical  nature  of  the  "targets,"  the 
random  phase  approximation  is  used  in  place  of  adding  the  in-phase  and 
quadrature  components.  This  changes  the  final  equations  somewhat,  but  in 
an  obvious  manner.  The  relevant  equations  are: 
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f<L) 

^2  ^(r;  - Ri>^ 


where  is  the  effective  off-axis  position  for  the  given  cloud  (as  com- 

puted by  subroutine  MULTOAR)  and  b Is  a correction  factor  related  to  the 
convolution  of  cloud  statistics  with  the  ambiguity  function  (see  Appendix  E) , 


The  phase  is  now  irrelevant.  MNOPLSM  forms  the  power  sums 

* 3 IX 


Du  3 E Ai  c“2  \ (U1  - uo> 


with  a similar  expression  for  Dv  . 


Then  it  turns  out  that 


uo  + P 


where 


0(u)  = arc  tan 


i 


• Du  i 


and  a similar  expression  for  v'  . 


As  before,  appropriate  actions  are  taken  for  targets  which  are  too 
small  or  too  far  out  of  the  mainlobes  to  contribute  anything  but  essen- 
tially random  noise. 
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APPENDIX  D 


THE  SPLIT-GATE  RANGE  TRACKER 


A simple  split-gate  range  tracker  is  used  in  ROSCOE,  matched  to  the 
single  approximate  pulse-compression  waveform.  It  is,  of  course,  applicable 
only  to  those  events  for  which  a range  gate  is  defined.  There  are  two  sub- 
routines which  embody  the  split-gate  tracker:  subroutine  SPLTGAT  for  use 

when  Np  * 1 , and  subroutine  MLTSPLT  for  use  when  Np  > 1 . These  were 
described  in  a brief  qualitative  way  in  Secs.  4.6.4  and  4.6.5.  Further 
details  will  be  given  below. 


The  basic  principle  of  the  split-gate  range  tracker  is  to  integrate 
the  video  signal  which  is  the  output  of  the  matched  receiver  filter  over 
two  range  regions  which  together  cover  the  total  gate  width.  If  W is  the 
range-gate  width,  R£  its  center  range,  and  F (R)  the  video  output  en- 
velope as  a function  of  range,  then  in  effect  the  system  forms  the  integrals 


/ 


F (R)  dR 


R -H 
c 2 


R + H 
c + 2 


F(R)  dR 


Then  the  sum  and  difference: 


T = I + I 


e 2 I - I 


are  taken.  In  the  form  used  in  the  ROSCOE  model,  the  error  signal,  4 , 
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is  related  to  the  ratio  e/T  by  a known  function,  by  means  of  which  the 
"true"  range  of  the  target  relative  to  R£  is  determined. 

If  there  is  only  a single  target  in  the  range  gate,  and  if  there  is 

no  dispersive  distortion  of  the  ambiguity  function,  then  the  form  of  F(R) 

is  simply  | x (R  - Rt)  | • where  Rt  is  the  true  range  of  the  target. 

Assuming  the  waveform  to  be  known,  one  can  then  in  principle  find  Rt  as 

a function  of  e/T  not  only  for  e/6R^  <<  1 , as  is  the  usual  approxima- 
tion, but  in  fact  for  all  values  of  e . It  should  be  noted,  however,  that 
the  apparent  range  given  by  this  model  will  always  lie  within  the  range 
gate,  wherever  the  true  target  may  be;  in  other  words,  even  if  there  is  no 
target,  or  there  are  many  targets,  within  the  range  gate,  the  algorithm 
will  produce  a value  of  Rt  lying  between  Rc  - (W/2)  and  Rc  + (W/2)  . 
This  is  a desirable  effect  for  the  model,  however,  since  it  is  equivalent 
to  what  a real  split-gate  tracker  would  do  under  the  same  circumstances. 
Similarly,  if  the  waveform  has  been  distorted  by  dispersive  effects,  the 
appropriate  incorrect  range  will  be  found. 

In  addition  to  Rt  , the  split-gate  algorithm  also  outputs  a "peak 
signal  power  factor,"  FM  , which  is  proportional  to  the  amplitude  of  the 
target  signal  in  the  case  of  a single  target  and  non-dispersive  propagation. 
In  other  cases,  an  appropriate  erroneous  value  of  Fjj  Is  found. 


Assuming,  therefore,  a single  point  target  located  at  R^  , and  as- 
suming further  that  the  ambiguity  function  in  range  has  the  approximate 
form 

sinOrR/SR^  l1,7 
x ^ " itR/^ 

as  discussed  in  Appendix  B,  it  can  easily  be  shown  that 


c 


- 6R.S 
m 1 
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and 


T - 


where 


S(x) 


• / M 


1.7 


dt 


Since  the  function  S(x)  can  be  tabulated  numerically,  its  inverse, 
x = QCS)  can  also  be  tabulated.  Then  if  we  define 

s = S(n)e/T 


the  error  signal,  A , can  be  found  at  once  using 


6R, 


Q(s) 


Whence  the  target  raige  is  given  by 
Rt  - Rc  + A 

The  inverse  function  Q(s)  is  computed  by  function  QINV  (Sec.  4.6.10);  it 
should  be  noted  that 


sign(A)  ■ sign(e) 

and  that  if  |s|  > S(ir)  then  |a|  h 6Rl  . The  peak  power  factor,  FM  , 
is  given  by 


F, 


M 


- T/[fSRisWJ 


179 


r 


fW 


The  algorithm  used  In  subroutine  SPLTGAT  for  the  cases  where  there 
Is  only  a single  ray  path  between  the  radar  and  each  real  target  (Np  ■ 1) 
begins  by  adding  the  in-phase  and  quadrature  components  of  the  returns  from 
each  Interfering  target  at  each  point  within  the  range-gate,  weighted  ac- 
cording to  the  appropriate  value  of  x at  the  given  point,  integrating  to 
find  I+  and  I , and  then  proceeding  as  described  above. 

The  algorithm  used  in  subroutine  MLTSPLT  for  the  case  of  refraction 
through  a striated  region  (Np  > 1)  uses  the  statistical  properties  of  the 
"spread  target  cloud"  as  found  by  subroutine  REFRCT  to  approximate  the 
video  signal  envelope  which  is  input  to  the  range-tracker.  Because  of  the 
statistical  nature  of  this  model,  in-phase  and  quadrature  components  are 
not  computed  for  each  target  image;  instead,  the  random-phase  approximation 
is  used  in  comh'ning  the  contributions  of  different  images  to  the  video 
signal.  See  a.  the  discussion  in  Appendix  E. 
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APPENDIX  E 

THE  SPREAD  TARGET  CLOUD  IN  THE  RADAR  BEAM 


When  the  ray  path  between  a target  and  the  radar  Intersects  a striated 
debris  region,  numerous  multiple  paths  are  generated  (Np  > 1),  giving  rise 
to  the  so-called  "spread  target  cloud"  of  images.  In  ROSCOE,  these  Images 
are  described  statistically,  with  the  parameters 


AR,  Au,  Av,  oAR,  alf  a2>  p,  a,  a ^ 


as  defined  in  Appendix  G and  discussed  in  Appendix  B.  For  a point  target 
whose  true  radar  coordinates  are  (R,u,v),  the  center  of  its  apparent  "cloud" 
distribution  is  at 

6R* 

Ri  * R + 1R  + V "B®  (£d  - fd> 


u^  • u + Au 


V 1 ■ v + Av 


where  IpC  *1  if  the  radar  is  operating  in  a pulse-compression  mode,  and 

I ■ 0 if  not. 
pc 

This  appendix  describes  the  mathematical  modeling  of  the  interaction 
between  a spread  target  cloud  so  defined  and  a radar  beam. 


E.l  INTERACTION  WITH  THE  RADAR  ANGULAR  BEAM 

In  sine-space,  the  cloud  can  be  represented  by  an  error  elLipse 
centered  at  (u^,v^)  with  semi-major  axis  of  length  and  semi-minor 

axis  of  length  (these  define  the  one-sigma  error  ellipse).  Its 

orientation  is  defined  by  stating  the  angle,  p , between  the  major  axis 
of  the  ellipse  and  the  u-axis  of  the  sine-space  coordinate  system  (to 
eliminate  possible  ambiguities,  we  select  0 <_  p < tt). 


l8l 


The  mainlobe  of  a circular  radar  beam  centered  at  (u  ,v  ) can  be 

o o 

described  by  Its  one-way  relative  gain  factor 

»«>  ' 

as  described  in  Appendix  C,  Sec.  C.l.  Here  we  need  only  recall  that  P 
is  a known  constant  and 


(u  - u )2  + (v  - v )2 
0 0 


In  order  to  describe  the  relation  of  the  cloud  error  ellipse  to  the 
radar  beam  pattern,  it  Is  convenient  to  use  a rotated  set  of  sine-space 
coordinate  axes,  (u',v*),  which  have  been  rotated  about  (U0>VQ)  through 
the  angle  y , so  that  u'  is  parallel  to  the  major  axis  of  the  ellipse 
and  v’  is  parallel  to  the  minor  axis.  Then  the  new  and  old  coordinates 
of  any  point  are  related  by: 


u1  - u ■ + (u  - u ) cos  y + (v  - v ) sin  y 

o o o 


v'  -v  ■ - (u  - u ) sin  y + (v  - v ) cos  y 

o o o 


Using  this  rotated  coordinate  system,  the  effective  one-way  relative 
power  received  from  the  cloud  can  be  approximated  by 


00 


F(U 


/in  a. 


exp 


(u1  - u[)2 

2a? 


/2t To, 


exp 


(v*  - v')2 

2a? 


du'dv' 


provided  that  o^  and  a^  are  sufficiently  small  that  the  integrand  has 
negligible  value  outside  regions  where  sin  x * x . 

This  integral  can  be  solved  by  using  the  exponential  approximation 
to  the  function  F(£)  ; i.e., 
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sin  x -Ax 

a e 

X 

where  the  constant  A can  be  chosen  to  force  agreement  at  any  particular 
point.  If  the  agreement  la  to  occur  at  the  half -power  point,  then  the 
constant  Is  A a 0.1788  . Making  use  of  the  fact  that 

£2  = (u*  - u )2  + (v'  - v )2 

O 0 

the  double  Integral  takes  the  approximate  form 


00 

;1  * 2^7  ff  “»  - 2Ap2[<“'  - “</  + <v’  - vo>2j 


(u»  - up2 


2 o 


(V1  - vp2 
2o? 


du’dv’ 


This  Integral  can  be  solved  exactly.  The  result  Is 


F1  “ 


(Kj  + 1)(K2  + 1) 


-1/2 


F(0  exp 


K1  - 

4<u’ 

°1 


.2  . 2 v 2 

uo>  + ~2  (v Vo) 
°2 


I! 


where 


K 


2 _ 


4AP2o2 


.»  = 


K7u  + uJ 
1 o 1 

2 

K*  + 1 


r*  = 


K2Vo  ^ V1 
K2  + 1 
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and 


l2  = (G*  - u )2  + (v'  - V )2 

O 0 

The  point  (u',v')  [which  can  be  directly  transformed  into  the  un- 
primed coordinate  system  as  (u,v)]  represents  the  effective  center  of  the 
received  power  from  the  target  cloud  as  mediated  by  the  antenna  beam  pat- 
tern. In  the  same  way,  £ represents  the  angle  off  beam  axis  (in  sine- 
space)  of  the  effective  center  of  the  spread  target  cloud,  and  is  the 

effective  one-way  relative  gain  factor  of  the  cloud.  In  this  way,  we 
essentially  have  reduced  the  cloud  to  a fictitious  point.  The  quantities 
(R,u,v),  £,  and  are  computed  in  subroutine  MULTOAR  (here  R = R^). 


E.2  INTERACTION  WITH  THE  AMBIGUITY  FUNCTION  IN  RANGE 

2 

We  have  postulated  that  there  are  Np  target  Images  in  the  cloud, 
with  a Gaussian  range  distribution  centered  at  r • R'  and  having  standard 
deviation  oAR  (which  in  the  present  subsection  will  be  renamed  o^  for 

notatlonal  simplicity).  It  Is  assumed  that  a « R'  , and  of  course 

2 r 

>>  1 . The  number  distribution  of  the  target  images  in  the  cloud  is 


N! 


N(r) 


^2tTo 


exp 


(r  - R 'Y 

2 a2 
r 
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It  is  further  assumed  that  the  mean  cross  section  of  the  individual 

target  Images  in  the  cloud  is  a , with  standard  deviation  ag  . The 

total  cross  section  of  the  entire  cloud,  if  observed  as  a whole  at  great 

2- 

distance,  is  o^  ■ Jo  . 

As  before,  let  Xjjfr  “ ?l)  be  the  waveform  ambiguity  function  * 

(normalized  to  unity  at  r ■ r^)  giving  the  output  amplitude  at  r due 
to  a unit  target  at  r^  . Thus  the  amplitude  distribution  at  a point  r 
due  to  the  target  Images  located  around  point  r^  is  given  by  ' 
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P(r)  - V°XN(r  - r1)N(r1) 

The  total  amplitude  at  point  r la  found  in  principle  by  integrating 
P(r)  over  all  values  of  r^  , taking  into  account  the  appropriate  phases 
of  each  contribution.  Because  of  the  statistical  nature  of  the  model, 
however,  the  random  phase  approximation  may  be  used;  here  this  means  that 
we  add  the  squared  amplitudes,  ignoring  phase,  and  take  the  square  root  at 
the  end.  Thus 

CO 

[A(r)]2  “ o J [xN(r  “ i^)]2  N(rx)  d^ 

r*o> 


It  will  be  assumed  that  the  ambiguity  function  in  range  can  always 
be  approximated  (where  it  is  of  significant  magnitude)  by  an  expression  of 
the  form 

2 2 

XN(x)  * exp[-a  x ] 

where  a is  a known  constant,  depending  on  the  particular  waveform  involved. 
With  the  approximation  the  integral  can  easily  be  evaluated  to  yield 


I 


A(r)  “ VS Npb"1^2  exp[-a2(r  - r')2/b2] 


where 


b2  = 1 + 4a2o2 


The  peak  value  of  A(r)  occurs  at  r ■ R'  and  has  the  amplitude 
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This  may  be  interpreted  to  mean  that  the  effective  range  of  the  cloud 
center  is  at 


R = R' 


and  the  effective  cross  section  of  the  target  cloud  as  mediated  by  the 
ambiguity  function  of  the  signal  waveform  is 

°E  “ NP  5 *Vb 

Notice  that  the  effective  cross  section  of  the  target  cloud  as  seen  by  a 
reasonable  radar  waveform  is  always  less  than  its  ideal  total  cross  section. 

For  the  two  kinds  of  waveforms  used  in  ROSCOE,  the  value  of  the  para- 
meter b will  now  be  calculated. 


Simple  unchirped  rectangular  pulse.  The  ambiguity  function  (in  range) 
for  this  waveform  is  shown  in  Appendix  B to  be  approximated  by 

XN(x)  ■ exp[-  0.884  (2x/6R^)2] 


whence  it  is  seen  that 

2 

b2  - i + M.w(jfr) 

0 

Chirped  pulse.  The  range  ambiguity  function  for  the  single  chlrped- 
pulse  approximation  used  in  ROSCOE  was  defined  in  Appendix  B to  be 


X 


N 


St 


8in(irx/6Rp 

ttx/5R^ 


1.7 


Making  use  of  the  exponential  approxiation  to  the  X function,  the 

value  of  b^‘  can  be  computed  to  be 


b2  > 1 + 12.00(^7) 

E.3  MULTIPLE  CLOUDS  IN  RANGE  AND  ANGLE 

If  there  Is  more  than  one  spread  target  cloud  In  the  radar  range  gate, 
and  their  combined  interference  function  must  be  found  (as  for  example  in 
subroutine  MLTSPLT,  SLNINTM,  and  MNOPLSM) , the  following  extension  of  the 
above  theory  is  used. 

Suppose  there  are  n overlapping  target  cloud  distributions,  all 

having  the  same  value  of  Np  and  or  but  each  ahvlng  its  own  value  of 

R|  and  0^  . When  combining  these  distributions  with  respect  to  the  same 

radar  waveform,  we  may  use  the  random  phase  approximation  as  above,  but 

2- 

modified  so  that  each  individual  value  of  Np0^  is  replaced  by  the  appro- 
priately weighted  value 

Xi  2 NP°iFi/(RI)4  i 

-4 

Here  the. factor  (R*)  takes  into  account  the  different  ranges  of  the 

1 2 

cloud  centers,  and  (F^)  is  the  relative  gain  factor  (two-way,  this  time) 
to  account  for  the  different  positions  of  the  clouds  with  respect  to  the 
antenna  beam  pattern.  Replacing  the  factor  Np^"o  in  the  equation  for 
A(r)  by  r/xj  , and  noting  that  the  exponential  exp[-a2(r  - r')^/b2]  in 
that  equation  is  merely  an  approximation  to  the  ambiguity  function 
XNI(r  - R' )/b]  , the  combined  amplitude  at  point  r is  given  by 

2)1/2  ! 


* 
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APPENDIX  F 


STOCHASTIC  REFRACTION  MODEL 


F.l  THIN  SCREEN  APPROXIMATION 

In  simulations  of  radar  propagation  through  a disturbed  environ- 
ment it  is  desirable  to  be  able  to  determine  the  statistical  descriptors 
of  the  beam  as  it  passes  through  a medium  whose  properties  vary  with 
position.  The  medium  will  generally  have  a smoothly  varying  background 
electron  density  plus  an  additional  fluctuating  density  which  varies  much 
more  rapidly  with  position.  If  the  fluctuating  part  of  the  electron 
density  does  not  lead  to  excessive  departures  of  a ray  in  the  beam  from 
its  nominal  smooth  trajectory,  it  is  possible  to  develop  a theory  which 
describes,  in  statistical  terms,  the  effefcts  of  the  electron  density 
fluctuations. 

The  thin  screen  approximation  provides  a description  of  the 
propagation  of  electromagnetic  radiation  through  an  ionized  region  of 
space,  provided  the  following  three  conditions  are  met : 

1.  Ray  Optics.  The  electron  density  gradients  must  be  suffi- 
ciently small  that  the  fractional  change  in  electron  density 
over  distances  comparable  to  a wavelength  is  small.  In 
addition,  the  attenuation  over  distances  comparable  to  a 
wavelength  must  be  small. 

2.  Straight-Line  Approximation.  The  ray  deflections  must  be 
sufficiently  small  that  electron  densities  and  density  gra- 
dients along  the  true  ray  path  may  be  estimated  by  consider- 
ing a nominal  straight-line  ray  path. 


jtU»3Mrt-av“wl 


3. 


Small  Electron  Densities.  The  electron  densities  must  be 


sufficiently  small  to  permit  a series  expansion  in  powers 
of  the  electron  density. 

The  name  "thin  screen  approximation"  is  a slight  misnomer 
since  the  basic  assumptions  may  remain  valid  even  though  the  path  length 
through  the  disturbed  region  is  fairly  long.  However,  for  a given  angu- 
lar deflection  there  la  a limiting  path  length  beyond  which  the  straight- 
line  approximation  breaks  down. 

He  consider  a region  with  a magnetic  field  in  the  z-direction 

and  an  electron  density  n (x,y)  which  is  independent  of  z . In  the 

6 

absence  of  the  electrons,  a wave  packet  of  radiation  moves  In  a straight 
line  with  the  speed  of  light,  and  we  choose  our  coordinate  system  so  that 
the  nominal  ray  lies  in  the  x , z-plane,  making  an  angle  y with  the 
z-axls. 


The  electric  field  then  has  the  form 
E ■ Eq  exp(lk4  sin  y + lkz  cos  y - lwt) 

■ Eq  exp (14  - iwt) 

where  4 Is  a function  only  of  x and  y , satisfying  the  equation 

m2  - 1 ^-5—  -l-o. 

ncr  sin42  y 

according  to  the  ray  optics  approximation.  Here  n is  the  critical 

C t 

electron  density  at  the  frequency  f , and  k - 2ir/\  . This  equation 
determines  the  electric  field  at  an  arbitrary  point  in  space  if  the 
value  of  4 is  known  on  the  plane  x ■ xq  . If  the  electron  densities 
are  sufficiently  small,  4 nay  be  expanded  in  powers  of  oft  so  that 
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1 


'I'  “ <»  + <l<1  + ^2  + • * * 


On  the  surface  x ■ xq  , t|/(xQ,y)  ■ xq  + ^(xQ,y)  , where  ^ is  assumed  to 

be  a first-order  quantity.  We  therefore  take  <j>Q(x,y)  ■ x , and  this  des- 


cribes the  nominal  straight-line  path. 


The  the  first-order  corrections  to  are  determined  by 


27\p  *7^.  ■ -c 
ro  Y1  e 


so  that 


1 

- 2 °e 


^(x,y)  - iKxo,y)  f °e(x,»y)  dx' 


Similarly,  the  second-order  terms  give 


2 5/  - -«* 


where 


is  the  angle  that  the  projection  of  the  ray  path  onto  the  xty-plane  makes 
with  the  x-axis  (to  first  order).  Thus 


m,*+n** 


m 

*2(x.y)  - - £ J 


Oe(x',y)  dx’ 


Jl 


6 (x',y)  dx' 


On  the  basis  of  these  approximations  to  the  electric 
field  it  is  possible  to  evaluate  some  physically  interesting  quantities. 
Thus  the  z-component  of  the  velocity  of  a wave  packet  is  c cos  y , and 
the  velocity  components  in  the  x,y-plane  are  given  by  c sin  Y Vi/>  . Con- 
sidering only  the  motion  in  the  x,y-plane,  it  is  clear  that  a ray  takes 
a time  T to  traverse  a given  path  which  is  given  by 


cT 


sin  y 


dx 

|3i|>/8x| 


and  associated  with  this  time  Interval  there  is  an  apparent  range  given 
by  R - cT  . It  is  convenient  to  roughly  separate  the  contributions  to 
the  apparent  range  into  a part  PL  , which  is  the  path  length  along  the 
ray  trajectory,  and  to  associate  any  additional  contributions  with  the 
group  velocity  of  the  wave  packet.  Thus 


R sin  y * cT  sin  y 


PL  + cT,  , sin  y 
del 


I 


4 


with 


PL-y  d>  i7*i 


and 

cTd.i sln  y -/wr -fil 

1/2 

Since  3iJ»/3x  is  known,  and  |v^|  - (1  - og)  , it  is  a simple  matter 
to  expand  these  expressions  in  powers  of  0£  . In  evaluating  the 
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expressions  obtained,  it  oust  be  remembered  that  y varies  with  x along 
a ray  trajectory  so  that 

. MlpL  s 8 

dx  h7*x 

The  results  of  such  an  expansion  are 


"•  • * - *0  + rdl 


where  the  "dogleg"  range  error  is  given  by 


i f dx’ 


cTd.l  ,tB  Y 


A 

"7  f °c<x'»y,) dx' +i  f dx' 


so  that  the  apparent  range  for  a ray  going  from  xq  to  x is 


R - x - x + -3 r,,  + cT.  , 

o sin  y dl  del 


In  applications,  the  second-order  contribution  to  ^ Is 
generally  negligible  compared  to  the  first-order  term,  and  may  be  ignored. 
However,  because  of  geometric  factors,  the  dogleg  term  r^  can  be  ap- 
preciable, and  should  be  retained. 


1 


The  phase  of  a ray  moving  through  the  medium  changes  so  that 


■ k sin  y + h cos  Y 
• k sin  y | Viji  | dl  + k cos2  Y 6T 

since  in  the  x,y-plane  6r  - V^6t/|Viji|  . It  Is  then  easily  seen  that  the 
phase  change  in  moving  from  xq  to  x is  given  by 

A*  - kR  - 2kcTdel  sin2  y 


In  deriving  the  last  equation,  an  approximation  has  been  made  in  that  some 

liberties  were  taken  with  the  second-order  contributions  to  T . , . 

del 

F.  2 STOCHASTIC  REFRACTION  MODEL 

(U)  In  the  last  section  an  approximate  expression  for  the  electric 
field  was  obtained  in  terms  of  the  electron  density  distribution.  In 
practice  the  electron  density  function  is  often  not  known  in  detail  but 
must  be  described  in  statistical  terms.  We  therefore  consider  a model  in 
which  the  electron  density  fluctuations  are  due  to  field  aligned  stria- 
tions  so  that 


ne(x,y) 


n 


cr 


v(x)  + ^vk  exp 
k 


2 


(x  - x^2  + (y  - yk)2 


where  the  location  of  the  strlation  axes  (xk»yk)  * the  peak  electron  den- 
sity of  the  strlation  v.  n and  the  strlation  size  a.  are  not  known 
J k cr  k 

precisely.  Instead  we  assume  are  randomly  distributed  in  the  x,y- 

plane  with  a two-dimensional  density  function  p(x)  , that  v.  and  a^ 

/ ni  n2\ 

are  also  randomly  distributed,  and  that  the  average  moments  \v  a / 
are  known  as  functions  of  x . It  is  then  desired  to  calculate  the  expected 
value  of  various  physically  interesting  quantities;  this  is  done  in  the 
following  subsections. 
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F.2.1  Integrated  Electron  Density 
The  quantity 


f n (x'.y) 
d(x,y)  - J dx'  — - 


cr 


Is  linearly  related  to  the  absorption  due  to  electron-neutral  collisions 
and  is  also  useful  in  that  several  other  quantities  are  readily  calcu- 
lated once  d is  known.  When  a ray  passes  the  kth  striatlon  there  is  a 
contribution  to  d(x,y)  given  by* 


•0 

dk(y)  - vk  exp|-(y  - yk)2/\]  / «*p(-x2/\)  to 


1/2  , .2 , 
vkakn  exp  -(y  - yfc) 


[-(y  - yk)2/l^] 


and 


d(x,y) 


•/ 


v(x')  dx*  + 


where  the  sura  over  j extends  over  all  strlations  in  the  region 
xQ  < Xj  < x and  -•  < y^  < «®. 

The  average  contribution  of  those  strlations  in  the  interval 
from  x'  to  x*  + dx  can  be  obtained  by  replacing  the  sum  over  j by 
an  Integral  over  the  possible  positions  of  the  striatlon  centers  so  that 

(U)  The  use  of  the  limits  -«  to  +»  for  the  x-integration  is  justi- 
fied by  the  assumption  that  a^  is  small  compared  to  other  scale  sizes 
of  the  problem. 
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dy 


4 


• / 


p(x')  dx 


va(w)^2e”y 


where  the  expectation  value  under  the  integral  sign  is  an  average 
over  the  distribution  in  v and  a . Thus 

x x 

D(x)  - ^d(x,y)^  - / v(x')  dx'  + ir  J p(x')^va2^  dx' 

X_  X 

o o 

F.2.2  Fluctuations  in  Integrated  Electron  Density 
Consider 


d(x,y) 


S dt(y)  Ya  dj(y> 

i j 


As  before  the  average  value  of  this  quantity  nay  be  obtained  by  replacing 
the  suns  by  Integrals.  However,  the  terns  with  i - j must  be  considered 
separately  since  they  also  give  rise  to  a finite  contribution.  Thus 


t 


1 


I 
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where  the  second  term  on  the  right-hand  side  comes  from  the 
terms  with  i ■ j . Thus  the  fluctuations  in  d(x,y)  are  given  by 


Op00  »/|d(x,y)  - D(x)j  J*  p(x')^2a^dx' 


F.2.3  Emergence  Angle 

Each  strlation  scatters  the  ray  through  an  angle 


vk*1/2(y  - yk)  f 2 

0k ~2 (y  ’ yk>  ' 

^ *ln  Y • 


so  that  the  ray  makes  an  angle  with  the  x-direction  given  by 


0(x)  - EfiJ 


which  is  easily  seen  to  have  zero  expectation  value.  The  square  of  this 
angle  has  expectation  value 

0*<x>  -<e2(,).<EE  (ej)2) 


and  replacing  the  sum  by  an  Integral  gives  the  result 


o2(x) 
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with 


1 

2 sin*  y 


F.2.4  Miss  Distance 

A ray  will  suffer  a displacement  In  the  y-directlon  given  by 
Ay(x)  - Z (x  - xk)e® 

which  has  zero  mean.  The  expected  value  of  Its  square  Is  readily  found 
and 


(x) 


(x  - x') 


2 def 

dx' 


dx' 


F.2.5  Range  Error 

A ray  in  propagating  from  the  kth  to  the  k + 1th  strlatlon 
makes  an  angle  6^  with  the  x-axls  given  by 


k 


e 


k 


£ 

J-l 


e 


8 

J 


and  the  distance  It  must  travel  exceeds  the 

2 

by  an  amount  (1/2)  (x^  - x^je^  . Thus  in 
Is  a "dogleg"  range  error  given  by 


nominal  distance  x^+^  - x^ 
traversing  the  medium  there 


rdl  " 2 ^ <*fcfl  " 
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which  has  expectation  value 


A A j 

Rdi  ■ 7 / °e<x')  *■'  • 7 / (x  - x,)  ST  dx' 

x. 

The  total  range  error  consists  of  this  term  plus  a tine  delay  term 

R - ~ — R..  + cT . . 
e sin  y dl  del 


where 


cT 


M 2.1n3y 


F.2.6  Range  Fluctuations 
Consider  first 


rd!  ■ { £ <Vl  ' x»)em  £ <Vl  - xn)6n 


It  is  easily  seen  that  the  terms  with  m • n give  a vanishing  contribu- 

2 

tlon  to  the  expectation  value  of  r^  so  that  it  is  sufficient  to  con- 
sider n.  > m and  multiply  by  2.  For  n > m , 

n n m m 


■ o, 


IEEE  6'e;eKe’t 

i-1  j-1  k-1  i-l  / 

+ 2[O0(X»)]2  + 2*<x„> 
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where  Che  first  tern  cones  from  1 ■ J,  k ■ t , end  the  next  term  comes 
from  1 ■ k,  J ■ t , and  1 ■ i,  J ■ k contributions,  and  i ■ J ■ k ■ 1 
gives  rise  to 

*<X)  ' <6‘)4)'  l(i)5/2  f dx> 


y Cx  - x')|[og(x')]  2 + d-Cx’)j  *«' 


The  cross  correlation  between  the  dogleg  and  line  delay  terms  is  given  by 


v< 


ctdelrdl)“  cTdelRdl 


3_  2/1) 
16  * \3J 


— V-  f (x 

sin8  y 


- x’)p(x') 


'«  * 7"7T“  Bt  + JI»  ‘ " , 

e 4 sin  y e 
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F.2.7  Beam  De focussing 

Consider  two  rays  which  initially  make  angles  with  the  x-axis 
which  differ  by  A0o  . After  propagating  through  the  medium,  two  such 
rays  will  be  separated  by  a distance 

6y(x)  - xAflo  + - xk) 

k 7 


The  sum  over  k has  zero  expectation  so 
(6y)2  . (Note  that  6y<3t^)  and  38*/3y 


If  we  let  S(x)  - then 


that  it  is  necessary  to  consider 
xr e statistically  Independent.) 


S(x)  - 


)S(x’)  dx' 


However,  because  of  the  presence  of  S on  the  right-hand  side  the  final 

result  for  S Involves  an  expansion  containing  all  powers  of  the  peak 

electron  density.  Since  in  the  underlying  theory  all  powers  of  the  peak 

electron  density  higher  than  the  second  have  been  Ignored,  these  higher 

prder  sums  in  S are  not  reliable.  Therefore,  in  applications  it  is 

2 

prudent  to  replace  S by  its  "Born  approximation"  S(x)  ■ x when  it 
occurs  under  the  Integral  sign. 

The  gain  G due  to  focussing  is  then 
C(x)  - [x2/S(x)]1/2 

It  must  be  emphasized  at  this  point  that  the  gain  (always  less 
than  unity)  given  by  this  equation  la  associated  with  a single  ray  path 
connecting  the  radar  to  the  target.  In  general,  there  may  be  many  such 
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one-way  images  and  the  total  energy  on  the  target  contains  a factor  NjG  . 
If  the  radar  beamwldth  is  sufficiently  large,  it  can  be  argued  that  N^G 
is  on  the  average  equal  to  1.  The  energy  in  the  beam  sweeps  through  any 
region  of  space  on  the  average  times  at  the  range  considered.  If  it 
is  assumed  the  electric  fields  associated  with  different  rays  are  inco- 
herent, then  the  energies  associated  with  different  rays  add  and  energy 
conservation  demands  that  N^G  * 1 . The  angular  separation  between  rays 
at  the  radar  associated  with  successive  images  can  also  be  estimated. 

As  before,  denote  the  miss  distance  by  Ay  . During  an  angular  Interval 
such  that  ^ > 0 a ray  will  pass  through  an  Interval  of  angle  such  that 
Ay  starts  at  a negative  minimum,  hits  the  target,  and  then  rises  to  a 
local  maximum.  If  the  root  mean  squared  angle  traversed  in  going  from  a 
minimum  to  a maximum  is  A8^  then 


Thus  if  the  antenna  beamwldth  is  y then  the  number  of  Images  is 


N 


1 


■jr.x. . . 

2o.  G 


and  the  energy  incident  on  the  target  is  multiplied  by  a factor 
(yx)/(2o  ). 

The  treatment  developed  above  has  some  advantages  over  previous 
formulations.  In  the  past,  expressions  for  the  dogleg  range  error,  its 
scintillation,  and  the  focussing  gain  have  been  obtained  by  heuristic 
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arguments.  The  formulas  for  this  paper  have  been  derived  more  systematically 
and  to  the  extent  that  the  underlying  thin  screen  assumptions  are  valid  they 
* are  correct.  If  there  are  several  targets  in  a radar  beam  it  is  possible, 

using  the  present  formulation,  to  integrate  up  to  the  first  target  and  then 
simply  continue  on  to  the  next  target.  In  the  older  formulations  it  was 
l necessary  to  redo  the  entire  calculation  for  each  new  target. 
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DEFINITIONS  OF  SYMBOLS 


Underlined  terms  are  lists  In  the  sense  explained  in  Vol.  1.  An 
asterisk  (*)  denotes  an  Input  parameter.  A double  asterisk  (**)  denotes 
a precomputed  quantity  derived  from  input  parameters  stored  in  some  dataset. 


Vi 


Dispersion  constant  u**ed  In  subroutine  REF1S, 

A_  - 0.0828  B f "3/2 
Dn  n o 


A 


x 


* 

B 


** 

B 

n 


B 

v 


C 


Dispersion  parameter  defined  In  REF1S 

Constants  relating  to  the  standard  deviations  of  the 
radar  measurement  errors ; these  give  the  signal-to- 
noise-independent  term  in  the  respective  expressions 

Multiplying  factor  used  in  POSSV,  POSLIST  to  Include 
possible  interfering  targets  within  a specified  angle 
of  the  desired  target 

Effective  signal  bandwidth 

Noise  bandwidth,  * 1/ (pulse  durat  ion) 

Constants  relating  to  the  standard  deviations  of  the 
radar  measurement  errors;  these  are  used  to  find  the 
signal-to-nolse-dependent  term  in  the  respective  ex- 
pressions 

Total  received  clutter  power  from  fireballs  and  dust 
clouds 


u 


I 


Dfi  Slope  of  the  angle  sidelobe  taper  (not  currently 

used) 

* 

{ ,fQ  Symbols  used  Interchangeably  for  the  radar  transmitted 

center  frequency 

F,F^  Symbols  used  Interchangeably  for  the  one-way  relative 

gain  due  to  off-axis  position  of  the  target  In  the  radar 
beam 


F^  Pulse  compression  ratio  for  the  given  waveform 

Expected  doppler  shift,  to  whl-h  the  receiver  filter 
is  matched 


fd  Actual  doppler  shift  for  a given  target 


FM.FMAX  Symbols  used  interchangeably  for  the  amplitude  factor 

from  which  the  output  signal  power  can  be  computed 


Altitude  below  which  It  Is  assumed  that  no  refraction 
takes  place  (used  In  REF1S) 


Hrad 


Current  altitude  of  the  radar  (used  in  REF1S) 


1 


E 


2 

Integrated  electron  density  (electrons/cm  ) along  the 
path  to  the  target  (computed  in  REFRCT) 


IFAIL  Integer  index  of  success  or  failure  of  threshold  cros- 

sing (see  RADIS,  etc.) 

IFDET  Integer  index  output  by  RADMODS;  IFDET  - 0.  indicates 

that  no  detection  has  occurred 
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ITO  Serial  number  in  POSLST  of  the  target  of  interest 

(computed  in  RADIS) 

Lp  Dispersive  loss  factor  (computed  in  D1SPERS) 

Lp  Faraday  rotation  loss  factor  (computed  in  REF1S) 

Two-way  path  absorption  to  the  target  point  (computed 
in  ABSORB) 


L 

om 

LISTA, LISTB 
MLST 


A maximum  absorption  parameter  (computed  in  RADIS) 

Lists  of:  interfering  targets  (see  Sec.  4.4.5) 

Number  of  targets  in  LISTA  or  LISTB  (computed  in  REFLST1 
or  REFLSTN) 


MGRP  List  of  targets  in  mutually  unresolvable  groups  (com- 

puted in  REFLST1) 

NLST  Number  of  targets  in  list  POSLST 

NGRPS  Number  of  close-target  groups  in  LISTA  (see  REFLST1) 

**N  Nominal  system  noise  power,  N ■ kB  T 

n 1 ' ’ n n n 

Np  Number  of  direct  propagation  paths  between  radar  and 

target  (computed  in  REFRCT) 

NTq  Identification  number  of  the  target  of  current  interest 

(used  in  POSLIST  and  POSSV) 

p Beamwidth  factor  (see  GOA,  etc.),  p ■ 1.392/6w 
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Possible-targets  list;  see  Sec.  4.2.3 


2 

Range  at  which  aim  target  will  just  exceed  threshold 
Range  to  the  center  of  the  receive  range  gate 
Range  to  a given  target  point 

t 

Range  to  the  particular  target  point  defined  in  RADIS 

Range  to  the  apparent  center  of  the  "spread  target 
cloud"  in  the  striated  case  (see  MULTOAR) 

Apparent  target  range  (an  ultimate  output) 

Apparent  target  range  when  used  as  an  input  to  the 
monopulse  routines  MNOPLS  and  MNOPLSM 

Range  to  the  highest  output  peak  in  SVPEAKS 

Additive  range  used  in  POSLIST  to  Include  possible 
interfering  targets  within  3000  feet  of  the  range 
gate 


Parameter  used  in  computing  received  signal  power, 
So  - <s/N>IiR08n 


Revised  parameter  used  in  computing  received  signal 
power;  after  REF1S  it  has  the  value 
2 


< 


I 


I 


f 


1 


< 


S One-way  angular  sidelobe  (power)  level,  relative  to 

s 

unity  on  the  bean  axis  at  borealght  (the  near-in  value 
for  tapered  aldelobe  model) 

*SgL  Constant  range  sidelobe  (power)  level,  relative  to 

unity  at  the  peak  of  the  undispersed  ambiguity  function, 
for  chirped  pulses  (used  in  AMBGN) 

(S/N),^  Pre-detection  signal-to-nolse  threshold  for  the  given 

event  and  waveform 


** 


(u,v) 


<U  ,V  ) 
o o 


(vV 


(u,v) 


Parameter  used  in  split-gate  range  tracking  (see 
SPLTGAT , MLTSPLT , and  QINV) 

Time  of  the  current  event 

Effective  receiver  temperature 

System  reference  temperature,  Tq  ■ 290*K 

Nominal  system  noise  temperature,  - Tq  + Tg 

Total  effective  outside  noise  temperature 

Sine-space  coordinates  of  a given  target  point 

Sine-space  coordinates  of  the  current  radar- beam  axis 

Sine-space  coordinates  of  the  particular  target  point 
defined  in  RADIS 

Sine-space  coordinates  of  the  apparent  center  of  the 
"spread  target  cloud"  in  the  striated  case  (see  MULTOAR) 
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(u'.v') 


W 

*T 

y 


z 

o 


a 

AR 

(AA.AE) 


(Au, Av) 


AXy 


** 


6R 

o 


Apparent  target  sine-space  coordinates  (an  ultimate 
output) 

Current  value  of  the  width  of  the  receive  gate 
Total  noise  plus  clutter  power 

Trial  signal-to-noise  ratio  which  Is  compared  with 
to  see  if  a threshold  crossing  is  possible,  in 
RADIS  and  REF1S 

Signal-to-noise-plus-clutter  ratio  at  the  output  of 
the  receiver  filter  (an  ultimate  output) 

Safety  factor  used  In  RADIS  (see  Sec.  4.3.1) 

Refractive  range  bias  error  (see  Appendix  B) 

* 

Refractive  angular  bias  errors  when  Np  ■ 1 (see 
Appendix  B) 

Refractive  sine-space  bias  errors  to  center  of  "spread 
target  cloud"  when  > 1 (see  Appendix  B) 

Equivalent  noise  Increment  due  to  range  sidelobe  target 
Interference 

Nominal  range-resolution  width  for  unchirped  pulse, 

6R  « ct/2 
o 

Radar  bias  error  in  the  range  coordinate 


210 


Nominal  range-resolution  width  for  a chirped  pulse, 

5R.  * 2oR  /F 
1 o pc 


Dispersion-distorted  range  resolution  widths, 
6R;  = Ld6Ro  and  6RJ  * 


Radar  bias  errors  in  the  u,v  coordinates 


One-way,  3-dB  half-beamwidth  in  sine-space, 

6w  ■ sin(66/2) 

One-way,  3-dB  angular  beamwidth  at  boreslght 

Angle  off  boresight  of  radar  beam  axis  (computed  in 
RADIS),  cos  0 ■ (l  - u2  - v2)1/2 

Slope-angle  in  sine-space  of  the  "spread  target  error 
ellipse,"  relative  to  the  u-axls  (see  Appendix  B) 

Angle  off  beam  axis  in  sine-space  for  a given  target, 

r r . , v2-,l/2 

C - [u  - uQ)  + (v  - vq)  J 

Angle  off  beam  axis  in  sine-space  of  the  apparent 
center  of  the  "spread  target  cloud"  in  the  striated 
case  (see  MULTOAR) 

Standard  deviation  of  the  refractive  range  bias  error 

Standard  deviations  of  the  refractive  angular  bias 
errors  when  ■ 1 (see  Appendix  B) 

Semi-major  and  semi-minor  axes,  respectively,  of  the 
sine-space  error  ellipse  in  the  "spread  target"  case 
(see  Appendix  B) 


o Mean  cross  section  of  the  Individual  targets  in  the 

"spread  target  cloud”  for  the  striated  case 

00  Standard  deviation  of  the  distribution  (about  a)  of 

the  individual  "spread  target"  cross  sections 

Xjg  Waveform  ambiguity  function  normalised  to  unity  at  its 

peak  (see  AMBGN) . For  chirped  pulses , its  argument 
already  contains  an  adjustment  due  to  range-doppler 
coupling  (see  REFLST1  and  MULT OAR) 

ftp  Faraday  rotation  angle  (vhera  applicable;  see  REF1S) 

t Pulse  duration 
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OTHER  GOVERNMENT  AGENCIE8 


Department  of  Commerce 
Office  of  Telecommunications 
Institute  for  Telecom  Science 
ATTN:  William  F.  Utlaut 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 


Aerospace  Corporation 

ATTN:  Norman  D.  Stockwell 

Brown  Engineering  Company,  Inc. 
ATTN:  Joel  D.  Blgtey 

Calapan  Corporation 

ATTN:  Romeo  A.  Dellberis 


ESL,  Inc. 

ATTN:  R.  K.  Stevens 
ATTN:  James  Marshall 
ATTN:  V.  L.  Mower 

General  Electric  Company 
TEMPO-Center  for  Advanced  Studies 
ATTN:  Tim  Stephens 
ATTN:  Warren  S.  Knapp 
ATTN:  DASLAC 

General  Research  Corporation 
ATTN:  John  Ise,  Jr. 

ATTN:  John  Boys 
ATTN:  Joe  Gabarlng 

Johns  Hopkins  University 
Applied  Physics  Laboratory 

ATTN:  Document  Librarian 

Lockheed  Missiles  4i  Space  Co.,  Inc. 
ATTN:  D.  R.  Churchill 

M.I.T.  Lincoln  Laboratory 

ATTN:  Lib  A-082  for  David  M.  Towle 

Martin  Marietta  Aerospace 
Orlando  Division 

ATTN:  Roy  W.  Heffner 

Mission  Research  Corporation 
ATTN:  Dave  Sowle 
ATTN:  W.  F.  Crevier 
ATTN:  Paul  Fisher 
ATTN:  Roy  Hendrick 
ATTN:  Russell  Christian 

Physical  Dynamics,  Inc. 

ATTN  : Joseph  B.  Workman 

Science  Applications,  Inc. 

ATTN:  D.  Sachs 
ATTN:  Raymond  C.  Lee 
ATTN:  Curtis  A.  Smith 

Science  Applications,  Inc. 

HuntsvU'e  Di  vision 

ATTN:  Noel  R.  Byrn 
ATTS'v:  JaleDivia 

Stanford  St  wjarch  Institute 
ATTN:  Jacqueline  Owen 
ATTN:  Walter  G.  Chesmit 

Visldyne,  Inc. 

ATTN-  J.  W.  Carpenter 
ATTN:  Henry  J.  Smith 
ATTN:  Chuck  H.  Humphrey 
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